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Abstract 

The  influence  of  stimulated  scattering  on  laser  intensity  in  fiber  optic  waveguides  is 
examined.  Stimulated  Brillouin  scattering  (SBS)  in  long,  multimode  optical  waveguides 
is  found  to  generate  a  Stokes  beam  that  propagates  in  the  fiber  LPoi  mode.  This 
characteristic  of  the  Stokes  beam  was  first  applied  to  beam  cleanup,  where  an  aberrated 
pump  generated  a  Gaussian-like  Stokes  beam. 

Additionally,  the  same  process  is  found  to  combine  multiple  laser  beams  into  a  single 
spatially  coherent  source.  The  mean  square  difference  between  the  two  beams  was  used 
to  measure  the  degree  of  spatial  overlap,  demonstrating  spatial  coherence  between  the 
Stokes  beams  even  when  the  pump  beams  are  not  spatially  correlated.  This  result  is 
obtained  regardless  of  whether  the  pump  beams  are  at  the  same  or  different  frequencies, 
producing  two  temporally  coherent  or  incoherent  Stokes  beams  respectively. 

Limitations  in  beam  cleanup  and  combining  are  also  examined  to  identify  ways  to 
overcome  them.  Output  couplers  are  designed  that  could  be  used  to  spatially  filter  the 
Stokes  beam  from  the  pump,  thus  increasing  the  number  of  beams  that  could  be 
combined.  The  combined  power  restriction  induced  by  second  order  Stokes  threshold  is 
examined  experimentally  and  theoretically  and  is  not  found  to  be  a  significant  limitation. 
Finally,  stimulated  Raman  scattering  (SRS)  beam  cleanup,  which  would  relax  the 
stringent  spectral  requirements  on  the  pump  beams,  was  also  investigated. 

The  last  portion  of  the  dissertation  theoretically  examines  suppression  of  stimulated 
Raman  scattering  in  fibers  to  eliminate  the  restriction  this  imposes  on  the  power  of  a  fiber 
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laser  or  amplifier.  The  suppression  was  modeled  using  both  a  holmium  dopant  and 
adding  a  long  period  grating  to  the  fiber.  Both  methods  were  shown  to  have  a  signifieant 
effect  on  the  SRS  threshold. 
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LASER  INTENSITY  SCALING  THROUGH  STIMULATED  SCATTERING 


IN  OPTICAL  FIBERS 


1.  Introduction 


1.1,  Motivation 

The  department  of  defense  has  identified  laser  systems  as  eritieal  to  the  modern 
military  [1],  They  are  used  for  numerous  applieations  ineluding  deteetion,  targeting, 
laser  radar,  remote  sensing,  munitions  guidanee,  and  even  as  weapons  in  their  own  right. 
To  be  effeetive  most  of  these  applieations  require  high  intensity  illumination  of  the  target. 
In  laser  radar,  for  example,  a  laser  pulse  is  sent  from  the  sensor  to  a  target.  The  sensor 
then  images  the  seattered  light  and  measures  the  time  of  flight  of  the  light  pulse  to  eolleet 
a  3-D  image  of  the  target.  Sinee  an  aberrated  laser  beam  diverges  faster  than  a  Gaussian 
beam,  any  laser  aberrations  reduee  image  resolution  and  lower  the  intensity  on  target, 
limiting  the  sensor  operational  range.  For  these  reasons,  great  eare  is  taken  to  ensure 
high  quality  laser  beams  are  used.  Unfortunately,  high  power  solid  lasers  tend  to  be 
highly  aberrated  [2],  thus  novel  eoneepts  must  be  developed  to  improve  the  beam  quality 
of  high  power  systems. 

Stimulated  Brillouin  seattering  (SBS)  in  long  optieal  fibers  is  one  way  to  generate  a 
high  quality  laser  beam.  This  method  has  many  benefits  ineluding  a  low  threshold,  a 
high  operational  dynamie  range,  its  simplieity,  and  its  sealability  to  multiple  laser  beams. 
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A  low  threshold  is  critical.  The  threshold  represents  light  that  is  lost  and  the  higher  the 
threshold,  the  lower  the  overall  efficiency  of  the  system.  SBS  is  typically  one  of  the  first 
nonlinear  processes  to  appear  in  fibers.  Silica  fibers  do  not  have  a  high  third  order 
nonlinear  coefficient,  but  the  guiding  nature  of  the  fiber  allows  for  extremely  long 
interaction  lengths  resulting  in  very  low  thresholds  [3]. 

Several  factors  can  limit  the  dynamic  range  of  an  SBS  system.  One  of  these  is  the 
second  order  threshold.  This  is  the  power  level  at  which  the  beam  generated  through 
SBS  (Stokes  beam)  is  strong  enough  to  excite  a  second  Stokes  beam.  The  reasons  stated 
above  may  lead  one  to  deduce  that  the  second  order  threshold  will  be  very  low. 
Fortunately,  however,  the  first  order  Stokes  beam  grows  nearly  exponentially  in  the  fiber. 
Thus,  there  is  only  a  strong  interaction  near  the  front  end  of  the  fiber,  resulting  in  a  short 
interaction  length  and  a  high  second  order  threshold.  For  this  reason,  the  second  order 
Stokes  threshold  is  not  a  significant  limiting  factor  in  the  generation  of  a  high-power, 
high-quality  beam. 

The  simplicity  with  which  SBS  can  generate  a  high-quality  beam  is  one  of  its  greatest 
advantages.  It  only  requires  that  the  aberrated  beam  be  coupled  into  a  long  optical  fiber. 
The  resulting  Stokes  beam  will  be  Gaussian-like  regardless  of  any  fluctuations  in  the 
transverse  mode  distribution  of  the  source.  Methods  to  produce  a  high  quality  beam 
through  phase  conjugation  can  be  more  difficult  to  implement  since  phase  conjugation 
merely  allows  one  to  recreate  the  source  transverse  profile  after  traveling  through  an 
aberrating  medium  [4-8].  If  the  source  is  not  diffraction  limited,  the  phase  conjugate 
beam  will  not  be  diffraction  limited. 
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Finally,  the  generation  of  a  Gaussian-like  beam  from  an  aberrated  souree  through 
SBS  is  easily  expandable  to  eombine  multiple  lasers.  Multiple  sourees  at  the  same 
frequeney  merely  appear  as  a  highly  aberrated  beam  eoupled  into  the  fiber  and  eombine 
to  form  a  single,  high-power,  Gaussian  beam  [9].  Multiple  sourees  at  different 
frequeneies  independently  generate  Stokes  beams,  eaeh  having  nearly  identieal  Gaussian 
profdes. 

Fiber  lasers  offer  a  seeond  method  to  generate  a  Gaussian-like  beam.  Unlike  their 
bulk  eounterparts,  the  eonlinement  of  the  light  in  a  single-mode  eore  guarantees  high 
spatial  quality  output  over  a  wide  power  range  [10].  Beeause  they  have  long  interaetion 
lengths  and  low  losses,  liber  lasers  ean  have  very  low  thresholds.  Often,  the  fiber  ends 
are  used  as  the  eavity  mirrors,  thus  fiber  laser  deviees  have  potential  for  exeellent 
meehanieal  stability.  In  addition,  the  small  eross  seetion  ean  reduee  thermal  loading  on 
the  fiber  [11],  eliminating  mueh  of  the  need  for  eooling. 

Unfortunately,  nonlinear  effeets  ean  be  devastating  for  liber  lasers.  Although  most 
fiber  lasers  have  a  wide  speetral  bandwidth,  whieh  eliminates  SBS  as  a  power-limiting 
faetor,  stimulated  Raman  seattering  (SRS)  beeomes  a  problem  with  relatively  modest 
average  powers  in  pulsed  laser  systems.  SRS  inereases  the  wavelength  of  the  souree  as 
well  as  heats  up  the  liber.  At  higher  laser  powers,  eliminating  or  redueing  SRS  beeomes 
inereasingly  important. 

1.2,  Overview 

The  main  objeetive  of  this  dissertation  researeh  was  to  investigate  methods  to 
inerease  laser  intensity  through  the  manipulation  of  stimulated  seattering.  The  researeh 
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was  broken  into  two  parts.  The  first  was  beam  cleanup  in  optical  libers.  Beam  cleanup  is 
the  process  described  above  where  stimulated  Brillouin  scattering  converts  an  aberrated 
input  pump  beam  into  a  Gaussian-like  output  Stokes  beam.  The  second  phase  of  the 
research  was  aimed  at  reducing  stimulated  Raman  scattering  in  liber  lasers. 

The  investigation  into  beam  cleanup  began  experimentally.  A  laser  beam  was 
aberrated  with  an  etched  glass  plate  and  coupled  into  a  liber.  The  liber  output  was 
monitored  to  determine  how  well  it  fit  a  Gaussian  as  well  as  to  determine  its  polarization 
properties  and  measure  the  efficiency  of  the  cleanup  process.  This  was  followed  by  an 
attempt  to  reduce  threshold  by  feeding  a  portion  of  the  Stokes  output  back  into  the  fiber 
and  again  measuring  the  ability  of  the  fiber  to  clean  the  pump  beam. 

The  study  continued  with  an  investigation  into  beam  combining.  Essentially  this  is 
the  simultaneous  cleanup  of  two  pump  beams  in  a  single  liber.  Experiments  showed  that 
two  pump  beams  generated  a  spatially  coherent  Stokes  beam  regardless  of  the  spatial  or 
spectral  characteristics  of  the  pump  beams. 

The  experimental  exploration  of  beam  cleanup  and  beam  combining  was  completed 
with  an  investigation  into  some  of  their  limitations.  In  this  segment,  ideas  are  presented 
that  would  enable  beam  combining  of  more  than  the  two  sources  as  well  as  an  increased 
collection  efficiency  of  the  scattered  radiation.  Additionally,  the  second  order  threshold 
was  measured  to  determine  the  useful  power  range  for  beam  combining  or  cleanup. 
Einally,  the  possibility  of  using  stimulated  Raman  scattering  instead  of  stimulated 
Brillouin  scattering  was  investigated  in  the  hope  that  the  requirements  on  the  spectral 
quality  of  the  laser  beam  could  be  relaxed. 
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A  theoretical  analysis  mirrored  the  experimental  research.  The  process  of  beam 
cleanup  was  modeled  through  an  investigation  of  the  relative  gain  experienced  by 
different  fiber  modes  for  a  given  pump  beam  transverse  spatial  distribution.  This  sheds 
light  upon  the  physical  mechanism  responsible  for  beam  cleanup  and  demonstrated  some 
limitations  to  beam  cleanup  that  could  not  be  demonstrated  experimentally.  Furthermore, 
an  unexpected  polarization  property  of  the  Stokes  beams  formed  through  beam  cleanup 
and  beam  combining  was  explained  theoretically.  Finally,  the  second  order  threshold  for 
stimulated  Brillouin  scattering  was  modeled  allowing  for  an  estimated  threshold  under 
circumstances  other  than  those  of  the  experiment  and  leading  to  a  method  to  increase  the 
second  order  threshold  to  more  than  80  W. 

The  second  phase  of  the  research  was  designed  to  reduce  stimulated  Raman  scattering 
in  fiber  lasers.  Two  methods  were  proposed  to  limit  the  growth  of  the  Stokes  beam. 

First,  absorb  the  Stokes  beam  with  an  intercavity  dopant.  Second,  couple  the  Stokes 
beam  out  of  the  fiber  using  a  long  period  grating. 

The  investigation  into  an  intra-cavity  dopant  to  increase  SRS  threshold  was 
completed  theoretically.  The  dopant,  holmium,  was  chosen  based  upon  the  existence  of  a 
strong  absorption  near  the  Stokes  frequency  and  little  absorption  near  the  laser  frequency. 
The  absorption  curve  of  holmium  in  a  silicate  glass  window  was  measured  using  an  FTIR 
spectrometer.  The  expected  Stokes  threshold  was  then  modeled  as  a  function  of  the 
dopant  level  in  a  single-pass  geometry. 

Long  period  gratings  (LPGs)  are  similar  to  Bragg  gratings  in  that  they  are  periodic 
structures  created  in  a  photosensitive  fiber.  The  result,  however,  is  that  instead  of  a 
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spectrally  narrow  reflection  (~1  nm),  the  grading  couples  light  from  core-guided  fiber 
modes  into  cladding-guided  modes.  The  cladding  modes  of  a  fiber  are  typically  very 
lossy  so  the  light  in  these  modes  is  rapidly  attenuated  and  does  not  contribute  to  the 
exponential  growth  of  the  Stokes  beam.  The  investigation  into  this  began  by  measuring 
the  Raman  gain  curve  of  the  fiber,  then  designing  an  appropriate  LPG  to  eliminate  the 
Stokes  beam  from  the  fiber.  The  SRS  threshold  was  then  theoretically  modeled  to 
determine  the  effect  the  LPG  had. 

1.3.  Document  Organization 

The  research  presented  throughout  this  document  is  original.  Some  of  the  significant 
accomplishments  are: 

1.  SBS  Beam  cleanup  in  a  50-micron,  multimode,  graded-index  fiber. 

2.  SBS  Coherent  beam  combining  in  a  50-micron,  multimode,  graded-index  fiber 
with  two  independent  laser  sources. 

3.  SBS  Incoherent  beam  combining  in  an  8.2-micron,  multimode,  step-index  fiber. 

4.  SBS  second  order  Stokes  modeling  to  include  second  order  Stokes  seeding  by 
four- wave  mixing. 

5.  SRS  beam  cleanup  in  a  50-micron,  multimode,  graded-index  fiber. 

6.  Modeling  of  SRS  suppression  using  an  intercavity  holmium  dopant  in  a  silicate 
fiber. 

7.  Modeling  of  SRS  suppression  using  a  long  period  grating  in  a  germanium 
enriched  fiber. 

These  research  efforts  are  explained  in  detail  after  the  introductory  background  found  in 
Chapters  2  and  3.  These  two  chapters  outline  the  tools  required  for  the  rest  of  the 
dissertation.  Chapter  4  describes  the  experimental  and  theoretical  research  into 
stimulated  Brillouin  scattering  beam  cleanup.  Beam  combining  is  thoroughly  examined 
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in  Chapter  5.  The  limitations  to  beam  combining  and  beam  cleanup  are  investigated  in 
Chapter  6.  Chapters  7  and  8  describe  the  experiments  designed  to  modify  the  threshold 
of  SRS  in  optical  fibers  using  an  intercavity  dopant  or  long  period  grating  respectively. 
Final  conclusions  are  made  in  Chapter  9. 
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2.  Propagation  of  Light  in  Optical  Fibers 


2.1,  Core  Guided  Fiber  Modes 

When  an  electromagnetic  field  travels  through  an  optical  fiber,  it  is  constrained  by 
boundary  conditions  at  the  core-cladding  interface.  Only  special  field  distributions  meet 
these  boundary  conditions  and  can  propagate  efficiently  through  a  fiber.  These  field 
distributions,  called  optical  modes,  form  the  basis  of  much  of  the  research  that  is 
presented  in  this  dissertation.  An  understanding  of  their  origin  and  characteristics  is 
critical  to  the  remaining  chapters.  This  chapter  presents  the  derivation  of  the  optical 
modes  that  propagate  in  step-index  optical  fibers.  Similar  derivations  are  found  in  most 
texts  about  optical  waveguides  [12,  13].  While  graded-index  fibers  are  also  used  in  the 
experiments  described  in  this  document,  the  derivation  of  the  modes  is  not  described 
here. 

The  optical  fibers  used  in  the  experiments  described  in  this  dissertation  are  made  of 
doped  SiOi,  a  dielectric  material  that  is  nonmagnetic  and  has  no  free  charges  or  currents. 
In  this  fiber.  Maxwell’s  equations  and  the  constitutive  relations  reduce  to 


\xE  = - 

dt 

(2.1) 

VxB  =  Mo¬ 
ot 

(2.2) 

V-D  =  0 

(2.3) 

v..g  =  o 

(2.4) 
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(2.5) 


D  =  SqE  +  P 

In  order  to  decouple  the  equations  for  the  magnetic  and  electric  fields,  take  the  curl  of 
Equation  (2.1)  and  substitute  Equation  (2.2)  in  the  right  hand  side.  Apply  the  identity 

VxVxE  =  V(V.E)-V"E  (2.6) 


and  use  Equation  (2.5)  for  the  final  result  of 


(2.7) 


The  bulk  of  this  paper  will  describe  the  effects  of  the  nonlinear  aspect  of  the  fiber 
optic  waveguide,  but  in  order  to  analytically  calculate  the  electric  field  associated  with 
different  fiber  modes,  it  is  necessary  to  assume  the  fiber  is  linear  and  isotropic.  This 
assumption  allows  Equation  (2.7)  to  be  simplified  to 


n 

Kcj 


-^E  =  0. 
dt^ 


(2.8) 


We  are  looking  for  solutions  that  propagate  along  the  fiber  axis,  thus  we  can  assume  a 
solution  to  the  electric  field  that  has  the  form: 

E  =  E^  {r)Eg  (6’)exp[^/(k^z-®t]Je  (2.9) 


where  k  is  the  guided  wavenumber.  In  cylindrical  coordinates  Equation  (2.8)  becomes 


E^  (r)  dr 


Ee  [0)  = 


2  A 

k  "  - 

n 

—  CO 

g 

V 

) 

(2.10) 
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This  equation  can  be  separated  into  a  function  of  6  alone  and  a  function  of  r  alone  using 
a  separation  constant,  the  square  of  m.  This  results  in  an  azimuthal  field  dependence 
given  by  Eg{0')  =  exp \im6^ .  The  boundary  condition  for  this  equation, 

Eg[6^  =  Eg[6  +  forces  m  to  be  an  integer. 

The  azimuthal  and  longitudinal  equations  are  simple  to  solve,  reflecting  the  fact  that 
the  fiber  properties  do  not  vary  along  these  two  axes.  Along  the  radial  direction, 
however,  the  fiber  appears  as  a  series  of  three  steps.  The  core  of  the  fiber  is  a  perfect 
cylinder  with  a  high  index  of  refraction.  The  cladding  covers  the  core  and  has  a  slightly 
lower  index.  A  protective  jacket  typically  surrounds  the  cladding.  Each  of  these  regions 
is  critical  in  the  solution  to  the  radial  equation. 

For  efficient  propagation  of  light  through  the  waveguide,  it  must  be  totally  internally 
reflected  at  either  the  core-cladding  interface  or  the  cladding-jacket  interface.  If  the  total 
internal  reflection  occurs  at  the  core-cladding  interface,  the  light  propagates  in  a  core- 
guided  mode.  If  the  total  internal  reflection  occurs  at  the  cladding-jacket  interface,  the 
light  propagates  in  a  cladding-guided  mode.  It  is  convenient  to  initially  solve  for  the 
core-guided  modes,  assuming  that  the  cladding  is  infinite.  With  these  solutions,  the 
cladding-guided  mode  is  easier  to  solve. 

A  core-guided  mode  is  only  possible  if  the  light  rays  are  incident  upon  the  core¬ 
cladding  interface  at  an  angle  greater  than  the  critical  angle,  sin(^^ )  =  Xliis 

restricts  the  z  -component  of  the  wavevector  to  n^, adding  f-^g-  ^core  f  ■  Therefore,  the 
radial  part  of  the  function  will  be  different  depending  on  if  it  is  solved  in  the  core  or  in 
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the  cladding.  If  a  is  the  radius  of  the  core  and  we  define  and  y  by 


j  -  kg  ,  y^  =  k^  -  ^  ^  then  the  radial  component  of  the  field  is 


r  —  E  {r\  +  r^ 
dr  ’ 

r  —  E  (r)  +  r^ 
dr  ^  ’ 


dr^ 

d_ 

dr^ 


Er{r) 


V' 

f  ;52 


+ 


2  /  ^  2 
r  k„  -m 


£■  (r)  =  0  r  <  : 


(2.11) 


E^{r^  -  r^y^ +m^  (r)  =  0  r 


>  a 


The  equations  describing  the  field  amplitude  in  the  fiber  core  and  cladding  are 
Bessel’s  equation  and  the  modified  Bessel  equation  respectively.  The  solutions  are  given 
in  [14]  and  consist  of  a  linear  combination  of  Bessel  functions,  Weber  functions,  and 
modified  Bessel  functions.  By  forcing  the  field  amplitude  to  remain  finite  at  both  r  =  0 
and  r  =  GO  ,  the  solution  is  simplifies  to: 


Er{r)  = 


CJ^{kj)  r<, 
Pn,K,n{yr)  r>, 


(2.12) 


where  Jm,  and  Km  are  the  m  order  Bessel  functions  and  modified  Bessel  functions 


respectively  and  and  are  constants  controlling  the  power  in  the  m  fiber  mode. 


The  radial  components  of  the  electric  displacement  and  magnetic  induction  must  be 
continuous  across  the  core-cladding  interface  to  satisfy  the  boundary  condition.  A 
considerable  simplification  can  be  made  in  weakly  guiding  fibers  such  that 
^core  ~  ^cladding  1  •  In  thls  approxlmatlou,  the  boundary  conditions  can  be  rephrased  as 

.  This  approximation  will  result  in 

r=d^ 


[r)  =£■  (r) 

V  /  r=a  ‘  ^  ^  r 


and 


SEr(r) 


dr 


_  dE^r) 


dr 
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linearly  polarized  fiber  modes  [15].  These  boundary  eonditions  require  that 


For  eaeh  fiber,  this  term  is  a  eonstant  at  a  given  wavelength  and  is  used  to  relate  to  / . 
When  the  fiber  parameter  is  used  along  with  a  short-hand  notation,  X  =  k^a  ,  Equation 
(2.13)  beeomes 


(2.15) 


For  a  given  fiber  parameter,  eaeh  m  results  in  several  fiber  modes.  A  sample  is  shown  in 
Figure  1.  This  figure  represents  the  right  and  left  hand  sides  of  Equation  (2.15)  when 
V  =  5  .  As  this  figure  shows,  there  are  several  interseetions,  eaeh  of  whieh  represents  a 
mode  of  the  waveguide.  These  modes  are  labeled  by  an  integer  representing  the 
azimuthal  dependenee  of  the  field,  m  ,  and  an  integer  that  represents  whieh  interseetion  is 
to  be  used  (a  one  represents  the  interseetion  elosest  to  X  =  0 ).  The  seeond  number  is  a 
measure  of  the  number  of  zeros  in  the  radial  field  distribution.  The  intensity  patterns  of 
the  modes  represented  by  Figure  1  are  shown  in  Figure  2. 
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(  ^ 

!  1 :  Normalized  propagation  constant  for  the  modes  of  a  fiber  optic  waveguide, 
,  s^-J  with  V  =  5  . 


# 

LPol 

• 

LPo2 

m 

# 

•  # 

# 

LPii 

LP21 

^  transverse  mode  patterns  of  a  fiber  optic  waveguide. 

When  light  travels  in  a  pure  fiber  mode,  the  cross  sectional  intensity  distribution  in 
the  liber  is  constant  along  the  length  of  the  liber.  On  the  other  hand,  if  the  light  travels  in 
a  linear  combination  of  modes,  beating  between  the  modes  leads  to  a  changing  intensity 
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distribution  along  the  length  of  the  fiber.  To  illustrate  this  point,  Figure  3  shows  the 
electric  field  distribution  of  the  pure  LPoi  and  LPn  modes  (a  and  b  respectively)  along 
with  the  transverse  intensity  distribution  of  the  superposition  of  the  two  modes  (c)  at  the 
front  end  of  a  representative  50-micron  fiber.  Half  the  beat  length  further  down  the  fiber, 
(^coi  “  ^cii )  ’  the  relative  phases  between  the  two  modes  will  shift  by  1 80  degrees. 

Figure  3(d-f)  represents  this  case.  The  peak  intensity  within  the  fiber  has  moved  from  the 
right  side  of  the  fiber  to  the  left. 


3:  (a,b)  The  transverse  electrie  field  pattern  of  the  LPqi  and  LPn  modes  at  the 
Jnt  end  of  a  representative  50  mieron  fiber,  (e)  The  intensity  distribution  due  to  a 
superposition  of  the  modes  in  (a,b).  (d,e)  The  transverse  eleetrie  field  pattern  of  the 

LPoi  and  LPn  modes  at  a  distanee  of  ^/(k^oi  “  ^eii)  fiber,  (f)  The 

intensity  distribution  due  to  a  superposition  of  the  modes  in  (d,e). 


The  V  parameter.  Equation  (2.14),  is  not  only  important  in  determining  the  intensity 
distribution  of  the  modes,  it  is  also  critical  in  determining  the  number  of  modes  that  can 
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propagate  down  the  fiber.  In  the  above  example,  four  different  mode  patterns  are 
supported.  Eaeh  of  these  patterns  results  from  at  least  two,  and  more  eommonly  four, 
orthogonal  polarization  modes.  Thus  in  the  above  example,  12  modes  are  supported  by 
the  fiber,  two  polarization  modes  in  eaeh  of  the  LPoi  and  LP02  modes  and  four  in  eaeh  of 
the  LPii  and  LP21  modes. 


Short  of  using  a  graph  such  as  Figure  1,  there  is  no  exact  method  to  calculate  the 
number  of  modes  supported  by  the  fiber,  yet  it  is  very  important  to  know  what  modes  are 
supported.  This  is  the  only  way  to  ensure  that,  if  one  is  transmitting  light,  the  energy  can 
be  carried  by  the  fiber  or,  if  one  is  generating  light,  the  output  pattern  is  acceptable.  Two 
limits  are  described  here.  First  is  the  case  where  only  a  single  mode  is  desired.  As 
shown  in  Figure  1,  for  all  values  of  n,  the  right  hand  side  of  Equation  (2.15)  terminates 
when  X  =  V  .  Therefore,  when  V  <  2.405  ,  Equation  (2.15)  will  only  be  satisfied  for  one 
value  of  X  .  Although  there  are  two  orthogonal  polarization  modes  associated  with  this 
value  of  X ,  the  fiber  is  considered  single  mode  since  only  the  fundamental  Gaussian-like 
FPoi  intensity  distribution  is  supported. 


The  opposite  extreme,  when  V  is  large,  is  also  significant  and  the  number  of  modes 
that  can  propagate  in  a  fiber  can  be  approximated  as  described  below.  In  this  limit,  the 
left  hand  side  of  Equation  (2.15)  can  be  approximated  as  shown  in  Equation  (2.16)  [14]. 


8X 


sm 


X-^(l  +  2m) 


n 


(2.16) 
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This  function  is  zero  when  X  -  f  (l  +  2m)  =  nn  where  n  is  any  integer  and  each  of 

these  zeros  ean  be  assoeiated  with  a  mode  pattern  as  long  as  the  right  hand  side  of 
Equation  (2.15)  is  real.  So  by  eounting  the  number  of  zeros  between  X  =  0  and  X  =  V , 
the  number  of  allowed  modes  for  any  given  m  can  be  estimated.  Then  allowing  for 
different  polarizations  for  each  m  and  summing  over  m  ,  the  number  of  modes  is  [13] 

(2.17) 

n 

This  number  is  eritical  in  the  demonstration  of  beam  cleanup.  If  the  number  of 
supported  modes  is  small,  a  beam  will  be  spatially  filtered  before  entering  the  fiber.  Thus 
beam  cleanup  using  stimulated  scattering  cannot  be  investigated.  Only  when  a  large 
number  of  modes  are  supported  by  the  fiber  can  beam  cleanup  be  attributed  with 
eonfidence  to  stimulated  scattering. 

2.2,  Cladding-Guided  Fiber  Modes 

In  the  above  analysis,  we  assumed  an  infinite  cladding,  but  it  was  also  mentioned  that 
modes  guided  by  the  cladding-jacket  interface  could  also  propagate  in  the  fiber.  This  will 
only  oecur  if  there  is  total  internal  refieetion  at  the  eladding-jaeket  border.  Physieally, 
this  requirement  on  the  propagation  constant  is 

7  ^  kgcl  ^  ^cladding  f  '  (2-18) 

The  solution  to  the  eleetric  field  intensities  is  found  using  the  same  method  as  described 
for  the  eore-guided  modes  with  the  exeeption  that  the  radial  portion  of  Equation  (2.10)  is 
solved  using  boundary  eonditions  at  the  core-eladding  and  at  the  eladding-jaeket 
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interface.  Because  the  index  step  between  the  cladding  and  the  jacket  is  relatively  large, 
the  LP  mode  approximation  cannot  be  used  and  the  more  exact  forms  for  the  fields,  the 
HE  and  EH  modes,  must  be  used.  This  leads  to  a  dispersion  relation  that  is  much  more 
complicated  than  the  infinite  cladding  case.  Erdogan  derived  the  form  of  the  dispersion 
relation,  the  final  result  of  which  is  given  in  Appendix  A  [12,  16].  The  right  and  left 
hand  sides  of  the  equation  are  plotted  below. 


The  cladding-mode  intensity  distributions  are  calculated  from  the  effective  indices  in 
the  same  way  that  the  core-guided  intensity  distributions  were  found  (see  Appendix  A). 
Eigure  5  shows  the  first  three  circularly  symmetric  cladding  modes.  The  lowest  order 
core-guided  mode  is  shown  for  comparison. 
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5:  (Top  and  Bottom  Left)  Cladding  mode  intensity  distributions.  Odd  modes 
have  a  local  maximum  centered  in  the  fiber  core.  Even  modes  have  a  local 
minimum  in  the  fiber  core.  (Bottom  Right)  The  lowest  order  core  guided  mode. 
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3,  Stimulated  Scattering  in  Optical  Fibers 


Equation  (2.8)  of  the  preceding  chapter  neglects  the  contribution  of  the  nonlinear 
electric  susceptibility  to  the  electric  field.  This  simplification  is  valid  for  low  electric 
field  strengths,  but  it  is  not  in  the  realm  of  interest  for  the  rest  of  this  dissertation.  In 
situations  where  the  incident  light  is  intense  enough,  stimulated  scattering  is  observed. 
Stimulated  scattering  is  a  nonlinear  effect  in  which  the  energy  in  a  pump  photon  is 
redistributed  such  that  a  lower  energy  photon  and  a  phonon  are  emitted.  The  origin  of 
this  effect  is  examined  in  this  chapter  along  with  critical  features  that  will  be  applied  to 
different  situations  throughout  the  rest  of  this  document. 

Equation  (2.7)  provides  a  starting  point  to  derive  the  equations  governing  the 
coupling  between  the  pump  and  Stokes  beams.  The  polarization  term  in  this  equation  can 
be  accounted  for  by  expanding  it  in  a  power  series  of  the  electric  field 

P  =  s,  •  E  +  x^"^  ■■  EE  +  : EEE  +  •  •  •)  (3.1) 


where  is  the  permittivity  of  free  space  and  the  symbols  • ,  : ,  :  represent  the  product 

between  the  susceptibility  tensor  and  the  electric  field  (e.g.  the  component  of  the 
product  between  the  third  order  susceptibility  and  the  electric  field  is  given  by 


j,k,l 


The  third  order  nonlinear  susceptibility,  x^^^  ?  is  responsible  for  stimulated  scattering. 
Although  higher  order  terms  exist,  under  the  conditions  examined  in  this  document,  they 
are  negligible.  Thus  the  series  can  be  truncated  after  the  terms  shown  in  Equation  (3.1). 
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The  expansion  eoeffieients  z^^^  are  2“‘*,  and  order  frequeney 

dependent  suseeptibility  tensors,  respeetively. 

In  an  isotropie  material  sueh  as  glass,  any  rotation  or  inversion  of  the  material  should 
result  in  the  same  polarization  response.  Applieation  of  this  requirement  to  the  nonlinear 

suseeptibility  reduees  to  a  sealar,  requires  z^^^  to  be  identieally  zero,  and  reduees  the 

81  elements  that  make  up  the  4**^  rank  tensor  z^^^  to  3  independent  elements  [17] 


(3) 

(3) 

(3) 

(3) 


-  y(3)  -  y(3) 

/  yyyy  —  ^  zzzz 

-  y(3)  _  y(3) 

-  y(3)  -  y{3) 

yzyz  —  z>;2y  —  Jl  zxzx 

-  y(3)  _  y{3) 

—  Z  ~  Z  zxxz 


(3)  ,  (3) 

—  Z  ^yy  +  Z  w  +  Z 

(3)  (3) 

—  z  ^-zzz  —  Z  ^^yy  ~  X 

(3)  (3) 

—  Z  ^zxz  —  Z  zyxv  —  Z 

(3)  (3) 

—  Z  zzzx  —  Z  zjyx  —  Z 


(3) 

(3) 

(3) 

(3) 


xyyx 


yyxx 

yxyx 


yxxy 


(3.2) 


Thus  in  fiber  optie  waveguides  made  of  Si02,  the  wave  equation  beeomes 

V  A((l  +  Z )  E  +  Z  :£££).  (3.3) 

Sinee  we  are  looking  for  the  interaetion  of  a  pump  beam  and  the  Stokes  beam,  the 
eleetrie  field  ean  be  represented  by 

^  =  e;,z4^(z)exp[/(k^z-®/)]  +  e,4(z)exp[/(±k,z-®^t)]  +  c.c.  (3.4) 


where  e.e.  represents  the  eomplex  eonjugate  of  the  previous  terms,  {k^)  is  the 
wavenumber  of  the  pump  (Stokes)  field,  0^(0^)  is  the  frequeney  of  the  pump  (Stokes) 
field,  Ap  (z)  ( A^,  (z)  )  is  the  eomplex  pump  (Stokes)  field  amplitude,  and  is  the 

eomplex  polarization  unit  veetor  of  the  pump  (Stokes)  field.  The  ±  before  the 
waveveetor  in  the  Stokes  term  integrates  the  effeet  of  the  direetion  of  the  Stokes  wave 
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into  the  equations.  This  representation  of  the  eleetric  field  ignores  all  transverse  (x  and  y) 
dependenee  of  the  field  amplitude.  The  transverse  term  is  assoeiated  with  the  fiber 
modes  and  has  already  been  investigated  in  Chapter  2. 

Equation  (3.4)  is  a  veetor  equation  that  describes  the  electric  field  in  the  fiber  and  its 
propagation  down  the  fiber  axis.  However,  the  equation  is  difficult  to  work  with  in  its 
current  form.  For  the  remainder  of  the  derivation,  a  specific  polarization  of  the  field  will 
be  assumed,  reducing  the  above  to  a  scalar  equation. 

The  wave  equation  contains  the  third  power  of  the  electric  field  which,  when 
expanded,  results  in  many  terms.  In  order  for  the  equation  to  hold  at  all  times,  the  terms 
that  oscillate  at  like  frequencies  can  be  set  equal  to  each  other.  The  terms  of  interest  in 
stimulated  Raman  scattering  are  those  that  oscillate  at  and  .  These  terms  result  in 

the  equation  shown  below. 


^2 

[o)p,s )  +  [o)p,s ,  )\Ap,,  f  +  6  | 


+ 


(3.5) 


Another  term  taking  four-wave  mixing  into  account  must  be  added  to  this  equation  if  the 
Stokes  beam  propagates  in  the  same  direction  as  the  pump  beam.  While  not  critical  to 
this  derivation,  it  becomes  very  important  when  second  order  SBS  is  analyzed  in 
Chapter  6. 

Since  the  x  -  and  y  -variation  in  the  electric  field  was  ignored,  the  Laplacian  operator 
reduces  to  the  second  derivative  with  respect  to  z .  In  nearly  all  applicable  situations,  the 
variation  in  the  amplitude  of  the  electric  field  will  be  small  in  comparison  to  the  variation 
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of  the  exponential  portion  of  the  field.  Assuming  that  ^y4(z)|  |A:^y4(z)| ,  evaluating 

the  Laplaeian,  the  time  derivative,  and  eaneeling  the  exponent  on  both  sides  results  in 

{l  +  )  +  •  •  • 

^  .  (3.6) 

•  •  •  +  )  1^;,,,  f  +  6z^'^  f }  =  0 

The  funetional  form  of  the  third  order  nonlinear  susceptibility  is  such  that  the 
magnitude  of  z^^^  ^  ~(^s,p  ’  )  will  be  greatest  when  there  is  an  excited  energy  level 

that  is  precisely  /z  )  above  the  ground  state.  Likewise,  there  will  be  a  resonance 

for  (the  susceptibility  associated  with  two-photon  absorption  [17]), 

when  there  is  an  energy  level  that  is  exactly  2/z®^  ^  above  the  ground  state.  We  have 

specifically  chosen  co^  such  that  the  first  condition  is  true  (in  silica  fibers  this 

corresponds  to  a  frequency  separation  of  about  13.2  THz  for  SRS  and  16  GHz  for  SBS), 
thus  it  is  unlikely  that  the  second  condition  will  also  hold.  The  two-photon  absorption 
term  will  therefore  be  negligible  compared  to  the  stimulated  Raman  scattering  term.  To 

further  simplify  the  equation,  multiply  Equation  (3.6)  by  -iA^  ^  and  add  the  complex 

I  |2  * 

51  zi  *  QA  dA 

conjugate  of  the  resulting  equation.  Noting  that  =  A^  ,  +  A^  ,  ,  this  yields; 
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Finally  it  should  be  noted  that  the  imaginary  part  of  the  linear  suseeptibility  is 
responsible  for  absorption  [3].  The  gain  eoeffieient,  g ,  is  due  to  the  imaginary  part  of 
the  third  order  suseeptibility  [17],  the  intensity  of  the  pump  and  Stokes  beams  is 
proportional  to  the  amplitude  of  the  E-field  squared,  and  the  imaginary  part  of  the  third 
order  suseeptibility  is  an  odd  funetion  of  detuning  between  the  frequeney  eomponents 

[IV]. 


^P,sC 


^  = - 

s^c 


Im 


Im 


(3.8) 

(3.9) 

1  |2 

1  \a 

'  p^s  p,s 

(3.10) 

(3.11) 

Where  a  is  the  absorption  eoeffieient,  and  I ^  ^  is  the  intensity  of  the  pump  or  Stokes 

beam.  Thus  the  intensity  of  the  Stokes  beam  in  a  fiber  is  given  by  the  following  set  of 
eoupled  equations  [3] 


Oh 

dz 

dz 


(3.12) 


Near  threshold,  when  the  pump  power  is  mueh  greater  than  the  Stokes  power,  it  ean 
be  assumed  that  there  is  no  pump  depletion  due  to  stimulated  seattering.  While  this 
violates  the  law  of  eonservation  of  energy,  it  allows  for  an  analytieal  solution  to  the 
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above  system  of  equations.  Under  this  approximation,  the  solution  for  the  intensity  of  the 
Stokes  beam  as  a  funetion  of  z  is 


T  If 

gjpo - 

1 

-az 


(3.13) 


where  1^^  is  the  intensity  of  the  Stokes  beam  ineident  on  the  fiber  faee.  If  the  Stokes 
beam  is  traveling  in  the  opposite  direction  as  the  pump  beam,  then  1^^  is  the  intensity 
incident  upon  the  rear  face  of  the  fiber.  If  the  Stokes  beam  is  propagating  in  the  same 
direction  as  the  pump  beam,  then  1^^  is  the  intensity  incident  upon  the  front  face  of  the 

fiber.  The  direction  of  the  traveling  beam  is  chosen  by  the  selection  of  either  the  plus  or 
minus  in  the  equation  for  the  Stokes  intensity.  The  plus  corresponds  to  a  wave  traveling 
in  the  same  direction  as  the  pump  beam. 

This  result  demonstrates  the  obvious  complication  that  SRS  can  impose  on  a  laser 
system.  The  exponential  increase  in  Stokes  power  will  quickly  deplete  the  energy  in  the 
pump  and  significantly  reduce  the  laser  efficiency.  By  their  nature,  the  cavity  of  a  laser 
system  will  have  extremely  high  intensities,  thus  even  when  there  is  no  injected  Stokes 
beam,  the  spontaneous  emission  of  light  through  blackbody  radiation  or  spontaneous 
Raman  scattering  can  be  strong  enough  to  start  the  stimulated  scattering  process  and 
reduce  the  laser  power  significantly. 

For  stimulated  Raman  scattering  to  have  a  significant  effect,  the  amount  of  power  at 
the  Stokes  beam  must  be  above  threshold.  Experimentally,  the  threshold  is  easily 
observed  as  a  sharp  knee  in  the  Stokes  output,  above  which,  the  Stokes  power  increases 
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linearly.  Theoretically,  it  is  much  more  difficult  to  identify  the  threshold  criteria.  Smith 
defines  it  as  the  pump  power  that  predicts  an  equivalent  Stokes  output  based  upon  the 
undepleted  pump  approximation  of  Equation  (3.13)  [18].  Although  the  no  pump 
depletion  approximation  is  obviously  invalid  at  this  point,  the  error  is  relatively  small 
when  estimating  the  threshold. 

Based  on  Equation  (3.13)  the  amount  of  Stokes  power  at  the  exit  end  of  the  fiber  is 
directly  proportional  to  the  amount  at  the  input  end  of  the  fiber.  Typically,  however,  the 
Stokes  power  builds  from  noise.  Eor  long  fiber  lengths,  when  ~  cCp,  the  noise  can  be 

modeled  by  injecting  a  fictitious  photon  at  a  certain  location  in  the  fiber  [18].  Eor 
forward  scattering,  the  photon  should  be  injected  at  the  front  end,  E  =  0  .  Eor  backward 

scattering  processes,  this  photon  should  be  injected  at  E  =  ^In^g^  ■ 


To  calculate  the  total  power  in  the  Stokes  beam,  the  Stokes  power  at  every 
wavelength  component  must  be  integrated  together.  In  the  forward  scattering  case  this 
integral  becomes 


Ps{L)  =  \\h(OQX^ 


(0) 


l-exp(-«^E) 


-aL 


•do) , 


(3.14) 


where  the  integration  is  performed  over  the  Raman  gain  bandwidth  and  the  pump  is 
assumed  to  be  monochromatic. 

In  order  to  solve  the  above  integral,  it  is  necessary  to  have  some  knowledge  of  the 
functional  form  of  the  gain  coefficient.  The  Raman  gain  spectrum  will  have  a 
predominant  peak  associated  with  a  vibrational  state  in  the  material.  The  spectrum  is 
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typically  a  complicated  function  of  the  frequency  difference  as  shown  for  Si02  in  Figure 
6.  It  would  be  extremely  diffieult  to  solve  Equation  (3.14)  analytically  if  the  exaet  form 
of  the  Raman  gain  was  used. 


Figure  6:  Raman  gain  spectrum  [19] 


Assuming  a  Lorentzian  gain  profile,  equal  absorption  at  the  pump  and  Stokes 
wavelengths,  and  long  fibers,  the  forward  stimulated  Raman  scattering  threshold  is 
approximated  by  [17,  18,  20]; 

ai5) 

SrO^eff 


Where 


t[l-exp(-«p^) 


is  the  effective  length  of  the  fiber.  The  result  for 


baekward  scattering  is  slightly  different  because  the  Stokes  photon  that  is  introduced  to 
simulate  spontaneous  emission  is  not  injected  at  the  end  of  the  fiber.  For  stimulated 
Raman  seattering,  the  backward  scattering  threshold  is  approximately 
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(3.16) 


A 


-4f 


SrO^eff 


From  these  two  equations,  it  is  apparent  that  forward  stimulated  Raman  seattering  will 
reaeh  threshold  before  baekward  SRS  and  is  mueh  more  likely  to  be  the  dominant  effeet. 

Stimulated  Brillouin  seattering  on  the  other  hand  is  only  possible  in  the  baekward 
direetion.  The  threshold  for  this  proeess  is 

(3.17) 

SrO^eff 


These  equations  are  approximate  and  ean  ehange  under  different  eonditions.  If  the 
absorption  at  the  Stokes  and  pump  wavelengths  are  different,  the  fiber  is  not  long  enough 
to  foree  a^L  »  1 ,  or  the  effeetive  length  of  the  fiber  is  different  at  the  Stokes  wavelength 

and  the  pump  wavelength,  the  derivation  may  break  down.  These  situations  will  be 
investigated  in  the  experiments  deseribed  in  Chapters  7  and  8. 
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4,  SBS  Beam  Cleanup 


4,1,  Background 

Recent  scientific  literature  describes  many  methods  to  clean  laser  beams  through 
stimulated  Brillouin  scattering.  Typically,  however,  they  refer  to  the  process  of 
eliminating  the  aberrating  effects  of  an  optical  component  such  as  an  amplifier.  In  these 
experiments,  after  the  beam  travels  through  the  aberrating  element,  it  is  phase  conjugated 
using  SBS  and  sent  back  through  the  aberrator  [21].  This  second  pass  removes  the  phase 
aberrations  from  the  wavefront  and  the  beam  returns  to  its  original  shape. 

Beam  cleanup  as  discussed  in  this  dissertation  is  significantly  different  in  that  it  refers 
to  exciting  an  LPqi  fiber  mode  through  SBS  using  an  aberrated  laser  beam.  The 
fundamental  difference  is  that  this  method  can  clean  any  aberrated  beam,  whereas  the 
phase  conjugation  experiments  described  above  cannot  clean  the  aberrations  inherent  to 
the  laser. 

The  process  of  exciting  low  order  fiber  modes  from  aberrated  sources  through 
nonlinear  processes  is  not  a  new  concept.  In  1981,  Nesterova  et  al.  showed  that  the 
transverse  Stokes  power  in  an  optical  fiber  was  redistributed  through  stimulated  Raman 
scattering  from  the  highly  aberrated  pump  beam  into  low  order  fiber  modes  [22].  Similar 
results  were  reported  by  Baldeck  et  al.  in  1987,  when  they  noted  that  an  aberrated  pump 
produced  a  ring-like  Stokes  beam  [23].  In  1988,  Grundinin  et  al.  found  that  solitons 
generated  from  SRS  in  multimode  fibers  propagated  in  the  fundamental  mode  [24]. 
Furthermore,  in  1992,  Chiang  noticed  that  when  multiple  Stokes  modes  were  excited 
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through  SRS  in  fiber  waveguides,  that  subsequent  higher  order  Stokes  beams  propagated 
in  fewer  optieal  modes  than  the  previous  one  [25]. 

In  1993,  Brusselbaeh  reported  similar  effeets  using  stimulated  Brillouin  seattering 
[26].  Later,  Rodgers  et  al.  showed  how  the  proeess  is  useful  for  beam  eleanup  and  beam 
eombining  [27].  Rodgers  demonstrated  SBS  beam  eleanup  by  transmitting  a  beam 
through  a  mieroseope  slide  etehed  in  aeid  and  effieiently  eleaning  the  beam  through 
stimulated  Brillouin  seattering  in  a  multimode  optieal  fiber. 

This  ehapter  further  examines  the  phenomenon  of  stimulated  Brillouin  seattering 
beam  eleanup.  Seetion  4.2  studies  the  theoretieal  souree  of  beam  eleanup.  Seetion  4.3 
experimentally  examines  the  properties  of  the  Stokes  beam.  Conelusions  are  made  in 
Seetion  4.4. 

4,2,  Theoretical  Beam  Cleanup 

As  described  in  Chapter  3,  the  interaction  of  the  pump  beam  and  the  Stokes  beam 
through  stimulated  Brillouin  scattering  lead  to  a  set  of  coupled  differential  equations 
given  by  Equation  (3.12).  These  are  repeated  here  with  the  approximation  that 
cOp!  0)^  —  1  and  explicitly  stating  the  dependence  of  the  intensity  on  z  and  the  transverse 

dimensions  . 

=  -gsls  ('k .  2) (l .  ('k . 

f  (4.1 

=  Sb^s  {ri,z)lp{r^,z)-  a  pip  (r^ ,  z) 


29 


Integrating  Equations  (4.1)  over  the  transverse  dimensions  and  letting 


z'^dr^  be  the  power  flowing  through  the  fiber  at  point  z  yields; 

=  -gBeff  {^)Ps  {^)Pp  (^)  (^) 

=  -gBeff  {^)Ps  {^)Pp  {^)-(^pPp  (^) 


gBeff  gB 


Ps{-)Pp{-) 


(4.3) 


Equation  (4.3)  defines  the  effective  gain  and  clearly  shows  its  dependence  on  the  overlap 
of  the  Stokes  and  pump  beams. 

The  pump  and  Stokes  electric  fields  can  be  expanded  in  terms  of  the  fiber  modes 
described  in  Chapter  2.  After  a  little  algebra  (see  Appendix  B),  the  effective  Brillouin 
gain  becomes 


gBeff  (^) 


-  s  \  ^  A  ( z  \  y  H- 

I  ^pj  \  )•  sm,sm,pj\pj  **• 

I  l^tnj 


g,(4goCn)' 

2  y  A  .  ■(z)y  .  .cos(aA:  zA-t^S  )  +  ... 

sm,sq,pj ,pj  V  //  sm,sq,pj,pj  y  smq  rsmq  J 


m,q<m,j 

^  ^  !  '^sm,sm,pj ,pv  y sm,sm,  pj,pv  cos(Akp^.,z  +  A^/i^^.,)  +  ... 


(4.4) 


m,j,v<  j 

2  A  (z)y  cos((AA:  .  -  Ak  )z  +  (A(zi.  -Ad  ))  +  ... 

sm,sq,pj,pv  \  // sm,sq,pj,pv  pjv  smq  J  y  '  P]'^  rsmq  Jj 

m,q<mj  ,v<  j 

^  ^  -I  ^sf^^sq,pj ,pv  sm,sq,  pj,pv  cos  (( Ak^^.,  +  Ak,„^ )  z  +  [A(I)^j^  +  A^,„^ ))  i 

m,q<mj,v<j  j 
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(4.5) 


rsn,,s,,pj,pv  =  j  [Wsm  (''l  )  ¥ sq  )  ¥  pj  (''l  )  ¥  p.  (''l  )} 

An,,s,,pj,p.  (^)  =  Am  (^)  A,  (^)  (^)  ^pv  (^) 

.  =  A:  ■  -  k 

pjv  pcj  pcv 

hJc  =  k  -k 

smq  scm  scq 


^(l>pjv  =  (t>pcj  -  (l^pcv 
^Amq  =  Acm  -  Acq 


where  the  summations  are  over  all  fiber  modes  and  the  m  *  mode  has  a  transverse 
distribution,  y/^  (r,^) ,  an  amplitude,  (z) ,  a  propagation  constant,  k^  ,  and  a  phase, 

(j)^  .  Unfortunately,  the  z  -dependence  of  the  modal  amplitude  coefficients  is  dependent 

upon  the  effective  Brillouin  gain  so  Equation  (4.4)  does  not  immediately  reveal  the  mode 
that  will  produce  the  greatest  gain.  If,  however,  it  is  assumed  that  all  pump  modes  decay 
in  the  positive  z  -direction  at  the  rate  (z) ,  and  all  Stokes  modes  decay  in  the  positive 

z  -direction  at  the  rate,  (z) ,  the  field  amplitude  can  be  written  as  a  product  of  a 

constant  and  the  decay  function.  This  assumption  is  only  valid  for  low  pump  powers, 
where  the  Stokes  gain  is  minimal.  Above  threshold,  experimental  evidence  shows  that 
the  Stokes  LPoi  modal  amplitude  grows  much  faster  than  the  others.  This  analysis, 
however,  still  conveys  insight  into  the  increased  gain  of  the  LPoi  mode  over  other 
intensity  distributions.  With  this  assumption  the  modal  amplitude  is  written  as 

Am,sq,pj,pj  (^)  =  d,  (zf  {zf  and  the  product  of  the  Stokes  and  pump  power 

becomes  {^)Pp  (^)  =  d^  (z)^  d^  (z)^  (O)  •  Substituting  this  into  Equation  (4.4) 

reduces  all  z  -dependence  in  the  gain  to  the  cosine  terms. 
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SBeff  (2)  = 


gB  [^s^cnf 


r.r. 


sm,sm,pj ,pj  t  sm,sm,pj,pj 


+  ... 


2  E  ^sm,sq,pj,pjT,m,sq,pj,pj  ^OS  [Ak^^^Z  +  )  +  . .  . 

m,q<mj 

2  ■  y  cosiAk  z-\-  Ad)  .  )  +  ...  (4-6) 

/  i  sm,sm,pj,pv/ sm,sm,pj,pv  \  pjv  '  Pjv  ) 

m,j^v<  j 

2  E  '^sm,.q,pj,pvr.m,sq,pj,p.  COS  (( -  Ak^^^ )  z  +  ( ))  +  ... 

m,q<m,j  ,v<  j 

2  E  ^sm,.q,pj,pvr.n,,sq,pj,pv  COS  (( AA:^.^  +  Ak^^^ )  z  +  ( +  A^,„^ ))  j 


m^q<m,j  ,v<  j 


Figure  7(a)  is  a  plot  of  the  effective  gain  along  the  length  of  a  fiber  when  the  pump 
and  Stokes  beams  are  made  of  LPoi  and  LPn  modes.  The  rapid  oscillations  are  due  to  the 
terms  in  Equation  (4.6)  fluctuating  at  spatial  frequencies  of  Ak  or  greater.  These 
fluctuations  will  quickly  average  to  zero  leaving  only  the  slowly  oscillating  terms  of 
Equation  (4.7). 


SBeJf{^)  = 


g^(4goCn)" 

2  E 


1  sm,sm,pj,pj/  s 


sm,sm,pj  ,pj 


+  ... 


m,j 


(4.7) 


m,q<m,j  ,v<  j 


sm,sq,pj,pvrsm,sq,pj,pv  COS  (( Ak^.^  -  Ak,„^ )  z  +  ( A^.^  -  A^,„^  )) 


This  is  plotted  with  a  solid  line  in  Eigure  7(b).  These  slow  oscillations  in  the  gain  are  due 
to  modal  dispersion.  As  described  in  Chapter  2,  the  transverse  power  distribution  in  a 
fiber  changes  along  the  fiber  axis  (see  Eigure  2).  Because  these  oscillations  are 
frequency  dependent,  if  the  pump  and  Stokes  beams  overlap  well  at  the  fiber  front  face, 

they  will  have  a  minimal  overlap  at  some  point  nji^Ak^  -  Ak^  j  in  the  fiber.  Pure  fiber 
modes  do  not  oscillate  in  the  fiber,  thus  their  overlap  does  not  change  along  the  fiber 
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length.  The  dashed  line  is  the  relative  gain  when  the  Stokes  beam  is  in  the  LPoi  pure 
mode. 


Figure  7:  Relative  effective  gain  along  the  length  of  the  fiber,  (a)  is  the  calculated  value, 
(b)  is  the  calculated  value  when  the  high  frequency  fluctuations  are  averaged  out. 


Comparing  the  solid  and  dashed  lines  of  Figure  7(b)  sheds  some  light  on  the 
mechanism  that  is  responsible  for  beam  cleanup  in  long  fibers.  In  a  short  fiber,  the 
Stokes  beam  would  propagate  in  the  same  linear  combination  of  modes  as  the  pump 
because  of  the  greater  gain.  In  a  longer  fiber  however,  the  Stokes  beam  propagating  in 
the  LPoi  pure  mode  has  a  greater  average  gain. 

Given  Equation  (4.7),  the  average  gain  over  the  fiber  length  is  a  simple 
computation. 


Sseff 


[l 

E,(0)/),(0)[2 


sm,s, 


.sm,pjypj/  sm,sm,pj ,pj 


+... 


y  ^ 

sm,s 


m,j 

y.r. 


sq,pj,pv/  sm,sq,pj,pv 


cos 


m^q<m,j  ,v<  j 

COS 


) sine (( )  c)  + . . .  (4.8) 


(( Ak^,v  -  )l) sin ( A^,  -  ) 
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The  Stokes  mode  that  will  be  exeited  for  a  given  pump  can  be  predicted  from  this 
equation.  The  pump  mode  distribution  fixes  the  values  for  all  /i p.  and  (/)pj .  Thus,  the 

only  variables  left  in  the  equation  are  the  amplitudes  and  phases  of  the  N  Stokes  modes 
excited  in  the  fiber. 

Solving  for  the  Stokes  field  in  this  manner  is  computationally  intensive.  A  typical 
50-micron  fiber  will  support  more  than  100  modes  with  a  1.064-micron  pump,  so  there 
are  hundreds  of  independent  variables  to  determine.  In  order  to  cut  down  on  the 
computation  time,  only  four  modes,  LPoi,  LPn,  LP21,  and  LP02,  were  allowed  to 
propagate  in  the  fiber  and  each  mode  was  assumed  to  have  a  phase  of  zero. 

The  results  of  the  analysis  for  an  infinitely  long  fiber  show  that  a  pure  mode  was 
excited  100%  of  the  time.  In  other  words,  when  the  input  was  a  linear  combination  of  all 
four  modes,  the  Stokes  beam  had  the  strongest  gain  when  it  propagated  in  one  of  the  pure 
fiber  modes.  The  LPoi  mode  had  the  greatest  gain  for  more  than  18%  of  the  possible 
pump  intensity  distributions. 

It  has  been  demonstrated  that  short  multimode  fibers  produce  phase  conjugate 
replicas  of  the  pump  beams  [28-30].  Thus,  in  the  limit  that  the  fiber  length  approaches 
zero.  Equation  (4.8)  should  predict  that  the  Stokes  beam  with  the  highest  gain  has  the 
same  modal  power  distribution  as  the  pump  beam.  Over  200  equally  spaced  pump  mode 
distributions  were  tested.  The  average  phase  conjugate  fraction  was  87%.  The  phase 
conjugate  fraction  is  defined  by  Equation  (4.9)  where  modal  amplitudes  are  normalized 
such  that  the  power  in  both  the  pump  and  Stokes  beams  has  a  magnitude  of  one.  This 
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ensures  a  phase  conjugate  fraction  of  one  when  the  Stokes  beam  has  the  same  transverse 
intensity  distribution  as  the  pump  regardless  of  the  total  power  in  the  Stokes  beam. 


(4.9) 


Figure  8(a)  is  a  plot  of  the  phase  conjugate  fraction  as  a  function  of  fiber  length. 
Figure  8(b)  illustrates  the  percentage  of  the  possible  pump  modes  resulting  in  a  Stokes 
beam  in  the  LPoi  mode.  Both  of  these  plots  show  the  expected  trend  of  an  increase  in  the 
beam  cleanup  and  a  reduction  of  the  phase  conjugate  fraction  as  the  fiber  length 
increases. 
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(a)  Average  phase  eonjugate  fraetion  as  a  funetion  of  the  fiber  length,  (b) 
—  JPereentage  of  pump  modes  that  lead  to  the  exeitation  of  the  LPoi  Stokes  mode 

versus  the  fiber  length. 


In  practice,  SBS  beam  cleanup  in  long  fibers  appears  to  be  much  more  robust  than 
predicted.  As  explained  above,  Equation  (4.8)  predicts  that  only  18%  of  the  possible 
pump  modes  generate  an  LPoi  Stokes  beam  in  an  infinitely  long  fiber.  In  the  laboratory, 
however,  the  LPoi  mode  was  excited  almost  100%  of  the  time  as  will  be  discussed  in  the 
remaining  sections  of  this  chapter.  There  are  several  possible  reasons  for  this  apparent 
disconnect. 
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First,  it  was  assumed  that  the  mode  strueture  of  the  pump  and  Stokes  beams  would 
not  evolve  along  the  fiber  length.  Even  ignoring  the  nonlinear  mode  coupling,  power  will 
flow  between  modes  because  the  fiber  is  not  a  perfect  cylinder.  The  fiber  is  wound 
around  a  spindle  with  an  8.5  cm  radius  of  curvature,  small  imperfections  in  the  index  of 
refraction  and  the  core  concentricity  are  unavoidable  in  fabrication,  and  microbends  all 
perturb  the  fiber  [31-36].  Kutz  analyzed  these  factors  and  found  that  the  imperfections  in 
the  core  concentricity  are  the  leading  cause  of  modal  power  diffusion  [33].  He  found  that 
in  a  graded-index,  multimode  fiber,  the  power  in  each  fiber  mode  obeys 

^  =  4  t  {n  +  l)(m+l)d,„[P„(z)-P,(z)l  (410) 

CIZ  ZK  rn=l[m^n'j 

where 


=  C  {gWn  ,  /  2)" 


(4.11) 


In  Equations  (4.10)  and  (4.1 1),  n  and  m  are  mode  numbers,  is  the  difference  in 
propagation  constants  between  the  two  modes,  is  the  power  in  the  n  *  fiber  mode,  g 

is  the  transverse  perturbation  to  the  index  of  the  fiber,  D  is  related  to  the  width  of  the 
Gaussian  noise  correlation  function,  C  is  a  constant,  and  s  is  the  magnitude  of  the 
perturbation.  Several  features  of  Equation  (4.10)  are  important  in  this  analysis.  Eirst  the 
coupling  coefficient,  ,  decays  exponentially  with  the  difference  between  the 

propagation  constants  of  the  two  modes.  In  other  words,  the  EPoi  mode  will  not  decay 
into  the  EPgi  mode  as  quickly  as  it  will  decay  into  the  EPn  mode.  Secondly,  the  higher 
order  modes  will  change  more  rapidly  than  the  lower  order  modes  due  to  the  degeneracy 
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term  in  Equation  (4.10),  (n  +  l)(m  +  l) .  Finally,  assuming  the  eoupling  eoefficient  is 

nonzero,  the  steady  state  eondition  exists  when  power  is  uniformly  shared  between  all 
fiber  modes.  In  other  words,  a  multimode  pump  beam  will  quiekly  distribute  power 
throughout  many  fiber  modes  ereating  a  uniform  pump  intensity  along  the  fiber  length. 
However,  the  Stokes  beam  that  is  generated  in  the  LPoi  mode  will  not  couple  to  higher 
order  modes  for  much  greater  distances.  This  asymmetry  in  the  mode  coupling  can 
explain  why  much  more  than  18%  of  the  pump  intensity  distributions  result  in  beam 
cleanup,  but  it  is  highly  dependent  on  the  index  distribution  in  the  fiber  [34,  35].  Mode 
mixing  in  step  index  fibers  is  strongest  in  low  order  modes,  thus  beam  cleanup  in  these 
fibers  may  be  more  difficult  to  attain  reliably. 

Similarly,  there  will  be  a  mixing  of  modes  induced  by  the  scattering  process.  In  the 
beam  cleanup  analysis,  it  was  assumed  that  all  pump  modes  decayed  at  the  same  rate.  In 
a  physical  fiber,  however,  the  modes  will  not  uniformly  decay,  but  will  change  relative  to 
each  other  since  the  Brillouin  reflection  is  not  uniform  across  the  transverse  fiber 
dimensions.  One  method  to  partially  lift  this  restriction  is  to  assume  that  the  pump  beam 
evolves  through  all  possible  mode  distributions  so  quickly  that  it  appears  to  be  in  all 
modes  at  once.  In  this  manner,  the  Stokes  mode  with  the  greatest  average  overlap  over 
all  possible  pump  modes  would  dominate.  The  gain  of  a  given  Stokes  mode  in  a  long 
fiber  under  these  conditions  is  shown  in  Equation  (4.12),  where  is  the  modal  volume 

and  the  integral  is  over  all  possible  combinations  of  pump  modes,  /4p  .  Again  only 
allowing  four  modes  in  the  fiber  and  assuming  the  phase  of  each  mode  is  zero,  the  LPoi 
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pure  mode  will  have  a  greater  gain  than  any  other  mode  distribution  in  an  infinite  fiber, 
regardless  of  the  input. 
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Additionally,  the  Brillouin  gain  in  the  fiber  was  assumed  to  be  eonstant.  This  is  a 
relatively  safe  assumption  along  the  transverse  direetion  in  a  step-index  fiber,  but  a 
graded-index  fiber  may  behave  remarkably  different.  The  index  is  eontrolled  by 
adjusting  the  eoneentration  of  dopants  in  the  eore  of  the  eladding.  For  example,  the 
addition  of  germanium  to  the  eore  of  a  fiber  inereases  the  index  of  refraetion  relative  to 
that  of  pure  Si02.  To  ereate  a  graded  index  fiber,  the  germanium  eontent  would  drop 
from  its  highest  value  at  the  eore  eenter  to  the  lowest  value  in  the  eladding.  A  side  effeet 
of  this  graduated  doping  level  is  that  sinee  the  Brillouin  gain  eoeffieient  is  greater  for 
Ge02  than  for  Si02,  modes  with  higher  intensities  at  the  eenter  of  the  fiber  will  be 
preferentially  exeited. 

Finally,  it  was  assumed  that  all  fiber  modes  had  the  same  attenuation.  In  reality,  this 
is  typieally  not  true.  The  higher  order  modes  have  more  power  near  the  eore-eladding 
interface  and  therefore  typically  have  a  higher  scattering  loss  [35,  37,  38].  Since  the 
Stokes  beam  grows  exponentially  in  the  fiber,  even  small  differences  in  the  loss  can 
significantly  change  the  overall  distribution  of  power  in  the  Stokes  modes. 
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In  summary,  a  theoretical  model  of  Brillouin  gain  was  presented.  The  power 
transmitted  through  the  fiber  was  decomposed  into  the  modal  amplitudes  and  the  gain  for 
each  Stokes  mode  and  combination  of  modes  was  calculated.  The  modal  distribution 
with  the  greatest  average  gain  was  assumed  to  be  the  dominant  Stokes  output.  This 
analysis  produced  the  trends  seen  in  the  laboratory.  Short  fibers  reflected  a  phase 
conjugate  replica  of  the  pump  whereas  long  fibers  produced  output  in  the  pure  fiber 
modes.  The  strength  of  beam  cleanup  was  not  fully  predicted,  but  limits  in  the  analysis 
were  presented  that  could  explain  the  differences  between  theory  and  experiment. 

4.3.  Experimental  Beam  Cleanup 

The  properties  of  the  cleaned  Stokes  beam  generated  in  a  multimode  fiber  are 
experimentally  examined  here.  The  experiments  described  in  Section  4.3.1  verify  that 
SBS  beam  cleanup  is  viable  in  a  50-micron  multimode  fiber.  This  is  followed  by  an 
investigation  into  the  polarization  characteristics  of  the  Stokes  beam  relative  to  the  pump 
beam  in  Section  4.3.2.  In  this  section,  the  polarization  is  found  to  be  conserved  upon 
Brillouin  reflection  both  experimentally  and  theoretically.  Section  4.3.3  examines  the 
conditions  in  which  higher  order  spatial  modes  are  excited  by  SBS  in  the  fiber.  Finally, 
in  Section  4.3.4  beam  cleanup  using  a  ring  cavity  is  investigated. 

4.3.1.  50-micron  Beam  Cleanup 

Rodgers  demonstrated  SBS  beam  cleanup  at  -830  nm  using  an  8.2-micron  fiber  [9]. 
The  experiment  clearly  shows  the  viability  of  SBS  beam  cleanup.  However  it  was 
limited  by  the  fact  that  only  12  fiber  modes  could  exist  in  the  Corning  SMF-28  fiber 
used.  Because  of  this,  the  coupled  energy  in  the  fiber  closely  resembled  the  lowest  order 
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fiber  mode,  indieating  that  some  of  the  beam  cleanup  was  due  to  spatial  filtering  rather 
than  SBS.  In  order  to  increase  the  number  of  modes  supported  by  the  fiber  and 
investigate  the  Stokes  beam  response,  an  experiment  similar  to  Rodgers’s  was  performed 
using  a  50-micron  fiber  core.  For  such  a  large  core,  hundreds  of  fiber  modes  can  be 
excited  as  predicted  by  Equation  (2.17). 

Experimental  Setup 

The  experimental  setup  used  to  demonstrate  SBS  beam  cleanup  is  shown  in  Eigure  9. 
The  output  of  a  cw  Nd:YAG  Eightwave  Electronics  laser  operating  at  1.064  microns  was 
used  as  the  pump  beam.  The  laser  had  an  instantaneous  linewidth  less  than  5  kHz,  far 
less  than  the  SBS  gain  bandwidth,  ~100  MHz.  A  four-pass  amplifier  was  used  to 
increase  the  beam  power  to  greater  than  a  watt.  The  output  of  the  amplifier  then  passed 
through  a  microscope  slide  that  had  been  etched  for  five  minutes  in  hydrofluoric  acid. 

The  acid  created  random  pits  in  the  glass  surface  that  aberrated  the  phase  front  of  the 
pump  beam.  This  aberrated  beam  was  then  transmitted  through  a  Earaday  rotator,  an 
angled  plane-parallel  plate,  and  finally  into  the  4.4  km,  50-micron,  graded- index  fiber. 

The  backward  traveling  Stokes  beam  was  collimated  by  the  microscope  objective  and 
a  small  portion  of  the  beam  reflected  off  the  back  surface  of  the  plane  parallel  plate  and 
was  coupled  into  a  CCD  camera.  The  remaining  Stokes  beam  was  coupled  back  through 
the  Earaday  rotator  to  be  detected  by  a  Eabry-Perot  for  spectral  characterization  and  a 
power  meter  for  a  measurement  of  the  SBS  efficiency. 

The  CCD  images  were  used  to  characterize  the  improvement  in  the  quality  of  the 
Stokes  beam  over  that  of  the  pump  beam.  Eirst,  a  set  of  images  of  the  reflected  beam  was 
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Figure  9:  Schematic  of  experiment  used  to  examine  beam  cleanup. 


collected  above  SBS  threshold  that  included  the  Fresnel  reflection  off  the  front  fiber  facet 
as  well  as  the  Brillouin  reflected  beam.  Secondly,  without  adjusting  the  beam  coupling, 
the  fiber  was  cut  about  two  meters  beyond  the  front  facet  and  another  set  of  CCD  images 
was  collected.  The  reduction  in  fiber  length  increased  SBS  threshold  well  beyond  the 
power  generated  by  the  laser  and  amplifier,  thus  the  second  image  was  purely  due  to  the 
Fresnel  reflection  off  the  front  fiber  facet.  Finally,  a  third  group  of  images  was  collected 
of  the  transmitted  pump  beam  at  the  end  of  the  two-meter  section  of  fiber.  The  power 
coupled  into  the  fiber  was  also  measured  at  this  point  to  determine  the  coupling 
efficiency  into  the  fiber. 
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The  effect  of  the  location  of  the  pump  beam  on  the  position  of  the  Stokes  beam  was 
also  explored.  To  perform  this  experiment,  the  amplifier  and  phase  aberrator  were 
removed  from  the  experimental  setup  and  the  pump  beam  was  again  coupled  into  a  4.4 
km  fiber.  The  pump  beam  position  was  controlled  with  the  mirror  indicated  in  Figure  10 
At  each  mirror  position,  a  set  of  images  was  collected  showing  the  position  of  the  pump 
and  Stokes  beam  and  the  power  of  the  transmitted  pump  and  reflected  Stokes  beam  was 
recorded. 


10:  Schematic  depicting  the  experimental  setup  to  investigate  the  dependence  of 
pumping  conditions  on  the  efficiency. 


Results 

The  measured  power  in  the  Stokes  beam  as  a  function  of  the  power  coupled  into  the 
fiber  is  shown  in  Figure  1 1 .  The  solid  line  indicates  the  best  fit  to  the  data.  It  results  in  a 
threshold  of  420  mW  and  the  slope  efficiency  of  52.4%. 

The  intensity  distributions  of  the  reflected  pump,  Stokes,  and  transmitted  pump 
beams  are  shown  in  Figure  12(a-c)  respectively.  Part  (a)  clearly  shows  that  a  highly 
aberrated  pump  beam  is  incident  on  the  front  fiber  facet.  Part  (b)  demonstrates  that  the 
Stokes  beam  appears  to  be  excited  in  a  Gaussian-like  beam.  Finally,  part  (c)  reveals 
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that  the  pump  beam  is  coupled  into  many  fiber  modes,  ensuring  that  the  fiber  did  not 
spatially  filter  the  pump  beam. 


While  qualitatively  satisfying,  Figure  12  does  not  give  a  quantitative  measure  of  the 
beam  cleanup.  Thus,  in  order  to  complete  the  analysis,  a  measurement  of  beam  quality 
had  to  be  developed.  Typical  measurements  require  beam  stability  over  the  time 
frame  of  the  experiment.  This  allows  the  beam  width  to  be  measured  at  several  different 
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points  as  the  beam  passes  through  foeus.  By  eomparing  the  divergenee  of  the  beam  with 
that  of  a  perfeet  Gaussian  with  the  same  beam  waist,  the  beam  quality  ean  be  redueed  to  a 
single  number.  Unfortunately,  due  to  instabilities  in  the  experiment,  this  measurement 
eould  not  be  undertaken.  There  are  a  eouple  possible  reasons  the  laser  beam  was 
unstable.  First,  sinee  there  is  no  feedbaek  in  a  single  pass  SBS  setup,  the  stimulated 
proeess  must  begin  from  spontaneous  seattering,  whieh  is  inherently  random.  Seeondly, 
the  power  eoupled  into  the  fiber  was  measured  and  found  to  fluetuate  by  as  mueh  as  10%. 
This  eould  be  due  to  pointing  instability  in  the  laser  or  vibrations  of  the  fiber  end.  In 
either  ease,  the  measurement  of  was  replaeed  by  a  simpler,  but  less  robust  method 
deseribed  below. 


A  perfeet  Gaussian  beam  is  diffraetion  limited.  Thus,  how  well  a  beam  matehes  a 
Gaussian  intensity  distribution  is  a  measure  of  its  beam  quality.  This  is  the  prineiple 
behind  the  beam  quality  measurements  taken  in  this  experiment.  First,  the  area  under 
eaeh  intensity  distribution,  was  normalized  to  one.  With  the  normalization  eomplete,  the 
nonlinear  fit  algorithm  in  Mathematiea  was  used  to  fit  the  data  to  a  Gaussian  with  the 
same  area. 
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where  b  is  the  beam  size  and  is  the  beam  eenter.  This  algorithm  finds  a  loeal 

minimum  in  the  least  square  error  between  the  data  and  the  fitting  funetion.  Sinee  a  loeal 
minimum,  not  a  global  minimum,  is  found,  it  is  eritieal  that  a  good  estimate  be  used  as  a 
starting  point  for  the  seareh.  The  initial  guess  for  the  beam  eenter  was  estimated  by 
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finding  the  centroid  of  the  intensity  distribution.  The  beam  diameter  was  then  estimated 
by  calculating  the  standard  deviation  of  the  intensity  distribution  from  the  centroid.  If  the 
beam  was  a  perfect  Gaussian,  this  would  return  the  correct  beam  center  and  waist.  For 
beams  that  are  not  Gaussian,  it  is  a  good  starting  point  for  a  fit  to  the  data. 

The  best-fit  output  was  then  compared  to  the  actual  data  using  the  estimated  variance. 

In  this  equation  g.  j  is  the  value  of  the  Gaussian  fit  at  point  (/,  j) ,  /.  .  is  the  intensity  of 
the  beam  at  (i,  j) ,  and  -  3  is  the  number  of  degrees  of  freedom  in  the  Gaussian  fit. 

This  analysis  results  in  the  beam  quality  factors  summarized  in  Table  1.  These  values 
confirm  the  qualitative  analysis  of  Figure  12.  The  Stokes  beam  is  found  to  have  an 
intensity  distribution  that  matches  a  Gaussian  beam  1 8  times  better  than  the  pump  input 
and  more  than  9  times  better  than  the  transmitted  pump  beam. 


■''i''  1 :  Beam  quality  factors  for  beam  cleanup  experiment 
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The  effect  of  the  location  of  the  pump  beam  on  the  Stokes  output  was  also  monitored 
with  the  CCD  camera.  The  results  are  shown  in  Figure  13.  This  data  shows  that  as  the 
pump  beam  is  moved  from  the  upper  right  to  the  lower  right  of  the  frame,  the  output 
position  of  the  Stokes  beam  is  constant.  This  unequivocally  demonstrates  that  the  Stokes 
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beam  is  not  a  phase  conjugate  replica  of  the  pump  beam.  A  phase  conjugate  beam  would 
move  along  with  the  pump. 


Pump  Beam 


Stokes  Beam 

13:  The  Stokes  beam  does  not  signifieantly  move  as  the  input  beam  position  is 
varied  from  the  upper  right  to  the  lower  left. 


Although  not  as  obvious,  Figure  13  also  hints  that  the  efficiency  of  the  Stokes 
reflection  increases  when  the  pump  beam  is  coupled  into  the  center  of  the  fiber.  The 
power  of  the  reflected  Stokes  beam  is  constant  in  each  frame.  However,  its  intensity 
relative  to  that  of  the  Stokes  beam  clearly  changes.  The  contours  in  the  second  frame 
illustrate  that  the  intensity  of  the  Stokes  beam  is  much  greater  than  the  intensity  of  the 
reflected  pump  beam.  The  first  and  last  frames,  however,  show  that  the  reflected  pump 
beam  intensity  is  nearly  as  strong  as  the  Stokes  beam. 

This  was  further  investigated  by  measuring  the  power  that  was  transmitted  through 
and  reflected  from  the  fiber.  Figure  14(a)  shows  the  relative  transmission  coefficient  as  a 
function  of  the  pump  power  and  the  beam  position.  On  the  far  edges,  the  transmission 
drops  to  zero  because  the  beam  is  clipped  by  the  edge  of  the  fiber  core.  However,  even 
when  the  beam  is  fully  coupled  into  the  fiber,  the  transmission  above  SBS  threshold  is 
lower  near  the  center  of  the  core  than  near  the  edge  due  to  an  increase  in  SBS  efficiency. 
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to  the  other. 


The  transmission  curves  have  the  advantage  of  clearly  showing  the  physical  edges  of 
the  fiber.  However,  the  Stokes  beam  is  the  true  signal.  Figure  14(b)  shows  the  measured 
Stokes  signal,  which  demonstrates  the  same  features  as  the  transmitted  beam.  SBS  is 
found  to  be  most  efficient  near  the  center  of  the  fiber.  These  results  are  consistent  with 
the  theory  presented  in  Section  4.2.  The  effective  Brillouin  gain  is  proportional  to  the 
overlap  of  the  pump  and  Stokes  beams  as  given  by  Equation  (4.3).  Since  the  Stokes 
beam  is  consistently  found  in  the  LPoi  mode,  when  the  pump  beam  is  coupled  into  the 
liber  edges,  the  overlap  will  be  diminished  and  the  effective  Brillouin  gain  is  lowered. 

Of  course  even  a  spot  that  is  coupled  into  the  edge  of  the  core  will  evolve  into  a 
superposition  of  many  fiber  modes,  thus  the  overlap  of  the  LPoi  mode  and  the  pump 
beam  at  the  front  fiber  facet  will  not  necessarily  approximate  the  overlap  of  the  two 
beams  along  the  entire  fiber  length. 

On  the  other  hand,  the  spatial  extent  in  which  SBS  is  observed  in  Figure  14  is  a  good 
measure  of  where  the  pump  beam  can  be  coupled  into  the  fiber,  while  still  contributing  to 
the  overall  Stokes  beam.  To  determine  this  spatial  extent,  the  efficiency  of  the  Brillouin 
reflection  was  assumed  to  be  proportional  to  a  Gaussian  function  concentric  with  the 
fiber  core. 
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where  rf  is  the  slope  effieieney  when  the  pump  beam  is  coupled  into  the  center  of  the 
fiber  core  and  \  is  the  radius  of  the  Gaussian  function.  The  efficiency  affects  the  overall 
Stokes  power  in  the  following  manner 


P.= 


0  ri{r)Pp<P,h 
l{r)Pp-P,h  l{r)Pp>P,h' 


(4.16) 


This,  in  turn,  can  be  used  to  estimate  the  transmission  coefficient  at  various  pump  power 
levels  and  to  find  . 


Figure  15  shows  several  slices  of  the  transmission  surface  at  various  power  levels 
along  with  the  best  fit  to  them.  In  part  (a),  the  transmission  below  threshold  is  shown. 
This  was  used  to  estimate  the  size  of  the  pump  beam  on  the  front  fiber  facet.  The  plots 
shown  in  (b,c)  demonstrate  the  transmission  above  threshold  at  two  power  levels.  The 
best  fits  were  found  with  different  values  for  .  In  (d),  the  size  of  b^  is  plotted  as  a 

function  of  the  pump  power.  Clearly,  this  is  not  a  constant  value  and  as  the  pump  power 
increases,  the  edges  of  the  fiber  become  increasingly  valuable. 

These  results  indicate  that  the  efficiency  of  SBS  beam  cleanup  will  be  a  function  of 
the  pump  beam  spatial  quality.  For  low  beam  quality,  when  much  power  is  coupled  into 
the  edges  of  the  fiber,  the  conversion  efficiency  will  be  low  and  threshold  will  be  high. 
Once  above  threshold  however,  the  entire  pump  beam  will  contribute  to  the  Stokes  beam. 
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This  is  unlike  a  spatial  filter,  whieh  bloeks  the  high  spatial  frequeney  eomponents  of  the 
beam  and  thus  has  a  very  low  efficieney  for  highly  aberrated  beams. 
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;  15:  Transmission  of  the  pump  beam  through  a  50-mieron  fiber  versus  pump  beam 
position  at  various  pump  power  levels,  (a)  Low  power  transmission  below  SBS 
threshold,  (b-e)  Transmission  above  threshold,  (d)  Effeetive  cleanup  diameter  as 
a  function  of  the  pump  power. 


4.3.2.  SBS  Polarization  Properties 

Rodgers  investigated  the  polarization  properties  of  the  Stokes  beam  generated 
through  stimulated  Brillouin  scattering  and  found  that  it  had  the  same  polarization  as  the 
pump  beam.  This  differs  from  a  standard  mirror  in  which  a  right-hand-circularly 
polarized  incident  wave  is  reflected  into  a  left-hand-circularly  polarized  beam.  It  also 
contradicts  the  prediction  of  ZeTdovich  [39],  and  therefore  was  investigated  further. 


Experimental  Setup 

Figure  16  shows  a  schematic  of  the  experimental  setup  used  to  investigate  the 
polarization  properties  of  the  Stokes  beam.  A  diode  laser  operating  above  SBS  threshold 
at  83 1  nm  was  transmitted  through  a  quarter- wave  plate  thus  converting  the  horizontally 
polarized  pump  beam  to  a  right-hand-circularly  polarized  beam.  The  light  was  then 
coupled  into  an  8.2-micron  step-index  fiber.  The  reflected  Stokes  beam  traveled  back 
through  the  quarter-wave-plate.  After  this  second  pass  through  the  quarter-wave-plate 
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portion  of  the  reflected  beam  that  was  right  (left)-hand-circularly  polarized  became 
vertically  (horizontally)  polarized.  A  portion  of  this  beam  was  then  reflected  off  an 
uncoated  plane -parallel  plate,  transmitted  through  a  Glan-Thompson  polarization 
analyzer  and  detected  by  a  Fabry-Perot.  The  relative  power  in  the  reflected  pump  and 
Stokes  beam  was  measured  by  the  Fabry-Perot  trace  and  recorded  at  various  analyzer 
angles. 
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Figure  16:  Schematic  diagram  of  polarization  characterization  experiment. 


Results 

Figure  17  shows  the  results  of  this  measurement.  The  power  in  the  reflected  pump 
beam,  represented  by  the  solid  line,  is  caused  by  the  Fresnel  reflection  off  the  front  fiber 
facet.  Fresnel  reflection  is  well  known  to  reverse  the  handedness  of  the  polarization. 
Conversely,  the  Stokes  beam,  represented  by  the  dashed  line,  is  found  to  have  a 
polarization  that  is  orthogonal  to  the  pump,  thus  showing  that  the  handedness  of  the 
polarization  is  conserved  when  the  beam  is  reflected  through  stimulated  Brillouin 
scattering. 
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-  ;  17:  Reflected  pump  (solid)  and  Stokes  (dashed)  power  as  a  tunction  of  the 

-  yzer  angle.  The  two  are  out  of  phase  by  90°  showing  that  the  polarization  of  the 

Stokes  beam  is  orthogonal  to  that  of  the  reflected  pump. 


Although  the  observed  polarization  properties  of  the  Stokes  beam  were  not  expected 
[39],  they  are  easily  explained  by  decomposing  the  pump  and  Stokes  electric  fields  into 
X  and  y  components  and  investigating  how  the  Brillouin  reflection  affects  the  phase  of 
each. 


Ep  =  Ap,  exp \i[k^z  - co^t  +  )] i  +  exp [/ [k 

Es  =  4.  exp  [/ {-k^z  -  aj  +  )]  i  +  4^  exp  [/ (-k^ 


-kz-(X>t  + 


(4.17) 


^px,sx  (  ^py,sy  )  thc  amplitude  of  the  pump  and  Stokes  fields  along  the  x-{y  -)axis  and 

^px  sx  ( ^px  sx )  is  iiie  phase  of  the  x  ( y )  component  of  the  field.  Assuming  the  Stokes 

beam  is  not  externally  seeded,  it  is  initialized  by  spontaneous  scattering  of  the  pump  off 
thermal  phonons  in  the  fiber  [40].  The  phonon  frequency  and  wavenumber  are  fixed  by 
material  properties  thus  the  density  fluctuations  in  the  fiber  is  assumed  to  have  the  form 
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(4.18) 


Pn  =  A  exp 


(j)^  is  the  random  phase  of  the  phonon  and  is  the  amplitude  of  the  wave. 


Stimulated  seattering  requires  that  the  reflected  Stokes  beam  interfere  with  the  pump 
beam  such  that  the  acoustic  wave  is  intensified  through  electrostriction.  In  other  words, 
the  high  intensity  regions  caused  by  the  interference  between  the  pump  and  Stokes  beams 
must  overlap  the  high  density  regions  in  the  fiber. 


exp {/[(k^  +  kj z  - t  +  (j)^^  -  (f)^^ ]| 
exp +  kj z  - t  +  (j)^^  - (j)^^ ]| 


=  exp{/[(k^+k^)z-(®^-(yJt  +  ^„]| 
=  exp{/[(k,+k^)z-(®^-®Jt  +  ^/i„]| 


(4.19) 


Both  of  these  equations  are  only  satisfied  if  the  difference  between  the  phases  goes  to 
zero.  Since  the  phase  of  the  thermal  phonon  and  both  components  of  the  pump  beam  are 
fixed,  the  phase  of  both  components  of  the  Stokes  beam  can  be  found. 


A=(l>p.-(l>n 

<t>sy=<t>py-<t>n 


(4.20) 


Thus  the  y  -components  of  both  the  pump  and  Stokes  beams  lag  the  x  -components  by  an 
equal  amount,  -  (j)^^  =  -  (j)^^ ,  and  the  handedness  of  the  polarization  is  preserved. 


Although  the  polarization  is  maintained  during  Brillouin  reflection  at  each  point  in 
the  fiber,  the  polarization  is  not  preserved  as  the  beam  propagates  through  the  fiber. 
Fiber  imperfections  result  in  a  small,  random  birefringence  that  leads  to  a  fluctuating 
elliptical  polarization  along  the  fiber  length.  The  birefringence  is  not  large, 
approximately  1  part  in  10  million  [41],  but  the  long  fiber  lengths  used  in  beam  cleanup 
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make  it  significant.  If  the  birefringence  was  oriented  in  a  constant  direction,  the  fiber 
could  rotate  the  polarization  of  a  guided  beam  by  90°  in  approximately  5  m.  In  practice, 
however,  fibers  that  are  thousands  of  meters  long  reflect  Stokes  beams  that  are  identical 
in  polarization  to  the  pump  beam.  If  the  Stokes  beam  propagated  at  the  same  frequency 
and  in  an  identical  fiber  mode,  it  would  pass  back  through  the  effective  waveplate 
produced  by  the  fiber  birefringence,  perfectly  retracing  the  polarization  of  the  pump 
beam.  The  Stokes  beam,  however,  is  not  at  the  same  frequency  as  the  pump  due  to  the 
Brillouin  shift  and  cleanup  modifies  the  beam’s  spatial  structure.  These  characteristics  of 
the  Stokes  beam  limit  the  useful  fiber  length  if  a  polarized  output  is  required. 

Since  the  birefringence  is  random,  the  relative  phase  shift  imparted  to  the  orthogonal 
polarization  modes  does  not  grow  linearly.  The  phase  difference  can  be  treated  like  a 
random  walk.  When  the  birefringence  is  aligned  such  that  the  fast  axis  coincides  with  the 
X  -direction,  the  random  step  advances  the  phase  of  the  x  -polarized  beam.  When  the 
birefringence  aligns  the  fast  axis  in  the  y  -direction,  the  random  step  advances  the  phase 
of  the  y  -polarized  beam.  For  simplicity,  the  birefringence  is  constrained  to  one  of  these 
two  axes.  Although  a  more  general  approach  would  include  random  steps  where  the  fast 
axis  is  aligned  at  random  angles  to  the  x  -  and  y  -axes,  the  approach  taken  shows  the 
features  of  polarization  mode  dispersion  needed.  The  birefringence  is  assumed  to  be 
constant  over  lengths,  ,  making  discontinuous  shifts  between  segments.  The  distance 
of  each  step  in  this  random  walk  is  the  magnitude  of  the  phase  shift  imparted  to  the  beam 
over  =  {^’y'2)Li,^n .  The  number  of  steps  taken  is  the  number  of  constant 

birefringence  zones  in  the  fiber,  N  =  ,  where  is  the  fiber  length. 
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The  probability  that  m  of  the  N  zones  are  oriented  such  that  light  polarized  along 
the  X  -axis  propagates  faster  than  that  polarized  along  the  y  -axis  is  given  by  the  binomial 
probability  distribution,  which  for  large  N  can  be  approximated  by  a  Gaussian 
distribution, 
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(4.21) 


If  a  large  number  of  fibers  were  tested,  it  would  be  found  that  the  average  number  of 
zones  with  the  fast  axis  oriented  in  the  x  -direction  is 
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(4.22) 


As  expected,  this  is  exactly  14  of  the  total  number  of  zones.  This  is  required  by  the  fact 
that  the  x  -  and  y  -axes  have  not  yet  been  defined  so  there  cannot  be  a  preferred  axis.  On 
the  other  hand,  this  is  only  an  average.  Some  of  the  fibers  would  be  found  to  have  a  net 
birefringence  along  the  x  -axis  and  some  with  a  net  birefringence  along  the  y  -axis.  The 
difference  between  the  number  of  zones  oriented  along  the  x  -axis,  m  ,  and  the  number 
oriented  along  the  y  -axis,  N  -m  ,is  the  number  of  zones  that  constructively  add  to 
create  the  total  birefringence  of  a  fiber,  b^=N  -  2m  .  The  average  magnitude  of  this 
number  is  proportional  to  the  average  magnitude  of  the  birefringence  of  a  fiber. 
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Thus  the  number  of  zones  that  eonstructively  add  together  is  proportional  to  the  square 
root  of  the  total  number  of  zones  and  the  average  relative  phase  shift  at  the  end  of  the 
fiber  is 


Afi,=Afi =  .  (4.24) 

This  shows  that  the  phase  differenee  grows  in  proportion  to  the  square  root  of  the  liber 

length  rather  than  linearly  with  the  length.  The  eolleetion  of  terms  ,  where  c 

is  the  speed  of  light,  is  known  as  the  eharacteristic  polarization  mode  dispersion  (PMD). 

It  is  an  important  measure  of  the  pulse  spreading  in  fiber  waveguides  beeause  a  high 
eharaeteristie  PMD  will  limit  the  rate  that  data  ean  be  sent  over  fiber  optic  cables.  It  is 
typically  on  the  order  of  0.1  [41]  and  is  introduced  here  because  it  can  be  used  to 

estimate  the  value  of  . 

As  the  pump  beam  travels  down  the  fiber,  the  phase  delay  of  orthogonal  polarization 
components  differ  by  an  amount  given  by  Equation  (4.24).  The  Stokes  beam  will 
experience  similar  phase  delays,  but  since  it  is  propagating  backward  through  the  fiber 
and  starts  with  the  polarization  that  is  conjugate  to  the  pump,  the  phase  delays  act  to 
reconstruct  the  original  pump  polarization.  Any  differences  between  the  phase  delay 
imparted  to  the  pump  and  the  phase  delay  imparted  to  the  Stokes  beam  will  generate  a 


55 


polarization  shift  between  the  input  pump  and  the  output  Stokes.  Setting  the  differenee 


between  the  two  phase  delays  equal  to  n  gives  an  estimate  to  the  allowed  fiber  length  that 
eonserves  polarization  during  Brillouin  refieetion. 


\  J 


=  n 
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Chromatie  dispersion  and  modal  dispersion  affeet  the  birefringenee  in  the  fiber. 
Chromatie  dispersion  is  relatively  small  over  the  wavelength  shift  indueed  by  SBS  and 
will  be  ignored.  Modal  dispersion,  however,  is  signifieant  since  beam  cleanup  requires 
that  the  Stokes  beam  propagate  in  the  LPoi  mode  and  the  pump  beam  to  propagate  in 
higher  order  modes.  Figure  18  shows  the  allowed  fiber  length  as  a  function  of  the  mode 
number  for  a  50-micron  step  index  fiber.  Clearly  the  higher  order  modes  require  a 
shorter  fiber  if  polarization  is  to  be  maintained.  However,  even  if  the  pump  is  coupled 
into  the  highest  order  mode  and  the  Stokes  beam  propagates  in  the  LPoi  mode,  the 
polarization  can  be  maintained  in  fibers  that  are  nearly  3  km. 

There  are  several  limitations  to  this  calculation.  First,  it  is  not  likely  that  the  pump 
beam  will  propagate  in  a  pure  fiber  mode,  so  the  actual  allowed  fiber  length  is  expected 
to  be  a  weighted  average  that  takes  into  account  the  percentage  of  pump  power  in  each 
fiber  mode.  Secondly,  the  50-micron  fiber  used  in  the  experiment  had  a  graded  index 
profile.  Graded-index  fibers  are  designed  to  have  much  less  modal  dispersion,  thus  it  is 
expected  that  the  fiber  length  could  be  much  greater  and  still  maintain  the  Stokes 
polarization.  Additionally,  the  calculation  ignores  the  fact  that  most  of  the  Brillouin  gain 
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occurs  in  the  front  end  of  the  fiber  so  it  is  more  appropriate  to  interpret  the  vertical  axes 
of  Figure  18  as  the  effective  length  rather  than  the  physieal  length  of  the  fiber. 


mode.  The  Stokes  beam  is  assumed  to  propagate  in  the  LPqi  mode. 


4.3.3.  High  Spatial  Order  Stokes  Modes 

In  the  preeeding  seetions  of  Chapter  4,  SBS  beam  cleanup  is  experimentally  shown  to 
produce  a  high  quality  Stokes  beam  in  the  LPoi  fiber  mode  under  several  different 
pumping  eonditions.  Flowever,  the  theoretieal  examination  of  beam  eleanup 
demonstrated  that  it  should  be  possible  to  exeite  other  pure  fiber  modes.  This  section 
describes  experimental  evidenee  of  this,  as  well  as  confirmation  of  mode  competition  in 
the  fiber. 


Experimental  Setup 

A  sehematie  of  the  experimental  setup  is  shown  in  Figure  19.  The  frequeneies  of  two 
Lightwave  Eleetronics  Nd;YAG  single-frequeney  lasers  were  tuned  to  within  the 
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Brillouin  gain  bandwidth  of  each  other.  The  beams  were  then  spatially  combined  using  a 
polarization  beamsplitter  and  coupled  into  a  4.4  km,  50-micron,  graded-index  fiber.  A 
portion  of  the  reflected  Stokes  beam  was  picked  off  using  the  Fresnel  reflection  off  the 
front  surface  of  an  uncoated  plane-parallel  plate.  It  was  then  viewed  with  a  CCD  camera. 
The  portion  of  the  beam  that  reflected  off  the  back  surface  of  the  plane-parallel  plate  was 
monitored  with  a  scanning  Fabry-Perot  interferometer  to  ensure  the  frequencies  of  the 
two  beams  did  not  drift  away  from  each  other. 


Optical 


Figure  19:  Multimode  Stokes  beam  experimental  schematic. 


Laser  2 
LWE  700  mW 


Laser  1 

LWE  700  mW 


The  coupling  of  the  two  beams  into  the  fiber  was  carefully  controlled.  With  the  first 
laser  blocked,  the  position  of  the  second  laser  beam  was  adjusted  to  maximize  the  Stokes 
power  generated.  Once  completed,  the  second  laser  was  blocked  and  the  coupling  of 
laser  1  was  adjusted  until  the  reflected  Stokes  beam  appeared  to  propagate  in  the  LPn 
mode.  It  was  not  a  simple  task  to  excite  this  spatial  mode.  Once  excited,  it  would  often 
spontaneously  shift  into  the  LPoi  intensity  pattern  and  again  have  to  be  adjusted  to  excite 


58 


the  LPii  mode.  CCD  images  of  the  Stokes  beam  excited  by  each  laser  operating 
independently,  as  well  as  both  lasers  working  together  were  collected.  Finally,  the 
frequency  of  the  second  laser  was  tuned  slightly  away  from  that  of  laser  1 .  Again  the 
intensity  of  the  Stokes  beam  generated  by  the  combined  action  of  both  beams  was  imaged 
on  the  CCD. 

Results 

The  Stokes  beams  generated  when  laser  1  and  laser  2  each  independently  pumped  the 
fiber  are  shown  in  Figure  20(a)  and  (b).  Laser  1  clearly  produced  a  two  lobed  output, 
unlike  the  beams  generated  through  SBS  with  a  highly  aberrated  beam.  This  output  is 
similar  to  the  LPn  optical  mode  in  a  fiber.  On  the  other  hand,  laser  2  generated  a  Stokes 
beam  with  a  Gaussian  intensity  distribution  consistent  with  the  LPoi  mode. 

Parts  (c)  and  (d)  of  Figure  20  show  the  intensity  distribution  of  the  Stokes  beam 
produced  when  the  fiber  was  pumped  by  both  lasers  simultaneously.  The  difference 
between  the  figures  results  from  the  spectral  characteristics  of  the  pump  beams.  Part  (c) 
shows  the  output  when  both  lasers  operate  at  the  same  frequency.  It  will  be  demonstrated 
in  Chapter  5  that  under  this  condition,  the  two  lasers  work  together  to  produce  a  single 
Stokes  beam  that  is  temporally  coherent.  In  contrast,  when  the  two  laser  frequencies  are 
separated  by  more  than  the  Brillouin  gain  bandwidth,  they  generate  independent  Stokes 
beams  and  do  not  interact  in  a  significant  way.  This  condition  is  shown  in  part  (d)  of 
Figure  20. 

When  the  two  pump  beams  work  together  (see  Figure  20(c)),  the  LPn  mode 
generated  by  laser  1  is  lost  and  a  single  Gaussian-like  Stokes  beam  is  generated.  In  order 
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• 

(a)  (b) 

• 

(c) 

(d) 

I;  20:  (a)  The  intensity  pattern  of  the  Stokes  beam  generated  by  laser  1.  (b)  The 
tensity  pattern  of  the  Stokes  beam  generated  by  laser  2.  (c)  The  intensity  pattern 
of  the  Stokes  beam  generated  when  both  lasers  pump  the  fiber  simultaneously  at  the 
same  frequeney.  (d)  The  intensity  pattern  of  the  Stokes  beam  generated  when  both 
lasers  pump  the  fiber  simultaneously  at  different  frequeneies. 


to  prove  that  the  LPn  mode  is  not  just  hidden  under  a  mueh  stronger  LPqi  mode,  the 
powers  of  the  two  independently  generated  Stokes  beams  were  measured  and  found  to  be 
nearly  the  same.  Thus,  if  these  two  beams  were  incoherently  added,  the  Stokes  intensity 
pattern  would  still  contain  a  strong  portion  of  the  two  lobes  from  laser  1 . 

It  could  also  be  argued  that  the  addition  of  laser  2  changed  the  optical  path  of  the  first 
laser  by  a  tiny  fraction.  This  could  be  a  result  of  heating  the  fiber,  a  change  in  path  due  to 
the  second  acoustic  wave  produced,  or  another  nonlinear  effect.  Since  the  generation  of 
the  higher  order  Stokes  beam  in  the  LPn  mode  was  extremely  sensitive  to  the  pump 
alignment,  any  of  these  minuscule  effects  could  change  the  coupling  enough  to  cause  the 
first  Stokes  beam  to  generate  an  LPoi  mode.  On  the  other  hand,  none  of  these  should  be 
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significantly  affected  by  a  shift  in  the  frequency  of  laser  2  by  less  than  IGHz.  Yet  when 
the  second  laser  was  tuned  outside  the  Brillouin  gain  bandwidth,  the  Stokes  beam 
intensity  distribution  reverted  to  that  shown  in  part  (d),  which  closely  matches  the 
incoherent  addition  of  the  intensity  patterns  shown  in  parts  (a)  and  (b). 

The  most  viable  explanation  to  the  observed  phenomena  is  that  the  spatial  modes  in 
the  fiber  eompete  for  the  gain.  When  laser  1  pumped  the  fiber  alone,  the  LPn  mode  had 
the  highest  gain,  resulting  in  an  output  that  closely  matched  it.  However,  when  both 
lasers  worked  together,  the  LPoi  mode  had  the  greatest  gain  and  beeame  dominant,  thus 
suppressing  the  higher  order  modes. 

4.3.4.  Ring  Cavity  Beam  Cleanup 

The  SBS  threshold  is  reduced  when  the  Stokes  beam  is  seeded  at  the  back  end  of  the 
fiber.  This  deereases  the  overall  loss,  thus  increasing  effieiency.  While  a  second  laser  at 
the  Stokes  frequeney  is  the  most  obvious  choiee,  in  practiee  it  is  difficult  to  control  the 
two  lasers  to  ensure  the  seed  beam  is  pumped  within  the  narrow  Brillouin  gain  bandwidth 
of  the  pump  beam.  The  approach  taken  by  Choi  in  his  research  at  AFIT  was  to  use  a 
portion  of  the  pump  to  excite  a  Stokes  beam  in  a  second  fiber.  This  beam  was  then 
coupled  into  the  back  end  of  the  fiber  used  for  beam  cleanup  [42].  In  this  geometry,  a 
single-mode  fiber  could  be  used  to  generate  the  seed  beam,  thus  threshold  for  the  seed  is 
very  low.  This  configuration  also  has  the  benefit  of  ensuring  that  as  the  frequency  of  the 
pump  beam  drifts,  the  seed  frequency  follows  it.  Unfortunately,  a  highly  aberrated  pump 
beam  will  not  effieiently  eouple  into  a  single  mode  fiber,  increasing  the  power  required  to 
generate  the  seed  and  reducing  the  overall  effieieney. 
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Another  approach  that  eliminates  these  problems  is  to  use  a  single  fiber  in  a  ring 
cavity.  This  has  the  benefit  of  reducing  the  SBS  threshold  [27],  but  does  not  waste  power 
due  to  low  coupling  efficiency.  This  approach  was  taken  in  the  research  described  here. 

Experimental  Setup 

The  experiment  was  completed  in  two  parts.  First,  the  Stokes  beam  shape  was 
observed  to  determine  if  seeded  beam  cleanup  is  possible.  Secondly,  the  Stokes  beam 
power  was  measured  at  different  levels  of  seeding  to  determine  how  much  was  required 
to  increase  SBS  efficiency.  These  experiments  used  the  setup  shown  in  Figure  21. 
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As  in  the  previous  beam  cleanup  experiment,  a  Lightwave  Electronics  700  mW 
source  was  aberrated  using  an  etched  glass  plate.  The  beam  then  traveled  through  a 
polarization  beamsplitter,  Faraday  rotator,  and  was  focused  into  a  4.4  km,  50-micron 
optical  fiber.  After  scattering  out  the  front  of  the  liber,  a  portion  of  the  beam  was 
reflected  off  the  front  and  back  surfaces  of  a  plane  parallel  plate.  The  reflection  off  the 
back  surface  was  imaged  with  a  CCD  camera.  The  reflection  off  the  front  surface  was 
measured  with  a  power  meter.  The  rest  of  the  beam  passed  back  through  the  Faraday 
rotator  and  was  transmitted  through  the  polarization  beam  splitter.  Another  half-wave 
plate  and  polarizer  were  used  to  control  the  amount  of  the  reflected  Stokes  beam  that  was 
coupled  into  the  back  end  of  the  fiber.  The  spectrum  of  the  transmitted  and  reflected 
beams  was  measured  using  the  reflection  off  a  microscope  slide  and  a  scanning  Fabry- 
Perot  interferometer. 

Results 

Figure  22  shows  the  intensity  distribution  of  the  Stokes  beam  for  a  variety  of  half¬ 
wave-plate  settings.  The  labels  indicate  the  ratio  of  the  power  that  is  reflected  by  the 
second  polarization  beam  splitter  to  the  power  incident  upon  it.  It  is  related  to  the  quality 
factor,  Q ,  of  the  ring  cavity,  and  is  directly  proportional  to  the  seed  power  coupled  into 
the  back  end  of  the  fiber.  The  proportionality  constant  is  determined  by  the  transmission 
of  all  the  optical  elements  as  well  as  the  coupling  efficiency  at  the  back  fiber  facet.  This 
coupling  efficiency  is  expected  to  be  low.  A  40x  microscope  objective  with  a  numerical 
aperture  of  0.4  was  used  at  this  end  of  the  fiber  to  ensure  the  seed  beam  filled  fiber  core 
( NA  =  0.2  ),  regardless  of  its  mode  structure. 
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22:  Stokes  beam  shapes  at  various  levels  of  Stokes  eoupling  into  the  baek  end  of 

the  fiber. 


For  high  cavity  quality,  the  ring  cavity  experiment  verified  the  results  of  Choi  [42] 
and  beam  cleanup  was  not  observed.  The  beam  quality,  as  given  by  Equation  (4.14),  is 
tabulated  below.  Both  the  intensity  distribution  and  the  tabulated  data  suggest  an 
improvement  in  the  quality  of  the  beam  as  the  amount  of  Stokes  seeding  is  reduced. 


■»2:  Beam  quality  factors  for  ring  cavity  SBS  experiment. 


SBS  Feedbal^ 

100% 

88% 

59% 

25% 

3% 

0% 

Variance  of  Gaussian 
Fit(10-')-y 

113±5 

62  ±20 

50±17 

17±7 

6.1F0.4 

4.3F0.3 

Relative  Fit 

26.3 

14.4 

11.6 

4.0 

1.4 

1 

When  very  little  of  the  Stokes  beam  was  coupled  back  into  the  fiber,  the  beam  spot 
appeared  to  have  a  Gaussian  intensity  distribution,  yet  the  benefit  of  that  slight  seeding  to 
the  overall  Stokes  power  is  not  readily  measured  from  this  data.  To  quantify  the  benefit 
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of  the  seed  beam,  the  power  in  the  Stokes  beam  was  measured  as  a  function  of  the  half¬ 
wave-plate  angle.  The  experiment  was  set  up  in  response  to  the  images  observed  in 
Figure  22  but  after  the  original  experiment  had  been  disassembled,  so  the  specific  images 
of  the  Stokes  beam  cannot  be  directly  correlated  to  the  effect  on  the  Stokes  power.  On 
the  other  hand,  the  experimental  schematic  was  virtually  identical  so  there  is  no  reason  to 
expect  that  the  Stokes  spatial  quality  would  be  significantly  different.  The  result  of  this 
measurement  is  plotted  in  Figure  23.  This  data  shows  that  even  a  relatively  small  Stokes 
seed  results  in  a  drastic  increase  in  the  Stokes  output,  but  it  contributes  no  information 
about  whether  the  improvement  is  due  to  an  increase  in  slope  efficiency  or  a  reduction  in 
threshold.  Comparing  it  to  the  spatial  quality  of  the  Stokes  beam  indicates  that  it  is  likely 
that  most  of  the  benefit  of  a  seed  could  be  obtained  while  maintaining  a  relatively  clean 
Stokes  beam. 
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23 :  Stokes  output  at  various  levels  of  Stokes  coupling  into  the  back  end  of  the 

fiber. 
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In  the  experiments  deseribed  above,  the  Stokes  beam  was  deliberately  eoupled  into 
multiple  fiber  modes.  A  representative  sample  of  the  coupled  seed  beam  is  shown  in 
Figure  24.  Since  it  is  unlikely  that  the  power  is  uniformly  distributed  over  all  of  the  fiber 
modes,  this  preferentially  enhances  certain  modes,  invalidating  the  beam  cleanup 
modeling  described  in  Section  4.2.  However,  if  the  ring  cavity  was  designed  such  that 
the  LPoi  mode  at  the  front  end  of  the  fiber  was  re-imaged  at  the  back  end  of  the  fiber, 
beam  cleanup  should  be  observed  regardless  of  how  much  of  the  Stokes  beam  is  coupled 
back  into  the  fiber.  Ensuring  the  Stokes  is  injected  into  the  fiber  in  the  pure  LPqi  mode 
can  be  done  using  a  number  of  techniques  [36,  43]. 


24:  Beam  shape  of  the  eoupled  Stokes  power  in  the  ring  cavity  setup. 


4,4,  Conclusions 

Throughout  this  chapter,  beam  cleanup  in  a  long,  50-micron,  graded  index  fiber  was 
analyzed.  As  described  in  Section  4.3.1  the  spatial  intensity  distribution  was  clearly 
altered  when  an  aberrated  pump  beam  was  coupled  into  a  fiber.  The  beam  quality  was 
quantitatively  analyzed,  confirming  that  the  Stokes  beam  quality  was  a  factor  of  18.2 
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better  than  the  incident  pump  and  9.1  times  better  than  the  coupled  pump  beam.  The 
beams  were  also  analyzed  qualitatively.  The  reflected  pump  and  transmitted  pump 
beams  were  highly  aberrated  with  clear  intensity  fluctuations  across  the  image,  whereas 
the  Stokes  beam  had  a  smooth  Gaussian-like  intensity  distribution. 

This  section  also  investigated  the  efflciency  of  the  beam  cleanup  process.  It  was 
found  that  the  efflciency  drops  dramatically  when  the  pump  beam  is  not  coupled  directly 
into  the  center  of  the  fiber.  Thus  a  highly  aberrated  pump  beam  is  not  converted  into  a 
Stokes  beam  as  efficiently  as  one  that  perfectly  matched  the  LPoi  fiber  mode.  The  results 
indicate,  however,  that  the  region  near  the  center  of  the  fiber  that  must  be  pumped  to 
efficiently  excite  a  Stokes  beam,  is  a  function  of  the  pump  power.  Thus  for  aberrated 
pump  beams  significantly  above  SBS  threshold,  the  efficiency  of  beam  cleanup  is 
expected  to  be  greater  than  that  of  a  spatial  filter,  which  blocks  the  aberrated  portion  of 
the  pump  beam. 

In  Section  4.3.2,  the  polarization  properties  of  the  Stokes  beam  were  analyzed.  In 
contrast  to  the  reflection  properties  of  a  planar  mirror,  it  was  experimentally  shown  that 
the  Stokes  beam  propagates  in  the  same  polarization  state  as  the  pump  beam.  This  was 
also  modeled  theoretically.  By  assuming  both  the  x  -  and  y  -components  work  together 
to  generate  a  single  density  perturbation  in  the  fiber,  the  observed  polarization  properties 
of  the  Stokes  beam  were  predicted.  Additionally,  the  model  allowed  for  an  estimate  of 
the  maximum  fiber  length  that  would  be  expected  to  maintain  the  same  polarization 
characteristics.  For  a  pump  beam  propagating  in  the  highest  order  fiber  mode  and  the 
Stokes  beam  propagating  in  the  LPoi  mode,  the  Stokes  polarization  could  be  maintained 
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for  nearly  3km.  However,  the  approximations  made  in  the  model  underestimate  the 
aetual  allowed  fiber  length,  thus  the  Stokes  beam  is  expeeted  to  maintain  the  polarization 
of  the  pump  beam  in  libers  that  are  in  fact  much  longer. 

The  modeling  described  in  Section  4.2  demonstrated  that  it  should  be  possible  to 
generate  a  Stokes  beam  in  higher  order  liber  modes  through  careful  pumping  conditions. 
The  results  presented  in  Section  4.3.3  show  that,  although  infrequent,  this  is  true.  A 
Stokes  beam  propagating  in  the  LPn  mode  was  observed.  When  a  second  laser 
simultaneously  pumped  the  fiber  at  the  same  frequency,  the  LPn  mode  was  lost  and 
power  from  both  lasers  was  observed  in  the  LPoi  mode.  This  provided  strong  evidence 
that  mode  competition  occurs  in  the  liber. 

Finally,  the  feasibility  of  using  a  Stokes  seed  to  increase  the  efficiency  of  SBS  beam 
cleanup  was  examined  in  Section  4.3.4.  Feeding  a  portion  of  the  Stokes  beam  back  into 
the  rear  end  of  the  liber  provided  the  seed.  When  the  Stokes  feedback  was  strong,  the 
beam  cleanup  characteristics  of  SBS  were  lost.  In  contrast,  when  the  feedback  was  kept 
low,  it  appeared  that  the  SBS  efficiency  could  be  increased  while  still  cleaning  the 
aberrated  pump  beam. 
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5.  SBS  Beam  Combining 


5.1.  Background 

SBS  beam  combining  has  been  investigated  for  many  years.  As  in  the  case  of 
beam  cleanup,  however,  these  approaches  typically  require  a  high  quality  source  that  can 
be  amplified.  In  one  scenario,  a  single  beam  is  split  into  several  arms.  Each  arm  is 
separately  amplified,  phase  conjugated,  sent  back  through  the  amplifier,  and  recombined. 
In  this  arrangement,  each  Stokes  beam  starts  from  a  different  noise  photon  thus  the  beams 
are  not  temporally  coherent  with  respect  to  one  another.  To  eliminate  this  effect,  the 
beams  can  be  conjugated  in  the  same  SBS  medium  so  that  they  grow  from  the  same  seed 
photon  [4-8]. 

Another  method  to  ensure  that  each  beam  is  coherent  with  respect  to  the  others  is  to 
generate  a  Stokes  beam  that  can  be  used  to  seed  each  arm  of  the  setup.  This  way  the 
different  arms  do  not  rely  upon  noise  photons  to  initialize  process  and  each  arm  will  be 
coherent  with  respect  to  the  others.  This  has  the  added  benefit  of  reducing  the  threshold 
for  SBS.  Unfortunately,  the  phase  conjugation  fidelity  is  reduced  if  the  seed  beam 
becomes  too  powerful  [44]. 

Rodgers  took  a  very  different  approach  to  beam  combining  that  is  akin  to  the  beam 
cleanup  found  in  the  preceding  chapter.  It  is  this  approach  that  forms  the  basis  of  the 
research  examined  in  this  chapter.  Since  a  single  pump  beam  produces  a  Stokes  beam  in 
the  LPoi  mode,  two  pump  beams  in  the  same  fiber  should  produce  two  Stokes  beams, 
both  propagating  in  the  LPoi  mode.  Thus  two  spatially  distinct  pump  beams  could 
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generate  two  spatially  coherent  Stokes  beams.  Rodgers  demonstrated  this 
experimentally.  He  split  a  pump  beam  such  that  half  the  power  traveled  one  path  and  the 
other  half  traveled  a  different  path  into  a  single  fiber.  The  resulting  Stokes  beams  were 
shown  to  be  spatially  coherent  [9].  This  experiment  was  limited  by  the  laser  source  used. 
The  overall  power  of  the  pump  laser  was  relatively  low,  thus  beam  combining  could  not 
be  examined  in  a  highly  multimode  fiber.  Additionally,  the  laser  frequency  drifted, 
precluding  any  significant  examination  of  combining  two  different  lasers  at  the  same 
frequency. 

This  chapter  builds  upon  the  results  of  Rodgers  and  Chapter  4  with  an  investigation 
into  spatially  combining  two  laser  beams  using  SBS.  In  Section  5.2,  two  beams  at  the 
same  frequency  are  combined.  Initially,  the  two  beams  originate  from  different  lasers  to 
demonstrate  that  they  could  work  together  to  produce  the  Stokes  beam.  This  was 
followed  by  experiments  in  which  the  two  beams  originate  from  the  same  laser, 
simplifying  the  investigation  into  the  characteristics  of  the  spatial  combining  by 
eliminating  the  added  complexity  of  ensuring  the  laser  frequencies  do  not  drift  apart.  In 
Section  5.3  the  results  of  experiments  combining  two  beams  at  different  frequencies  are 
explained.  Similar  to  Rodger’s  research,  most  of  this  experiment  is  limited  to  fibers  with 
small  core  diameters,  but  the  results  still  clearly  show  that  the  two  beams  can  be  spatially 
combined. 

5,2.  Coherent  Combining 

Coherent  combining  is  the  process  of  combining  two  beams  at  nearly  the  same 
frequency.  It  requires  that  both  beams  work  together  to  generate  a  single  density 
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perturbation  in  the  fiber,  thus  their  frequeneies  must  fall  within  a  single  SBS  gain 
bandwidth,  -'100  MHz  [45].  Sinee  both  pump  beams  work  together,  the  overall 
effieieney  ean  be  higher  than  that  of  two  beams  eombined  ineoherently.  Additionally,  the 
Stokes  beam  is  mueh  more  speetrally  pure  and  will  therefore  have  a  mueh  longer 
eoherenee  length.  A  number  of  different  experiments  are  deseribed  in  this  seetion  that 
examine  the  properties  of  two  pump  beams  that  are  eoherently  eombined.  In  Seetion 
5.2.1,  the  eoupling  between  the  two  pump  beams  are  examined  to  demonstrate  their 
mutual  interaetion  that  results  in  a  temporally  eoherent  Stokes  beam.  The  polarization 
properties  of  the  Stokes  beams  are  examined  in  Seetion  5.2.2,  verifying  the  predietions  of 
the  model  deseribed  in  Seetion  4.3.2.  In  Seetion  5.2.3,  a  single  laser  is  split  into  two 
paths  and  reeombined.  The  quality  of  this  eombining  is  eharaeterized  both  qualitatively 
as  well  as  quantitatively.  The  seetion  also  deseribes  how  the  beam  eombination  is 
affeeted  by  the  position  of  the  pump  beams  on  the  front  fiber  faeet.  Finally,  eonelusions 
about  eoherent  eombining  of  two  laser  beams  are  diseussed  in  Seetion  5.2.4. 

5.2.1.  Two  laser  coherent  combining 

Experimental  Setup 

Coherent  beam  eombining  was  investigated  using  the  setup  shown  in  Figure  25.  The 
power  generated  by  two  Lightwave  Eleetronies  Nd;YAG  lasers  was  independently 
eontrolled  with  half-wave-plates  and  polarization  beam  splitters.  One  of  these  beams  was 
then  amplified  by  two  passes  through  a  Nd:YAG  rod.  The  two  beams  were  then  nearly 
spatially  eombined  using  a  polarization  beam  splitter  and  eoupled  into  a  50-mioron,  4.4 
km,  graded  index  fiber.  A  portion  of  the  seattered  beam  was  refieeted  off  a  plane  parallel 
plate  oriented  at  near  normal  ineidenee  to  the  beam.  A  seanning  Fabry-Perot 
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interferometer  monitored  the  frequencies  of  the  two  beams  to  ensure  they  were  close 
enough  for  coherent  combining.  The  power  of  the  Stokes  beam  was  also  monitored  as  a 
function  of  the  power  of  both  pump  beams. 


Beam 


Optical 

Isolator 


The  mutual  interaction  of  the  two  pump  beams  in  a  single  fiber  is  clearly  shown  using 
the  output  of  both  the  Fabry-Perot  and  the  power  meter.  Two  Fabry-Perot  scans  are 
shown  in  Figure  26.  The  top  trace  shows  the  Stokes  output  generated  by  the  two  lasers 
when  their  frequencies  are  separated  by  approximately  400  MHz.  This  is  more  than  four 
times  the  Brillouin  gain  bandwidth,  thus  two  individual  Stokes  beams  are  generated.  The 
scaling  along  the  vertical  axis  is  arbitrary.  However,  it  is  a  good  measure  of  the  relative 
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power  in  each  of  the  Stokes  beams  and  does  not  change  from  the  top  Fabry-Perot  scan  to 
the  bottom  scan.  The  heights  of  the  two  beams  seen  in  top  image  are  30  and  31.  In  the 
bottom  image,  the  frequencies  of  the  pump  beams  were  tuned  to  overlap,  thus  only  a 
single  Stokes  beam  is  observed.  Additionally,  the  Stokes  beam  relative  power,  100,  is 
found  to  be  more  than  the  incoherent  addition  of  the  two  individual  pump  beams,  61. 
This  synergy  is  a  result  of  the  two  beams  working  together  to  form  the  Stokes  output. 


b  26:  (Top)  Scanning  Fabry-Perot  interferometer  output  when  the  pump  beams  are 
I  separated  by  more  than  the  Brillouin  gain  bandwidth.  Relative  peak  output  of  the 
Stokes  beams  is  ~30.  (Bottom)  Interferometer  output  when  the  pump  beams  are 
tuned  to  the  same  frequency.  Peak  relative  output  is  100,  more  than  twice  that  of 

each  laser  alone. 


Coherent  combining  was  also  examined  using  the  measured  Stokes  power  as  a 
function  of  the  input  beam  powers.  This  data  is  shown  in  Figure  27,  illustrating  that  the 
two  pump  beams  collaborate  to  produce  a  single  Stokes  beam.  Figure  28  shows  the 
Stokes  power  generated  along  the  path  indicated  in  bold  on  the  surface  plot  of  Figure  27. 
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Initially,  the  power  of  the  second  laser  is  held  at  zero,  while  the  first  laser  is  increased  to 
its  maximum  power.  For  the  second  leg  of  the  path,  the  first  laser  is  held  at  maximum 
power  and  the  second  laser  is  increased  to  its  maximum  power.  This  path  is  insightful 
because  it  shows  that  while  the  first  laser  had  to  overcome  a  threshold  before  it  began  to 
contribute  to  the  Stokes  beam,  the  second  laser  does  not  and  contributes  to  the  Stokes 
beam  immediately.  This  is  a  direct  result  of  the  mutual  interaction  between  the  pump 
beams.  The  difference  in  slope  efficiency  for  the  two  lasers  will  be  explained  in  Section 
5.2.3. 
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Input  PowCT  (mW) 

Figure  28:  Stokes  power  generated  along  the  bold  path  shown  in  Figure  27. 

5.2.2.  Polarization  properties 
Experimental  Setup 

The  polarization  properties  of  the  eoherently  combined  Stokes  beam  were 
investigated  using  an  experimental  setup  very  similar  to  that  shown  in  the  previous 
section.  The  only  significant  variation  is  in  the  tools  used  to  analyze  the  reflected  beam. 
The  power  meter,  50/50  beamsplitter,  and  Fabry-Perot  were  replaced  by  a  quarter-wave- 
plate,  half-wave-plate,  analyzer,  and  a  power  meter  arranged  in  the  configuration  shown 
by  Figure  29.  From  this  configuration,  a  complete  description  of  the  beam  polarization 
was  ascertained  through  a  series  of  measurements  at  various  quarter-wave-  and  half¬ 
wave-plate  angles.  These  angles  correspond  to  measuring  specific  polarization 
components  of  the  Stokes  beam  and  are  tabulated  below  (Table  3).  From  these 
measurements,  the  Stokes  parameters  were  calculated  using  the  relationships  in  Table  4 
[46].  The  measurements  were  repeated  for  each  power  level  shown  along  the  bold  path 
of  Figure  27. 


75 


Power  Meter 


Stokes  beam 
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29:  Polarization  analysis  setup. 


Table  3:  Measured  polarization  components. 


X/2  angle  (°) 

X/4  Angle  (°) 

Polarization 

Component 

Label 

0 

0 

Florizontal 

H 

45 

0 

Vertical 

V 

0 

45 

Right  Hand  Circular 

RHC 

45 

45 

Left  Hand  Circular 

LHC 

22.5 

45 

Linear  @  45° 

L" 

-22.5 

45 

Linear  @  -45° 

U 

Table  4:  Conversion  from  the  measured  values  to  Stokes  parameters. 


Stokes 

Parameter 

So 

Si 

S2 

S3 

Equation 

+  V  +  RHC  +  LHC  +  n  +  r) 

H-V 

H-r 

RHC -LHC 

Results 

The  experimental  results  were  analyzed  to  determine  if  the  coherently  combined 
beams  follow  the  modeling  described  in  Section  4.3.2  for  a  single  pump.  In  this  model, 
the  X  -  and  y  -components  of  the  Stokes  beam  were  treated  independently  except  to 
ensure  that  they  reflect  off  the  same  density  perturbation  in  the  fiber.  The  experimental 
setup  is  well  suited  for  this  analysis  because  the  two  pump  beams  are  orthogonally 
polarized.  Thus  the  polarization  model  predicts  that  the  horizontally  polarized  pump 
would  only  contribute  to  the  horizontally  polarized  Stokes  beam  and  vice  versa.  The 
measurements  are  consistent  with  this  theory  and  are  shown  in  Figure  30.  While  the 
power  of  the  first  laser  is  increased,  the  Stokes  beam  is  predominantly  horizontally 
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polarized.  Likewise  when  the  power  of  the  seeond  laser  is  inereased,  the  Stokes  beam 
gains  a  strong  vertical  component. 


L2  =  0  L2  =  0  L2  =  0  L2  =  0  L2  =  60  L2  =  160  L2  =  260 

Input  Power  (mW) 

'  f  30:  (solid)  Total  Stokes  power  as  a  function  of  the  pump  power  along  the  path 
Idicated  in  bold  on  Figure  27.  (dotted)  Horizontally  polarized  Stokes  power  as  a 
function  of  pump  power,  (dashed)  Vertically  polarized  Stokes  power  as  a  function 

of  pump  power. 


While  the  magnitudes  of  the  vertical  and  horizontal  components  of  the  beams  are 
easily  observed,  they  do  not  form  a  complete  description  of  the  polarization  of  the  Stokes 
beam.  A  complete  description  is  only  determined  when  the  phase  difference  between  the 
two  components  is  known.  However,  since  the  two  pump  beams  are  created  by  separate 
oscillators,  each  with  a  linewidth  of  ~5  kHz,  the  polarization  of  the  pump  beam  is 
expected  to  randomly  shift  in  time.  Even  if  the  center  frequencies  of  the  two  beams  are 
perfectly  aligned,  there  will  be  a  change  in  polarization  on  a  time  scale  of 

YskHz  ~  •  The  coherence  time  will  be  even  shorter  if  the  pump  beam  frequencies 

do  not  overlap  perfectly.  The  degree  of  polarization,  V  ,  takes  into  account  the  random 
shifts  in  polarization.  It  is  defined  in  terms  of  the  Stokes  vector  as  [46] 
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(5.1) 


V  = 


If  the  polarization  of  the  beam  does  not  change  in  time,  this  equation  is  identically  equal 
to  one.  On  the  other  hand,  if  there  are  random  fluctuations  in  the  polarization,  the  degree 
of  polarization  drops. 


Using  the  definitions  for  the  elements  of  the  Stokes  vector  in  Table  4  and  assuming 
that  the  phase  difference  between  the  horizontal  component  (laser  1)  and  vertical 
component  (laser  2)  of  the  total  pump  beam  is  random,  the  degree  of  polarization  for  the 
pump  beam  is 


V  = 


H-V 
H  +  V  ■ 


(5.2) 


This  theoretical  value  is  represented  by  the  dotted  line  in  Figure  31. 

As  described  in  Section  4.3.2,  the  phase  difference  between  the  x  -  and  y  - 
components  of  the  pump  beam  is  conserved  after  Brillouin  reflection.  Thus  the  Stokes 
degree  of  polarization  should  mimic  that  of  the  pump.  The  solid  line  of  Figure  31  shows 
the  degree  of  polarization  of  the  Stokes  beam  as  calculated  using  the  measured  Stokes 
vectors.  As  soon  as  the  reflected  power  in  the  first  laser  is  detectable,  the  degree  of 
polarization  is  relatively  high.  However,  when  the  second  laser  is  turned  on,  the  degree 
of  polarization  drops,  consistent  with  that  of  the  pump  beam. 

All  of  the  experimental  evidence  agrees  with  the  trends  predicted  by  the  polarization 
model  described  in  the  previous  chapter.  However,  the  experiments  described  in  this 
section  cannot  definitively  prove  that  the  polarization  of  the  Stokes  beam  was  identical  to 
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3 1 :  The  degree  of  polarization  (solid  line)  falls  off  rapidly  as  the  second  laser 
begins  to  pump  the  fiber,  contributing  to  the  overall  Stokes  output.  This  is 
relatively  consistent  with  the  predicted  degree  of  polarization  of  an  idealized  pump 
beam  (dotted  line).  The  total  reflected  Stokes  power  rises  linearly  above  threshold 
regardless  of  the  polarization  states  of  the  pump  beams  (dashed  line). 


that  of  the  pump  beam.  In  fact,  when  only  a  single  laser  pumped  the  fiber,  the  pump 
beam  was  nearly  100%  polarized,  but  the  Stokes  beam  was  found  to  be  approximately 
80%  polarized.  It  is  believed  that  this  discrepancy  was  caused  by  power  fluctuations  in 
the  Stokes  beam.  The  measurement  of  the  Stokes  parameters  required  six  independent 
power  measurements.  Any  change  in  the  Stokes  power  from  one  of  these  measurements 
to  another  would  result  in  a  reduced  degree  of  polarization.  This  limitation  in  the 
experiment  could  be  remedied  by  a  better  sample  of  data  points.  Unfortunately,  that 
would  require  laser  frequency  stability  over  a  longer  time  frame. 


5.2.3.  Spatial  characteristics 
Experimental  Setup 

Two  lasers  were  critical  to  determine  the  characteristics  of  the  Stokes  beams 
described  in  Sections  5.2.1  and  5.2.2.  However,  ensuring  coherent  combination  was 
difficult  because  the  laser  frequencies  slowly  drifted  out  of  resonance  with  each  other. 
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To  eliminate  this  problem,  the  remaining  experiments  were  eompleted  with  a  single  laser 
split  into  two  independent  paths.  This  allowed  for  a  elear  demonstration  of  beam 
eombining  and  an  aeeurate  measurement  of  the  eombining  effieieney  as  a  funetion  of  the 
pump  spatial  position. 

A  sehematic  of  the  experiment  used  to  investigate  the  spatial  charaeteristies  of 
eoherent  eombining  is  shown  in  Figure  32.  The  beam  from  a  700  mW  laser  with  an 
instantaneous  speetral  width  ~5  kHz  was  amplified  by  a  Nd:YAG  amplifier  and  split 
using  a  polarization  beam  splitter  into  two  paths.  The  beams  were  then  nearly  spatially 
eombined  using  a  seeond  polarization  beam  splitter,  and  eoupled  into  a  50-mioron,  4.4 
km,  graded-index  fiber.  To  demonstrate  beam  eombining,  it  was  eritieal  that  the  pump 
beams  were  not  spatially  eombined  with  the  polarization  beam  splitter.  Thus  the  beams 
were  spatially  separated  by  as  mueh  as  possible  while  still  ensuring  that  they  were 
eoupled  into  the  fiber.  After  refieetion  by  the  fiber,  a  portion  of  the  beam  was  observed 
with  a  CCD  eamera  while  a  Fabry-Perot  and  a  power  meter  were  used  to  determine  the 
speetrum  and  relative  power  of  the  Stokes  beams. 

Several  measurements  of  the  refleeted  beams  were  taken  with  the  CCD  eamera  to 
determine  the  degree  of  spatial  eoherenee  in  the  Stokes  beams.  First,  images  of  the 
Stokes  beams  generated  by  each  path  were  collected.  These  depict  the  relative  position 
of  the  Stokes  beams  generated  by  the  two  pump  beams.  Then,  without  adjusting  the  front 
fiber  facet,  the  fiber  was  cut  to  be  ~2  m  long.  The  shortened  interaction  length  increased 
SBS  threshold  well  above  the  power  delivered  by  the  two  beam  paths.  Two  more  CCD 
images,  one  for  each  pump  beam,  were  taken  of  the  Fresnel  reflection  off  the  front  facet 
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Figure  32:  Two  beam  coherent  combining  schematic. 


of  the  fiber.  These  images  were  used  to  subtract  the  contribution  of  the  Fresnel  reflection 
from  the  images  of  the  Stokes  beam,  as  well  as  allowed  for  a  measurement  of  the  degree 
of  spatial  coherence  of  the  two  pump  beams.  Finally,  to  ensure  that  the  beams  were 
spatially  combined  through  SBS  rather  than  some  other  process,  such  as  spatial  filtering, 
the  CCD  camera  was  moved  to  the  output  end  of  the  2  m  fiber  and  the  two  transmitted 
pump  beams  were  observed. 

The  beam  combining  efficiency  was  also  investigated.  This  was  completed  in  a 
similar  manner  to  the  efficiency  measurements  described  in  Chapter  4.  While  path  one 
was  held  constant,  the  second  beam  was  translated  along  the  front  fiber  facet  by  rotating 
the  mirror  in  this  path.  At  each  point,  the  Stokes  power  was  measured  at  the  CCD 
camera  and  an  image  of  the  Stokes  beam  was  collected.  The  Stokes  beam  was  much 
stronger  than  the  reflected  pump  beam,  thus  the  pump  beam  was  not  noticeable  in  the 
images.  In  order  to  determine  the  location  of  the  pump  beam  on  the  fiber,  the  long  fiber 
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was  again  cut  to  meters  and  the  position  of  the  pump  beam  observed.  This  also 
allowed  for  a  measurement  of  the  transmitted  power  and  hence  a  determination  of  the 
fiber  core  size. 

Results 

The  results  of  the  coherent  beam  combining  experiment  are  shown  in  Figure  33- 
Figure  35.  In  each  figure,  part  (a)  is  an  image  of  the  relative  beam  intensity  due  to  the 
beam  injeeted  from  the  first  path.  In  this  image,  darker  shades  indicate  higher  intensities. 
The  total  area  under  the  intensity  distribution  is  normalized  to  one,  whieh  eliminates 
fluctuations  from  one  image  to  the  next  due  to  differenees  in  the  power.  The  second 
image  in  each  figure,  part  (b),  shows  the  relative  beam  intensity  due  to  the  second  path. 

In  this  image,  lighter  shading  indicates  higher  intensities  and  again  the  area  is  normalized 
to  account  for  power  fluctuations.  Finally,  part  (c)  is  the  combined  image  of  parts  (a)  and 
(b).  It  represents  the  pixel-by-pixel  subtraction  of  the  intensity  of  one  path  from  the 
other.  Keeping  the  convention  of  parts  (a)  and  (b),  dark  regions  indicate  a  higher 
concentration  of  power  due  to  the  first  path  and  brighter  regions  indicate  a  higher 
intensity  due  to  the  second  path.  Gray  areas  result  when  there  is  equal  power  from  each 
path. 

Figure  33  indicates  the  relative  positions  of  the  pump  beams.  Interpreting  part  (c)  as 
described  above,  clearly  shows  that  the  overlap  of  the  two  beams  is  insignificant.  It  can 
be  inferred  from  this  picture  that  the  pump  beams  are  not  coaxial  before  entering  the 
fiber.  The  relative  positions  of  the  transmitted  beams  are  shown  in  Figure  34.  The 
structure  in  this  figure  demonstrates  that  even  though  there  is  some  overlap,  the  beams 
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(a)  (b)  (c) 

Figure  33:  Images  of  the  pump  beams  from  each  path,  a)  The  pump  beam  from  path  1 
(high  intensities  are  indicated  by  dark  regions),  b)  The  pump  beam  from  path  2 
(high  intensities  are  indicated  by  bright  regions,  c)  Image  of  both  pump  beams. 


(a)  (b)  (c) 

Figure  34:  Images  of  the  transmitted  beams  from  each  path,  a)  The  transmitted  beam 
from  path  I  (high  intensities  are  indicated  by  dark  regions),  b)  The  transmitted 
beam  from  path  2  (high  intensities  are  indicated  by  bright  regions,  c)  Image  of 

both  transmitted  beams. 


(a)  (b)  (c) 

:  35:  Images  of  the  Stokes  beams  from  each  path,  a)  The  Stokes  beam  from  path  I 
(high  intensities  are  indicated  by  dark  regions),  b)  The  Stokes  beam  from  path  2 
(high  intensities  are  indicated  by  bright  regions,  c)  Image  of  both  Stokes  beams. 
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are  not  spatially  combined  upon  entering  the  fiber.  Finally,  Figure  35  shows  the  Stokes 
beams  generated  by  each  path.  Part  (c)  of  this  figure  shows  very  little  structure, 
representing  the  fact  that  the  two  paths  generate  nearly  identical  Stokes  beam  shapes. 

A  more  formal  analysis  of  the  beam  combining  results  was  developed  to  quantify 
how  well  the  beams  were  combined.  This  analysis  uses  the  mean  square  error  between 
the  normalized  intensity  from  path  1  and  the  normalized  intensity  from  path  2.  This  is 
defined  by 

(5.3) 

Ij  and  Ig  represent  the  intensity  of  each  path  over  the  range  of  pixels  in  the  image,  and 
N  is  the  total  number  of  pixels.  The  results  of  these  calculations  are  shown  in  Table  5. 
The  Stokes  beams  are  found  to  have  a  misalignment  that  is  about  15  times  lower  than  the 
transmitted  beams  and  over  500  times  lower  than  the  pump  beams. 


Table  5:  Misalignment  calculation 


Pump 

SBS 

Transmitted 

3270 

5.98 

91.8 

Relative  Misalignment 

547 

1 

15 

A  similar  analysis  was  completed  using  aberrated  pump  beam  for  each  input.  In  this 
experiment,  both  pump  beams  propagated  through  a  piece  of  Saran  wrap  before  coupling 
into  the  fiber.  The  results  support  the  same  conclusions  as  above,  but  were  not  as 
dramatic  for  several  reasons.  First,  scattering  off  the  Saran  wrap  reduced  the  power 
coupled  into  the  fiber,  thus  making  the  ratio  of  Stokes  power  to  pump  Fresnel  reflection 
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lower.  Secondly,  as  will  be  shown  in  the  following  pages,  the  lowered  power  required 
the  pump  beams  to  be  closer  together  to  excite  the  Stokes  beams.  Finally,  the  aberrations 
on  the  beams  increased  their  physical  size,  causing  greater  overlap  of  the  input  beams. 

When  the  pump  beams  are  not  coupled  into  the  center  of  the  fiber  core,  the  combining 
efficiency  suffers.  Figure  36  shows  the  relative  Stokes  power  as  the  path-two  pump 
beam  was  moved  from  one  side  of  the  fiber  to  the  other.  The  top  of  the  figure  shows  the 
position  of  the  pump  relative  to  the  Stokes  beam  for  three  of  the  data  points.  As  the 
figure  shows,  the  efficiency  has  a  form  similar  to  that  observed  for  beam  cleanup.  This 
implies  that  although  the  contribution  from  power  near  the  edge  of  the  fiber  will  be  low, 
all  power  coupled  into  the  fiber  will  contribute  to  the  Stokes  beam. 


36:  (Top)  Stokes  beam  (black)  as  a  function  of  the  position  of  the  path-two  pump 
(white).  (Bottom)  Relative  SBS  efficiency  as  a  function  of  the  path-two  beam 
position.  The  power  coupled  into  the  fiber  is  shown  by  the  solid  line. 


5.2.4.  Conclusions 

As  the  data  has  shown,  coherent  beam  combining  is  possible  using  stimulated 
Brillouin  scattering.  In  Section  5.2.1,  the  two  pump  beams  were  shown  to  mutually 
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excite  the  fiber  when  their  frequencies  overlapped.  This  resulted  in  a  Stokes  beam  at  a 
single  frequency  and  a  higher  Brillouin  reflection  coefficient  than  when  the  pump  beams 
operated  independently. 

The  polarization  characteristics  of  the  Stokes  beam  were  examined  in  Section  5.2.2. 
The  results  of  this  measurement  indicated  that  the  pump  beams  retain  their  original 
polarization  after  Brillouin  reflection.  This  was  demonstrated  through  a  measurement  of 
the  components  of  the  Stokes  beam  that  were  co-  and  cross-polarized  with  respect  to  each 
pump  beam.  Additionally,  the  total  pump  beam  degree  of  polarization  was  found  to 
match  that  of  the  Stokes  output  relatively  well.  This  was  expected  based  upon  the 
polarization  model  described  in  the  previous  chapter. 

Finally,  in  Section  5.2.3,  it  was  shown  that  when  two  beams  were  coupled  into  the 
fiber  at  different  positions  and  angles,  the  resulting  Stokes  beams  have  the  same 
frequency  and  are  spatially  combined.  The  efficiency  of  the  process  was  shown  to  mimic 
that  of  a  single  aberrated  beam  with  two  lobes.  When  coupled  into  the  center  of  the  fiber, 
the  efficiency  was  greater  than  when  the  pump  was  coupled  into  the  edge  of  the  fiber. 

5.3.  Incoherent  Combining 

As  described  in  the  previous  chapter,  in  order  for  two  beams  to  coherently  combine, 
the  frequencies  of  both  beams  must  fall  within  a  single  Brillouin  gain  bandwidth.  Even 
in  a  laboratory  setup  under  ideal  conditions,  this  can  be  difficult  to  obtain  for  two 
independent  lasers.  Scaling  this  solution  to  even  more  beams  could  become  prohibitively 
complex.  On  the  other  hand,  if  each  pump  laser  operates  at  a  different  frequency,  they 
will  excite  independent  Stokes  beams  in  the  fiber  [45,  47-49].  Each  of  these  Stokes 


86 


beams  should  be  a  eleaned  version  of  the  pump  and  is  expeeted  to  propagate  in  the  LPoi 
mode.  Thus  the  spatially  ineoherent  pump  beams  would  produee  spatially  eoherent 
Stokes  beams.  Beeause  the  Stokes  beams  are  eaeh  exeited  at  different  frequeneies,  they 
are  temporally  incoherent  with  respect  to  one  another.  Thus  this  is  referred  to  as 
incoherent  beam  combining. 

Incoherent  beam  combining  can  be  used  to  scale  beam  power  to  higher  levels  than 
those  achievable  through  coherent  beam  combining.  If  one  were  to  coherently  combine 
many  lasers  into  a  single  beam,  the  resulting  Stokes  wave  would  eventually  reach  second 
order  SBS  threshold,  resulting  in  decreasing  efficiency  with  increasing  input  power. 
While  this  is  not  expected  to  be  a  severe  limitation,  at  some  point,  its  effect  would  have 
to  be  considered.  Incoherent  combining  would  relax  this  limitation  significantly.  For  the 
incoherent  process,  each  Stokes  beam  is  produced  independent  of  the  others,  thus  the 
second  order  threshold  is  dependent  on  the  power  in  each  Stokes  beam,  not  the  total 
Stokes  power  and  second  order  Stokes  threshold  will  be  greater  than  that  for  a  single 
coherent  beam  by  a  factor  of  the  number  of  lasers  being  combined  [45]. 

Moreover,  this  technique  is  completely  general  in  that  lasers  of  any  wavelength  can  be 
combined  as  long  as  their  wavelengths  he  in  the  transmission  window  of  the  fiber 
materials.  This  kind  of  beam  combining  technique  would  be  useful  for  any  applications 
that  require  high  brute-force  laser  power  in  a  single  beam  such  as  in  material  processing. 
It  would  also  be  appropriate  for  applications  that  require  multiple  discrete  wavelengths  in 
a  single  beam,  such  as  in  multi-spectral  electro-optic  countermeasures  or  hyper-spectral 
imaging.  The  spatial  coherence  that  this  technique  provides  makes  it  possible  for  the 
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combined  beam  to  be  direeted  into  a  diffraetion-limited  spot,  be  it  a  tight  focus  or  long¬ 
distance  targeting. 

The  investigation  of  ineoherent  beam  eombining  was  eompleted  in  two  parts.  In 
Seetion  5.3.1,  a  set  of  experiments  are  deseribed  that  were  designed  to  investigate  the 
mutual  interaction  of  the  two  lasers  in  a  single  fiber.  The  goal  of  the  experiments  was  to 
demonstrate  that  sinee  the  frequeney  of  the  two  beams  are  separated  by  greater  than  the 
Brillouin  gain  bandwidth,  they  exeite  independent  Stokes  waves  in  the  fiber.  In  Seetion 
5.3.2,  the  spatial  overlap  of  the  Stokes  waves  generated  by  the  two  laser  beams  is 
eharaeterized.  In  this  seetion  it  is  shown  that  although  the  pump  beams  do  not  overlap, 
the  two  Stokes  beams  have  exeellent  spatial  overlap.  The  spatial  eoherenee  of  the  Stokes 
beams  is  examined  in  Section  5.3.3.  In  this  experiment,  the  Stokes  beam  is  split  and 
travels  two  paths  to  a  CCD  eamera.  The  interferenee  of  the  two  beams  is  then  analyzed 
as  the  paths  change.  Finally,  conclusions  concerning  ineoherent  beam  combining  are 
presented  in  Seetion  5.3.4. 

5.3.1.  Independent  excitation 
Experiment 

The  setup  used  to  demonstrate  ineoherent  beam  combining  is  shown  in  Figure  37. 

Two  SDL  single  mode  (longitudinal)  lasers  were  collimated  using  20x  mieroscope 
objeetives.  Laser  A  emitted  light  at  853  nm  while  laser  B  operated  at  832  nm.  The 
eollimated  beams  passed  through  40-dB  Faraday  optieal  isolators  and  were  then  spatially 
superimposed  with  two  half-wave  plates  and  a  polarizing  beam  splitter.  At  the  beam 
splitter,  lasers  A  and  B  are  polarized  horizontally  and  vertieally,  respeetively.  After  the 
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beam  splitter,  the  light  passed  through  a  quarter-wave  plate  oriented  to  convert  the 
crossed  linear  polarizations  of  the  two  beams  into  orthogonal  circular  polarizations. 
Finally,  using  a  lOx  microscope  objective,  the  beams  were  focused  into  a  4.4  km,  step- 
index,  Corning  SMF-28  fiber  with  an  8.2-micron  core  and  a  numerical  aperture  of  0.12. 
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^37:  Two  beam  incoherent  combination  experimental  setup  using  a  multimode 

fiber. 


The  light  that  was  reflected  off  the  surface  of  the  fiber  by  Fresnel  reflection  as  well  as 
the  light  generated  through  SBS  was  collimated  as  it  passed  back  through  the  microscope 
objective  and  the  quarter  wave  plate.  A  portion  of  this  reflected  off  both  the  front  and 
back  surfaces  of  an  uncoated  fiat  window.  The  fraction  of  the  beam  that  was  reflected  off 
of  the  rear  surface  of  the  flat  passed  through  a  half-wave  plate,  a  Gian- Thompson 
polarizer  and  neutral  density  filters  before  being  captured  by  a  CCD  camera.  The  portion 
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of  the  beam  that  was  refleeted  from  the  front  surfaee  of  the  flat  was  picked  off  using  a 
mirror  and  directed  through  a  Gian- Thompson  polarizer  into  either  a  Coherent  power 
meter  or  a  Fabry-Perot  interferometer  for  spectral  characterization.  After  monitoring  the 
position  of  each  beam,  the  fiber  was  cut  -'2  m  beyond  the  front  fiber  facet  and  more 
intensity  profiles  were  collected.  These  images  were  used  to  extract  the  intensity  of  the 
pure  Fresnel  reflection  from  that  of  the  Stokes  beam  since  the  shortened  fiber  increased 
SBS  threshold  well  beyond  the  power  produced  by  the  laser. 

Results 

Figure  38(a)  shows  a  2-D  surface  plot  displaying  the  power  of  the  horizontally 
polarized  Stokes  beam  as  it  enters  the  Fabry-Perot  over  all  possible  input  powers  for  laser 
A  and  laser  B.  The  separation  of  the  Fresnel  reflection  from  the  Stokes  beam  was 
facilitated  by  the  fact  that  the  polarization  of  the  Fresnel  reflection  is  orthogonal  to  that  of 
the  Stokes,  as  discussed  in  the  previous  chapter.  Obviously,  as  the  current  on  laser  A 
increases  (thus  laser  power),  the  horizontally  polarized  Stokes  beam  (co-polarized  with 
the  laser  A  input)  increases  as  well,  once  threshold  is  passed.  There  is  no  significant 
change  to  the  Stokes  power  due  to  laser  A  as  the  power  coupled  into  the  fiber  from  laser 
B  is  increased.  Figure  38(b)  demonstrates  that  the  converse  is  also  true.  As  the  current 
driving  laser  B  is  increased,  the  power  in  the  vertically  polarized  beam  (co-polarized  with 
the  laser  B  input)  increases,  independent  of  the  laser  A  power  coupled  into  the  fiber.  The 
combined  output  of  both  polarizations  is  shown  in  Figure  38(c).  ff  one  were  to  trace  out 
the  path  along  this  surface  shown  by  the  bold  line,  distinct  threshold  powers  are  seen  for 
each  laser.  This  path  is  plotted  in  Figure  39,  demonstrating  that  each  laser  must  exceed 
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threshold  in  order  to  excite  its  own  Stokes  beam.  This  behavior  contrasts  with  that  seen 


during  incoherent  combining  when  a  combined  threshold  was  observed  for  both  lasers 
(Figure  28). 


;  38:  Polarimetric  Stokes  Power  vs.  Pump  Power  in  mW.  (a)  SBS  Power  refleeted 
into  a  horizontally  polarized  beam,  (b)  SBS  Power  reflected  into  a  vertically 
polarized  beam,  (c)  Total  SBS  Power. 


Pump  Power  (mW) 

39:  The  indicated  path  of  Figure  38(c)  shows  that  the  two  input  beams  must  each 
overcome  a  threshold  before  contributing  to  the  Stokes  beam. 


Several  slices  of  Figure  38(c)  were  taken  along  both  axes,  to  determine  the  possible 
dependence  of  the  SBS  threshold  for  one  laser  on  the  power  of  the  other.  The  slices  were 
normalized  such  that  the  left  hand  side  of  each  slice,  corresponding  to  no  cross  laser 
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input,  was  set  to  zero.  These  sliees  are  all  plotted  using  a  common  axis  in  Figure  40.  The 
similarity  of  the  curves  indicates  that  the  slope  efficiency  and  threshold  for  each  laser  is 
independent  of  the  power  in  the  other  laser.  Investigating  the  slices  along  the  laser  B  axis 
further  (Figure  40(a))  shows  that  the  average  slope  efficiency  of  these  curves  is  68  ±  8%, 
the  error  margin  being  one  standard  deviation.  The  standard  deviation  associated  with 
the  calculation  of  each  slope  efficiency  is  on  average  7%.  The  average  threshold  is 
17  ±  2mW .  The  standard  deviation  associated  with  the  calculation  of  the  threshold  is  on 
average  5  mW.  This  data  shows  that  the  threshold  and  slope  efficiency  are  nearly 
identical  within  the  experimental  uncertainty  over  the  full  range  of  laser  A  input  powers, 
indicating  the  complete  independence  between  the  Stokes  beam  generated  by  laser  B  and 
the  power  of  laser  A.  The  same  analysis  was  completed  for  the  Stokes  power  generated 
by  laser  A,  demonstrating  no  dependence  of  its  threshold  and  slope  efficiency  on  the 
power  of  laser  B. 

The  frequency  of  all  the  beams  was  continuously  monitored  throughout  these 
measurements  with  a  Fabry-Perot  spectrum  analyzer.  A  sample  spectrum  containing  the 
pumps  and  their  Stokes  components  is  shown  in  Figure  41 .  It  was  observed  that  the 
Stokes  shift  remained  constant  as  the  power  of  both  lasers  was  varied. 

The  observations  described  above  collectively  lead  to  the  conclusion  that  when  two 
laser  beams  of  arbitrary  frequencies  simultaneously  excite  SBS  in  a  common  fiber,  they 
individually  excite  their  corresponding  Stokes  beams  without  interfering  with  each  other. 
The  simultaneous  presence  and  generation  of  the  Stokes  beams  by  two  incoherent  pump 
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Relative  Power 


beams  in  a  single  fiber  have  no  effect  on  the  wavelength,  threshold,  or  efficiency  of  the 
Stokes  beam  generated  by  each  pump. 


(a)  (b) 

40:  Slices  through  Figure  38(c)  at  different  values  of  laser  B  (A)  and  laser  A  (B). 
ach  slice  is  normalized  at  the  left  side.  The  similarity  of  each  curve  indicates  that 
the  lasers  produce  independent  Stokes  beams. 


Frequency 

41 :  Dual  pump  and  SBS  spectrum  generated  from  a  Fabry-Perot  interferometer 
with  a  free  spectral  range  of  8  GFIz. 
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5.3.2.  Spatial  Overlap 
Experimental  Setup 

A  set  of  experiments  was  performed  to  investigate  the  spatial  coherenee  properties  of 
the  Stokes  beam  generated  through  the  incoherent  beam  combining  technique.  The  goal 
of  the  experiments  was  to  demonstrate  the  near-perfect  spatial  coherence  of  the  Stokes 
beams  generated  in  the  SBS  process  by  two  arbitrary  lasers  in  a  multi-mode  fiber. 

The  experimental  setup  was  identical  to  Figure  37  except  the  quarter-wave  plate  was 
removed.  Under  such  conditions,  all  light  from  a  single  pump  (the  Stokes  beam  and  the 
Fresnel  reflection)  had  the  same  polarization  and  could  be  kept  separate  from  the  second 
pump.  The  Corning  SMF-28  fiber  used  had  a  fiber  parameter  of  4.2,  indicating  that  four 
modes  (FPoi,  FPn,  FP02,  and  FP21)  could  be  supported  by  the  fiber  [50]. 

The  spatial  coherence  properties  were  analyzed  in  the  same  manner  as  those  for 
incoherent  combining.  First,  each  laser  was  turned  on  well  above  their  respective  SBS 
thresholds.  The  half-wave  plate  located  in  front  of  the  beam  profiler  was  tuned  for 
maximum  transmission  of  laser  A  and  the  resulting  Stokes  beam.  At  that  point,  eight 
intensity  profiles  were  captured  and  averaged.  The  averaging  was  necessary  due  to  the 
temporal  fluctuations  that  accompanied  the  Stokes  beam.  This  averaged  profile  consisted 
of  energy  in  the  Stokes  beam  created  by  laser  A  as  well  as  energy  that  reflected  off  the 
front  fiber  face.  The  half-wave  plate  was  then  rotated  45°  and  similar  images  were 
collected  showing  the  Stokes  beam  intensity  profile  due  to  laser  B  and  its  associated 
pump  beam. 
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Upon  completion  of  the  first  set  of  images,  the  fiber  length  was  eut  to  two  meters  from 
the  entranee  faee.  This  was  done  without  altering  the  position  of  the  front  faee  of  the 
fiber  or  adjusting  the  eamera  loeation.  Two  more  sets  of  intensity  profiles,  identieal  to 
the  first  were  then  eolleeted.  However,  sinee  the  fiber  was  only  two  meters  long,  the  SBS 
threshold  was  raised  well  above  the  operating  power  of  the  two  lasers.  This  ensured  that 
the  profiles  ineluded  light  due  only  to  Fresnel  refleetion  off  the  air-fiber  interfaee.  These 
profiles  show  the  relative  position  of  the  two  input  beams  and  were  subtraeted  from  the 
first  set  to  obtain  an  intensity  profile  due  purely  to  the  Stokes  beam. 

The  last  set  of  profiles  was  taken  at  the  exit  faee  of  the  fiber.  Eaeh  laser  was 
independently  powered  when  the  images  were  eolleeted.  The  two  profdes  elearly  show 
that  the  input  beams  were  not  spatially  filtered  into  the  fundamental  fiber  mode. 

Results 

Both  the  qualitative  measure  of  the  beam  overlap  and  the  numerieal  figure  of  merit 
used  to  analyze  eoherent  beam  eombining  were  applied  to  the  experimental  data.  The 
qualitative  measure  was  a  series  of  grayseale  plots  that  show  the  relative  position  and 
intensity  of  eaeh  laser.  The  brightness  of  the  plots  was  ealeulated  using  the  following 
algorithm.  The  intensity  image  from  eaeh  laser  was  separately  seanned  through  a  1 0 
pixel  X  10  pixel  mean  filter.  The  size  of  the  filter  was  ehosen  to  be  large  enough  to 
remove  the  random  noise  from  the  intensity  data  but  small  enough  to  retain  real  intensity 
variations  present  in  the  images.  The  power  of  the  beam  in  eaeh  intensity  image  was  then 
separately  normalized  to  one.  The  results  of  this  analysis  are  shown  in  Figure  42-Figure 
44.  In  eaeh  of  these  figures,  (a)  shows  the  intensity  profile  of  laser  A,  where  a  brighter 
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(a)  (b)  (c) 

Figure  42:  Reflected  beam  overlap,  (a)  Image  of  Laser  A.  Lighter  shades  indicate 
higher  intensity,  (b)  Image  of  Laser  B.  Darker  shades  indicate  higher  intensity, 
(c)  Combined  laser  beams  showing  significant  mismatch  between  the  reflected 
positions  of  laser  A  and  laser  B. 


(a)  (b)  (c) 

Figure  43:  Transmitted  beam  overlap,  (a)  Image  of  Laser  A.  Lighter  shades  indicate 
higher  intensity,  (b)  Image  of  Laser  B.  Darker  shades  indicate  higher  intensity, 
(c)  Combined  laser  beams  showing  significant  mismatch  between  the  transmitted 
positions  of  laser  A  and  laser  B. 


(a)  (b) 

44:  SBS  beam  overlap,  (a)  Image  of  Laser  A. 
intensity,  (b)  Image  of  Laser  B.  Darker  shades 
Combined  laser  beams  showing  the  two  SBS 


(c) 

Lighter  shades  indicate  higher 
indicate  higher  intensity,  (c) 
beams  are  nearly  coaxial. 
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output  indicates  a  higher  intensity.  Part  (b)  of  these  figures  shows  the  intensity  profile  of 
laser  B.  These  images  are  shaded  sueh  that  darker  regions  indieate  higher  intensities. 
Finally,  part  (e)  is  the  addition  of  both  the  images  from  (a)  and  (b).  This  image 
eorresponds  to  a  pixel-by -pixel  image  subtraetion  of  beam  B  from  beam  A. 

Figure  42  shows  the  relative  positions  of  the  Fresnel  refleetion  off  the  front  fiber  faeet. 
This  image  was  eolleeted  from  the  fiber  after  it  was  cut  to  approximately  two  meters. 

The  loeation  of  the  refleeted  beams  on  the  eamera  indieates  of  the  relative  position  of  the 
input  beams  on  the  front  faee  of  the  fiber.  As  ean  be  seen  by  the  obvious  strueture  in 
Figure  42(e),  there  is  a  signifieant  mismateh  between  the  loeations  of  the  output  from 
laser  A  and  laser  B.  This  elearly  indieates  that  the  two  lasers  were  not  eoaxially  aligned 
before  entry  into  the  fiber.  The  transmitted  intensity  profiles  are  shown  in  Figure  43.  ft 
is  apparent  by  the  relative  displaeement  of  the  two  beams  that  spatial  filtering  due  to  the 
fiber  was  not  responsible  for  the  alignment  of  the  laser  beams.  Finally,  Figure  44  shows 
the  Stokes  beams  generated  by  the  two  lasers.  Both  beams  appear  to  be  emitted  from  a 
Gaussian-like  LPqi  fiber  mode  and  as  ean  be  seen  from  Figure  44(e),  the  two  beams  have 
nearly  identieal  spatial  intensity  distributions. 

While  visually  gratifying,  the  above  analysis  does  not  offer  a  quantitative  measure  of 
the  two-beam  alignment.  In  order  to  define  a  measurable  quantity  that  would  indieate 
how  well  the  two  laser  beams  were  aligned  under  the  different  seenarios,  the  images  were 
again  separately  normalized  but  then  the  mean  square  deviation  between  the  two  images 
was  ealeulated.  This  deviation  is  defined  as  follows. 
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(4) 


where  I^  g  (x,.)  indicates  the  intensity  of  the  pixel  at  x,.  associated  with  laser  A,B  and  N 

is  the  total  number  of  pixels  in  the  image.  The  resulting  number  is  a  measure  of  the 
quality  of  the  beam  alignment. 

Table  6  summarizes  the  degree  of  misalignment  of  the  various  beams.  In  this 
calculation  N  is  57,600  and  the  pixel  intensity  differences  are  those  seen  in  Figure  42- 
Figure  44.  The  results  in  the  table  reinforce  what  was  clear  based  on  the  images.  The 
mean  square  deviation  described  indicates  that  the  Stokes  beams  overlap  nearly  7  times 
better  than  the  transmitted  and  input  beams.  Since  the  output  power  of  a  single  laser  or 
amplifier  is  inherently  limited  by  beam  degradation,  thermal  loading,  or  pump  coupling, 
this  method  of  incoherent  beam  combining  can  be  used  to  scale  diffraction-limited 
powers  beyond  the  limit  that  a  single  emitter  can  produce. 


Table  6:  Quantitative  measure  of  beam  overlap 


Pump 

SBS 

Transmitted 

7.07 

1.05 

6.86 

Relative  Misalignment 

6.73 

1 

6.53 

5.3.3.  Spatial  coherence 
Experimental  Setup 

The  previous  experiments  into  SBS  incoherent  combining  were  completed  with 
relatively  low  power  diode  lasers,  thus  the  fiber  diameter  had  to  be  kept  small  to  observe 
SBS.  At  the  time  of  the  experiment,  those  were  the  only  appropriate  lasers  for  the 
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investigation.  However,  two  single-frequeney  Lightwave  Eleetronics  Nd:YAG  lasers 
were  aequired  about  18  months  later.  Eaeh  laser  eould  be  tuned  several  GHz  thus  their 
frequeneies  eould  easily  be  separated  by  more  than  the  Brillouin  gain  bandwidth.  This 
made  these  sourees  perfect  for  further  investigation  into  incoherent  beam  combining. 

Although  the  previous  experiments  investigated  the  spatial  overlap  of  the  pump  and 
Stokes  beams,  they  do  not  necessarily  prove  spatial  coherence.  A  Gaussian-like  Stokes 
output  could  be  observed  if  the  beam  traveled  in  a  superposition  of  fiber  modes  or  if  the 
beam  spatial  pattern  fluctuated  fast  relative  to  the  CCD  frame  rate.  This  Gaussian-like 
output  would  not  have  the  high  spatial  quality  of  the  EPqi  mode.  The  following 
experiment  was  designed  to  further  strengthen  the  argument  that  not  only  do  the  Stokes 
beams  overlap,  but  they  are  also  spatially  coherent. 

A  schematic  of  this  experimental  setup  is  shown  in  Eigure  45.  Two  Eightwave- 
Electronics  Nd;YAG  lasers,  operating  at  slightly  different  frequencies,  were  nearly 
spatially  combined  with  a  polarization  beam  splitter.  Both  beams  were  then  coupled  into 
a  4.4  km,  graded-index,  50-micron  fiber.  The  portion  of  the  scattered  beam  that  reflects 
off  the  back  surface  of  an  uncoated  plane-parallel-plate  was  split  using  a  50/50  beam 
splitter.  The  two  beams  then  traveled  two  paths  of  equal  length  to  a  second  50/50  beam 
splitter  where  they  were  nearly  spatially  re-combined.  The  interference  pattern  generated 
was  then  imaged  with  a  CCD  camera  at  several  points  as  the  angle  of  the  mirror  indicated 
in  the  interferometer  was  adjusted. 
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Results 

Figure  46  shows  the  interference  pattern  generated  by  the  two  paths  at  several  mirror 
angles.  As  the  mirror  is  rotated  and  the  two  paths  overlap,  a  very  distinct,  linear 
interference  pattern  becomes  visible.  These  interference  lines  are  visible  as  the  moving 
beam  is  swept  through  both  sides  of  the  stationary  beam,  which  is  strong  evidence  that 
the  beams  are  spatially  coherent.  A  spatially  incoherent  beam  may  a  produce  linear 
interference  pattern  when  the  beams  overlap  perfectly.  However,  the  interference  would 
wash  out  or  become  distorted  with  any  spatial  misalignment  between  the  two  paths. 


Since  the  separation  between  two  adjacent  constructive  interference  points  is 
frequency  dependent,  it  was  critical  in  this  experiment  that  the  pump  beams  operated  at 
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Figure  46:  Interference  patterns  of  the  Stokes  beam  demonstrating  its  spatial  coherence. 


nearly  the  same  frequeney.  If  they  were  not,  the  interference  lines  generated  by  the  two 
frequency  components  of  the  Stokes  beam  would  become  much  more  difficult  to  interpret 
even  if  both  beams  propagated  in  the  LPoi  mode. 


5.3.4.  Conclusions 

In  this  chapter  a  method  for  incoherently  combining  multiple  laser  beams  into  a  single 
beam  of  high  spatial  coherence  through  stimulated  Brillouin  scattering  in  a  long  multi- 
mode  liber  is  examined.  The  experiments  demonstrate  the  independence  of  the  Stokes 
wave  excitation  for  pump  beams  with  different  wavelengths.  The  quality  of  beam 
overlap  has  been  quantified  using  the  mean  square  error  between  the  beam  profiles. 
Using  this  measure,  it  was  shown  that  the  two  multi-mode  input  pump  beams  that  were 
not  incident  coaxially  were  converted  into  a  beam  of  single  spatial  mode.  In  addition,  it 
was  shown  that  the  two  Stokes  beams  resulting  from  the  separate  lasers  had  essentially 
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identical  spatial  intensity  profiles.  Finally,  observing  the  interference  between  different 
portions  of  the  beam  further  supported  the  spatial  coherence  of  the  beams. 
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6,  SBS  Beam  Cleanup  and  Combining  Limitations 


There  are  many  potential  limitations  to  the  sealability  of  the  beam  eombining  and 
beam  eleanup  experiments  described  in  Chapters  4  and  5.  Some  of  these  boundaries  and 
schemes  to  overcome  them  are  described  here.  In  Section  6.1,  a  scheme  to  increase  the 
number  of  beams  that  can  be  combined  in  a  single  is  presented.  This  is  followed  by  a 
discussion  of  the  second  order  Stokes  threshold,  and  the  chapter  is  concluded  in  Section 
6.2  with  a  description  of  stimulated  Raman  scattering  beam  cleanup. 

6,1.  Multiple  Beam  Combining 

The  beam  combining  schemes  of  the  previous  chapter  were  limited  by  the  output 
coupling  of  the  Stokes  beam.  In  all  of  the  cases  described,  the  useful  output  was 
collected  by  the  Fresnel  reflection  of  an  uncoated  plane -parallel-plate.  This  is 
approximately  4%  of  the  total  Stokes  power  and  makes  the  true  efficiency  of  the  beam 
combining  very  low.  Coatings  on  the  surface  of  the  flat  could  increase  this  reflection 
coefficient,  but  since  the  Stokes  shift  is  so  small,  this  will  come  at  the  expense  of  the 
pump  power  and  nothing  is  gained.  Neither  can  polarization  be  used  since  the  pump 
beams  are  orthogonally  polarized.  If  incoherent  beam  combining  is  used,  the  beams  can 
be  nearly  spatially  combined  with  a  dichroic  mirror  and  polarization  can  be  used  to 
extract  the  Stokes  beam.  However,  scaling  this  to  more  than  a  couple  lasers  would  be 
extremely  complicated. 

A  proposed  solution  to  this  problem  is  to  pump  the  edges  of  the  fiber  and  extract  the 
Stokes  beam  from  the  center  of  the  fiber.  Although  this  spatial  output  coupling  was  not 
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completed,  the  ideas  examined  during  the  eourse  this  researeh  are  presented  here.  In  eaeh 
of  these  designs,  the  ideas  are  readily  expanded  to  more  pump  beams. 

The  most  eoneeptually  simple  experiment  along  these  lines  is  to  use  four  independent 
mirrors  equally  spaeed  about  a  eentral  hole  (see  Figure  47).  Four  pump  beams  ineident 
upon  these  mirrors  eould  then  be  eoupled  into  the  fiber  at  relatively  high  angles  while  the 
Stokes  beam  would  eseape  out  the  eentral  hole.  In  praetiee  however,  this  experiment 
would  be  very  diffieult  to  implement.  The  multimode  fiber  used  for  the  majority  of  the 
beam  eleanup  and  eombining  experiments  has  a  eore  diameter  of  50-miorons.  Simply 
imaging  the  mirror  assembly,  diameter  =  J  ,  on  the  fiber  faeet  with  a  single  lens,  foeal 
length  =  fi ,  would  require  that  the  distanee  between  the  mirror  and  the  lens  be 
approximately 


s 


fid 

50x10  ■ 


(6.1) 


A  lens  with  a  short  foeal  length  and  a  small  mirror  assembly  is  preferable  to  reduee  the 
distanee  from  the  mirrors  to  the  lens.  Assuming  a  standard  mieroseope  objeetive  with  a 
numerieal  aperture  that  matehes  the  numerieal  aperture  of  the  fiber  {NA  =  0.201),  the 
foeal  length  is  expeeted  to  be  approximately  16mm.  Using  a  relatively  eompaet  mirror 
assembly  sueh  that  d  =  25mm ,  results  in  5  =  8m  .  Sueh  a  large  distanee  between  the 
mirrors  and  the  lens  will  result  in  a  system  that  is  very  sensitive  to  vibration  and 
unreliable. 
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Figure  47 :  Multi-beam  coupling  for  beam  combining. 


A  more  elegant  solution  to  the  problem  would  be  to  use  a  short  bundle  of  fibers.  In 
this  design,  six  fibers  eould  easily  be  arranged  about  a  eentral  output  fiber,  see  Figure  48. 
Sealing  to  larger  bundles,  sueh  as  adding  another  ring  of  12  input  fibers,  would  be  a 
relatively  simple  task.  Immediately  after  the  fiber,  eaeh  of  the  diverging  waves  would  be 
collimated  with  a  lens  array  or  fiber  pigtail  collimators  and  the  overall  beam  would  be 
focused  into  the  long,  multi-mode  fiber. 


Figure  48:  Fiber  bundle  multiple  beam  combining  schematic  diagram. 


There  are  several  constraints  on  the  pump  beam  and  fiber  array  required  to  ensure 
efficient  beam  combining.  First,  each  pump  beam  must  have  a  spot  size  that  fits  within 
the  active  fiber  core.  Second,  the  numerical  aperture  of  the  combined  pump  beams  must 
be  less  than  the  numerical  aperture  of  the  fiber.  Finally,  the  spot  size  of  the  reflected 
Stokes  beam  must  be  efficiently  coupled  into  the  output  fiber.  The  effect  of  each  of  these 
requirements  is  analyzed  below. 

Assuming  each  beam  exiting  the  fiber  bundle  has  a  Gaussian  intensity  distribution 
and  is  collimated  by  the  lens  array  such  that  the  beam  diameter  is  ,  each  beam  will  be 
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2^  f 

focused  by  the  eoupling  lens,  foeal  length  =  /  ,  to  a  spot  size,  w  =  — — .  Any  aberration 


^  Tid 


on  the  pump  beam  will  inerease  this  spotsize  by  a  eonstant  faetor.  Since  this  spot  size 
must  fit  within  the  fiber  eore,  the  first  eonstraint  on  the  pump  beams  beeomes 


2  2 

w  =  M  w  <  —a  . 
^  3 


(6.2) 


where  is  the  beam  quality  faetor  reviewed  by  Siegman  in  [51]  and  a  is  the  radius  of 
the  fiber  eore.  The  faetor  of  2/3  is  ineluded  in  the  inequality  to  ensure  that  more  than 
99%  of  the  pump  beam  is  enelosed  by  the  fiber  eore. 

The  seeond  eonstraint  on  the  set  of  pump  beams  is  that  it  has  a  numerieal  aperture, 
NA^ ,  less  than  that  of  the  aetive  fiber,  NA  .  The  effeetive  numerieal  aperture  of  the  pump 

depends  on  the  total  diameter  of  the  pump  ring  after  the  collimating  lens,  ,  and  the 
foeal  length  of  the  eoupling  lens. 

NA=-^<NA  (6.3) 

P  2/ 

Additionally,  the  baekward  traveling  Stokes  beam  must  be  effieiently  eoupled  into  the 
output  fiber.  Assuming  a  graded  index  fiber,  the  LPoi  mode  is  elosely  approximated  by  a 

Gaussian  beam  with  a  spot  size  of  /  [52].  If  the  waist  is  loeated  at  the  front  fiber 

V  tiNA 

faeet,  the  spot  size  just  after  the  eoupling  lens  is 
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(6.4) 
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To  reach  the  approximate  solution  it  is  assumed  that 


(NAf^ 


»  1 .  The  output  coupler 


\  a  J 

must  be  at  least  1 .5  times  the  size  of  the  Stokes  beam  to  collect  more  than  99%  of  the 
power.  As  can  be  seen  by  Figure  49,  the  total  diameter  of  the  ring  of  pump  beams  must 
therefore  be 
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(6.5) 


Plugging  this  into  Equation  (6.3)  and  combining  the  result  with  Equation  (6.2)  yields 
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Thus 
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Solving  for  results  in 


M"< 


2 

6 


2  NAna  iNAna  ^ 

V  ;l  ^ 


(6.8) 


Using  this  equation,  the  maximum  value  ean  be  ealeulated  as  a  funetion  of  the  fiber 
parameters,  NA ,  a  ,  and  the  wavelength  of  the  light  used,  A  .  The  multi-mode  fiber  used 
in  the  beam  eombining  experiments  of  the  previous  ehapter  had  a  numerieal  aperture  of 
0.206  and  a  fiber  radius  of  50  mierons.  If  pumped  at  1.064  mierons,  the  maximum  value 
of  is  2. 1 .  Thus  the  fiber  bundle  must  deliver  good  quality  pump  beams  for  effieient 
beam  eombining.  This  eonstraint  eould  be  signifieantly  relaxed  using  a  fiber  with  a 
larger  eore  diameter,  higher  numerieal  aperture,  or  using  a  shorter  wavelength. 

Equation  (6.6)  ean  be  used  to  determine  the  best  set  of  values  for  dp  and  / . 

However,  geometrieal  eonstraints  will  also  have  to  be  taken  into  aeeount.  Convenienee 
may  dietate  that  eaeh  fiber  in  the  bundle  be  equally  spaeed,  whieh  would  foree  the 
spaeing  between  fiber  eores  to  be  the  greater  of  d ^  and  .  This  will  inerease  the 

eonstraint  on  and  for  some  fiber  parameters  and  pump  wavelengths  it  may  not  be 
possible  to  obtain  effieient  beam  eombining. 

Newer  pumping  sehemes  may  also  be  employed  to  simplify  beam  eombining.  Lueent 
Teehnologies  has  developed  a  tapered  fiber  bundle  for  use  in  pumping  fiber  lasers  and 
amplifiers  [53].  A  sehematie  diagram  of  the  fiber  is  shown  in  Figure  50.  In  eontrast  to 
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the  design  above,  the  pump  fibers  are  physieally  attaehed  to  the  aetive  fiber,  thus  the 
system  would  be  very  robust. 


- 


Figure  50:  Tapered  fiber  bundle  designed  by  Lucent  Technologies  [53]. 


6,2.  Second  Order  Threshold 

The  seeond-order  SBS  proeess  eould  potentially  limit  the  maximum  power  that  ean 
be  eombined  or  eleaned  in  a  single  fiber.  Above  the  seeond-order  threshold,  any  inerease 
in  the  pump  power  will  eontribute  primarily  to  the  seeond-order  beam,  resulting  in 
deereasing  effieieney  with  inereasing  pump  power.  Other  proeesses  that  rely  upon  SBS 
in  libers  sueh  as  phase  eonjugation  [54]  will  suffer  from  similar  limitations.  Thus, 
knowledge  of  this  threshold  is  eritieal  to  understanding  the  potential  benefits  and  limits  of 
these  applieations.  In  this  seetion,  the  seeond-order  Stokes  threshold  is  measured  and 
analyzed. 

This  investigation  was  eompleted  in  two  parts.  First,  a  long  fiber  with  a  small 
diameter  was  used  to  measure  the  threshold  experimentally.  Seeond,  a  set  of  the  eoupled 
differential  equations  for  the  pump,  first  Stokes,  and  seeond-Stokes  beams  was 
numerieally  solved  using  appropriate  eoupling  terms  and  fit  to  the  experimental  data  for 
validation.  The  eoupled  equations  eould  then  be  used  to  determine  the  seeond-Stokes 
threshold  for  fibers  of  arbitrary  length  and  size. 
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6.2.1.  Experimental  Setup 

The  experiment  was  performed  using  the  optieal  setup  sehematieally  shown  in  Figure 
51.  A  4.4  km  Corning  SMF-28  fiber  was  used  as  the  SBS  medium.  It  was  pumped  at 
1.064  pm  with  a  diode-pumped,  single-frequeney  Nd;YAG  laser  produeing  up  to 
600  mW.  After  aeeounting  for  two  passes  through  an  amplifier  as  well  as  eoupling 
losses,  more  than  800  mW  was  measured  in  the  fiber.  As  the  pump  power  was  varied, 
the  power  in  the  refleeted  Stokes  beam  was  measured  and  the  relative  power  transmitted 
through  the  fiber  was  monitored  with  a  seanning  Fabry-Perot  speetrum  analyzer.  The 
speetrum  analyzer  output  was  used  to  determine  seeond-order  threshold.  Above 
threshold,  the  output  elearly  showed  three  distinet  peaks  (Figure  52,  inset).  The  height  of 
the  peak  furthest  to  the  left  indieates  the  relative  intensity  of  the  transmitted  residual 
pump  beam  that  was  not  absorbed  or  eonverted  into  the  Stokes  beam.  The  seeond  peak 
shows  the  portion  of  the  first  Stokes  beam  that  was  internally  refleeted  off  the  front  fiber 
faeet  and  then  traveled  forward  through  the  liber  to  be  deteeted  by  the  Fabry-Perot.  The 
third  peak  is  due  to  the  seeond-order  Stokes  beam.  The  Stokes  shift  of  16.1  GHz  is 
nearly  twiee  the  free  speetral  range  of  the  Fabry-Perot  resonator,  8  GHz,  resulting  in  the 
near  overlap  of  the  three  peaks. 

The  height  of  these  peaks  as  a  funetion  of  the  input  pump  power  is  shown  in  Figure 
53(a-e).  These  plots  reveal  a  strong  eoupling  between  the  transmitted  beams.  At  first 
order  threshold,  1 8  mW,  the  growth  of  the  transmitted  pump  beam  was  strongly 
attenuated  and  the  forward  propagating  first  order  Stokes  beam  began  to  grow  linearly. 
Similarly,  at  the  seeond  order  threshold,  273  mW,  the  growth  of  the  forward  propagating 
Stokes  beam  was  suppressed  and  the  seeond  order  Stokes  beam  began  to  inerease.  This 
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Figure  5 1 :  Schematic  of  the  experimental  setup  used  to  measure  the  second  order  Stokes 

threshold. 


trend  was  expected  since  power  is  transferred  to  lower  energies  through  SBS.  However, 
at  the  same  point,  the  transmitted  pump  beam  also  began  to  increase.  This  behavior 
indicated  that  four-wave  mixing  (FWM)  between  the  forward  propagating  Stokes  beam, 
the  pump  beam  and  the  second-order  Stokes  beam  was  responsible  for  part  of  the  transfer 
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of  energy  to  the  pump  and  seeond-order  beam.  At  even  higher  powers,  -500  mW,  the 
transmitted  pump  and  seeond-order  beams  tapered  off.  This  behavior  is  found  to 
eoineide  with  the  appearanee  of  a  third-order  Stokes  beam  in  the  transmitted  speetrum 
and  is  also  believed  to  be  due  to  four- wave-mixing. 

The  power  in  the  baekward  traveling  Stokes  beam  is  plotted  as  a  funetion  of  the  pump 
power  in  Figure  53(d).  This  further  reinforees  the  elaim  that  four-wave  mixing  is  at  least 
partly  responsible  for  the  growth  of  the  seeond-order  Stokes  beam.  The  slope  effieieney 
below  seeond  order  threshold  was  found  to  be  88  ±  1% .  Above  seeond  order  threshold 
but  below  the  appearanee  of  the  third-order  Stokes  beam,  the  slope  effieieney  was  found 
to  be  8 1  ±  2%  .  Sueh  a  small  ehange  in  slope  effieieney  reveals  that  the  eoupling  between 
the  baekward  traveling  Stokes  beam  and  the  seeond-order  Stokes  beam  is  at  best  very 
weak. 

To  test  this  theory,  the  normal  eleave  at  the  front  liber  faeet  was  replaeed  by  an  angle 
eleave  and  the  transmitted  and  refleeted  beams  were  monitored  in  the  same  way  as 
before.  This  data  is  shown  in  Figure  54.  The  angle  eleave  prevented  any  of  the  Stokes 
beam  that  was  internally  refleeted  by  the  front  fiber  faeet  from  being  guided  by  the  fiber. 
Thus  Figure  54(b)  shows  no  transmitted  Stokes  beam.  Additionaly,  without  the  forward 
propagating  Stokes  beam  in  the  fiber,  FWM  was  suppressed,  eliminating  the  seeond  order 
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CoupledPump  Power  (mW) 
(c) 


CoupledPump  Power  (mW) 
(d) 


]  p  53:  (a-c)  Power  transmitted  through  the  SMF-28  fiber  as  a  funetion  of  the  input 
l  .~  jpower  at  the  pump,  Stokes  and  second  order  Stokes  frequencies.  The  solid  lines 
indicate  the  best  theoretical  fit.  The  dashed  lines  show  the  best  fit  ignoring  either 
second  order  SBS  contributions  or  four  wave  mixing  effects,  (d)  The  first  Stokes 
power  as  a  function  of  input  power. 


Stokes  beam.  Unfortunately,  the  angle  eleave  also  redueed  the  coupling  efficiency  into 
the  fiber  and  the  maximum  power  coupled  into  the  fiber  was  only  500mW  in  the 
transmitted  spectrum.  This  is  nearly  twice  the  original  second-order  threshold.  However, 
no  second  order  Stokes  beam  was  observed.  Unfortunately,  the  slope  efficiency  was  also 
reduced  to  66  ±  1% .  This  reduction  from  the  high  slope  efficiency  of  the  normal  cleave 
data  is  believed  to  be  caused  by  a  reduction  in  the  overlap  of  the  pump  beam  and  the  first 
order  Stokes  beam  and  not  a  direct  result  of  the  angle  cleave. 
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CoupledPump  Power  (mW)  CoupledPumpPowa-(mW) 

(c)  (d) 

54:  (a-c)  Power  transmitted  through  the  SMF-28  fiber  as  a  funetion  of  the  input 
liwer  at  the  pump,  Stokes  and  second  order  Stokes  frequencies  when  the  front  fiber 
facet  was  cleaved  at  an  angle.  The  solid  lines  indicate  the  best  theoretical  fit.  The 
dashed  lines  show  the  best  fit  ignoring  either  second  order  SBS  contributions  or 
four  wave  mixing  effects,  (d)  The  first  Stokes  power  as  a  function  of  input  power. 


6.2.2.  Modeling 

The  plane  wave  solution  to  Maxwell’s  equations  including  the  effects  of  SBS,  four- 
wave  mixing  and  absorption  on  the  coupling  between  the  pump,  ,  backward  traveling 

Stokes,  ,  forward  traveling  Stokes,  A^^ ,  and  second  order  Stokes,  A^.^ ,  electric  field 
amplitudes  is  given  by 


^  =  -gs^aCnA^  |4  f  +  ils^n^canA/ A J  exp  [z  (2^^  - 

^  4  =  -gs^oCnA^  |4  f  +  gsi^oCnA^i  |4  f  +  y  4 
^  4/  =  '4^o«2®«4/442  exp  [z  (4  +  ^2  -  )  z]  -  y  4^ 

^  42  =  gsi^oCnA^i  |4  r  +  «2^o'^2®«4/"4*  [''  (24/  -  4  )  ^]  ■ 


(6.9) 
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In  these  equations,  So,  ,  and  Uj  are  vaeuum  permittivity,  the  Brillouin  gain,  and 

nonlinear  index,  respeetively.  Other  nonlinear  eoupling  terms  sueh  as  self-  and  eross- 
phase  modulation  [3]  did  not  affeet  the  results  in  a  meaningful  way  and  were  therefore 
negleeted. 

To  aeeount  for  the  faet  that  the  fields  travel  in  fiber  modes  and  thus  eannot  be  plane 
waves,  the  gain  eoeffieients  of  the  various  eomponents  in  Equation  (6.9)  must  be 
adjusted.  Additionally,  sinee  the  depletion  of  the  baekward  traveling  Stokes  beam  due  to 
seeond  order  SBS  was  relatively  minor,  this  term  was  negleeted  in  the  equations.  Stokes 
depletion  was  partially  aeeounted  for  with  the  initial  eonditions  that  will  be  diseussed 
later.  Finally,  beeause  the  SBS  Stokes  shift  is  so  small,  the  phase  mismateh  was  assumed 
to  be  negligible.  Thus  Equation  (6.9)  beeomes 


^  4  =  -^SBSs  gB^oCn4  |4  f  +  f  4 

■^4/  ~  442  “^4/ 

~^42  ~  ^SBS2Sb2^0^^42\4\  4  ~  ^ 


(6.10) 


The  magnitude  of  eaeh  of  these  fields  at  the  front  fiber  faeet  was  used  as  the  initial 
eondition  to  solve  these  eoupled  equations  numerieally.  The  Stokes  beam  amplitude  is 
initialized  based  upon  the  measured  refieeted  Stokes  beam.  Below  first-order  threshold, 
the  amplitude  is  zero.  Between  the  first-  and  seeond-order  thresholds,  the  Stokes 

intensity  is  rj^  “4i)  >  where  is  the  measured  slope  effieieney  of  this  region,  88%, 
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Ip  is  the  intensity  of  the  pump  beam  and  is  the  first-order  Stokes  threshold.  Above 


second-order  threshold,  ’  the  Stokes  beam  intensity  was  initialized  to 

V\  {^thi  ~  ^th\ )  +  ^2  “  ^thi )  ’  where  rj^  is  the  slope  efficiency  of  this  region.  Based  solely 

on  the  reflected  Stokes  beam,  this  efficiency  would  be  81%.  However,  this  number  was 
used  as  a  parameter  to  improve  the  fit  to  the  transmitted  data  resulting  in  //j  =  84% .  The 
forward  propagating  Stokes  beam  was  initialized  by  assuming  that  3.6%  of  the  Stokes 
beam  was  internally  reflected  at  the  front  fiber  facet,  which  was  all  coupled  back  into  the 
fiber.  Finally,  the  second-order  Stokes  beam  was  initialized  by  injecting  a  single 
spontaneous  photon  per  mode  at  the  front  liber  facet  [18]. 

Using  Equations  (6.10)  and  the  appropriate  initial  conditions  on  the  fields  at  the  front 
liber  facet,  the  standard  numerical  differential  equation  solver  in  Mathematica  generated 
the  expected  transmitted  powers.  Along  with  the  slope  efficiency  of  the  Stokes  beam 
above  second  order  threshold,  were  adjusted  to  minimize 

the  mean  square  error  between  the  relative  transmitted  power  in  each  frequency 
component  and  the  predicted  transmission  based  upon  Equation  (6.10).  The  fit  is  based 
solely  on  the  region  below  the  appearance  of  the  third  order  frequency  component, 

<500  mW  where  the  model  is  expected  to  perform  reasonably  well.  The  best  fit  values 
are  summarized  in  Table  7  and  the  resulting  curve  is  shown  by  the  solid  lines  in  Eigure 
53. 


Table  7:  Fit  parameters  for  second  order  Stokes  threshold  model. 


Vi 

^SBSp 

VsBSs 

VsBSl 

^FWM 

Eit  Value 

0.84 

0.67 

0.58 

2.7 

1.5 
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The  effect  that  each  of  these  parameters  has  on  the  overall  transmission  curves  is 
rather  complicated.  However,  the  numerical  factors  can  be  informative.  The  SBS  gain  in 
a  non-polarization  preserving  fiber  is  expected  to  be  a  factor  of  1 .5  less  than  that  of  a 
polarization  preserving  fiber  [55].  This  explains  the  factor  of  0.67  for  . 

The  difference  between  and  is  responsible  for  the  portion  of  the  pump 

beam  that  is  Brillouin  reflected  into  a  radiation  mode  and  is  therefore  not  guided  by  the 
fiber.  Although  not  modeled  in  this  manner,  the  presence  and  excitation  of  these  modes 
is  explained  by  Garth  et  al.  [56].  The  ratio  of  these  numbers  is  closely  related  to  the 
slope  efficiency  of  the  first-order  Stokes  beam  ^sbss/^sbsp  ~  0-86  . 

The  second  order  Brillouin  gain  is  enhanced  in  this  model  by  a  factor  of  •  This 

factor  is  due  in  part  to  the  increased  effective  gain  that  results  from  the  enhanced  overlap 
of  the  Stokes  and  second  order  beams  caused  by  beam  cleanup.  Since  in  Chapter  4  it  was 
demonstrated  that  an  aberrated  pump  beam  excited  a  Stokes  beam  in  an  LPoi  fiber  mode, 
it  is  reasonable  to  assume  that  the  second  order  Stokes  beam  would  also  propagate  in  the 
LPoi  mode.  Thus  each  of  these  modes  is  confined  to  the  core  more  strongly  than  an 
aberrated  beam  and  the  effective  intensity  is  increased  over  that  of  a  uniformly 
illuminated  fiber  core.  This  accounts  for  part  of  the  enhanced  Brillouin  gain.  However  it 
is  unlikely  that  it  completely  explains  the  magnitude  of  this  factor.  Another  contribution 
involves  the  neglected  depletion  of  the  Stokes  beam.  Depletion  of  the  Stokes  wave  due 
to  the  second  order  beam  would  slow  the  growth  of  the  Stokes  beam  in  the  backward 
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direction.  Because  of  the  forced  condition  on  the  Stokes  power  at  the  front  fiber  facet, 
this  elevates  the  intensity  of  the  Stokes  beam  over  a  longer  portion  of  the  liber,  thus 
increasing  the  fiber’s  effective  length  and  improving  conversion  to  the  second  order 
beam. 

Similar  arguments  could  be  made  to  explain  the  enhancement  to  the  four-wave¬ 
mixing  gain.  However,  it  is  complicated  by  the  inclusion  of  three  field  distributions  and 
the  fact  that  the  different  terms  were  not  allowed  to  vary  independent  of  one  another. 
Initially,  these  were  allowed  to  independently  change  to  enhance  the  fit  between 
experiment  and  theory.  However,  the  benefit  of  the  two  extra  degrees  of  freedom  was 
minimal. 

To  model  the  effect  that  the  angle  cleave  had  on  the  transmitted  beams,  the  initial 
conditions  for  both  the  backward  and  forward  propagating  Stokes  beams  had  to  be 
modified.  The  initial  condition  for  the  backward  propagating  Stokes  beam  at  the  front 
fiber  facet  was  modeled  by  the  reflected  Stokes  power  shown  in  Figure  54(d).  Thus 
below  first  order  threshold,  the  intensity  was  set  equal  to  zero.  Above  threshold,  the 

slope  efficiency  was  constant  and  the  input  intensity  was  equal  to  ,  where 

//j  =  66%  .  On  the  other  hand,  the  forward  propagating  Stokes  beam  was  set  equal  to 
zero  at  the  front  fiber  facet  for  all  pump  powers,  accounting  for  the  assumption  that  none 
of  the  Stokes  beam  that  is  internally  reflected  at  the  front  liber  facet  is  coupled  back  into 
a  guided  mode.  The  results  of  these  calculations  are  shown  by  the  solid  lines  in  Figure  54 
and  demonstrate  a  strong  agreement  with  the  measured  data. 
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Since  second  order  threshold  without  the  forward  propagating  Stokes  beam  was  not 
observed,  the  threshold  for  this  condition  had  to  be  calculated  theoretically  with  no 
experimental  verification.  When  =  0  ,  Equation  (6.10)  becomes 

\^p f  +  W (6-11) 

oz  '  2 

OZ 

Which  can  be  written  in  terms  of  the  beam  intensities 


=-gBlpIs-(^Ip 

oz 

~^^s2  ~  gai^sl^s  ~^^s2 
OZ 


(6.12) 


Since  the  depletion  of  the  backward  propagating  Stokes  beam  was  ignored,  the  intensity 
distribution  of  the  second  order  beam  can  be  solved  explicitly  as  a  function  of  the  Stokes 
beam  intensity. 


^2(^)  =  ^2(0)exp 


gB2\ls{^)(i^ 
V  0 


-a 


(6.13) 


According  to  Smith  [18],  first-order  threshold  for  SBS  occurs  when  the  gain  in  the 
exponent  is  equal  to  21.  Assuming  this  numerical  factor  is  the  same  for  the  second-order 
threshold  results  in 
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P  =21 

^  sthl 


§B2^eff 


(6.14) 


where  P^,^2  is  the  Stokes  power  at  threshold,  is  the  effective  cross  sectional  area  of 


the  fiber,  and  the  effective  length,  ,  is  given  by 


(6.15) 


The  Stokes  intensity  distribution  can  be  solved  exactly  if  the  pump  and  Stokes 
absorption  is  ignored.  This  leads  to 


/,(0)exp{g,[/^(0)-/,(0)]z}-/,(0) 


(6.16) 


and  the  effective  length  is  analytically  solved  to  be 


Kff  - 


!.i0) 


if  (0) -/,(«)] +  ^  Ini 

§B 


exp 


{i®.f,(0)-i,(0)]}i,(0)-i.(0) 


MO)-l(O) 


.(6.17) 


With  absorption  taken  into  account,  the  Stokes  intensity  distribution  in  the  fiber  can 
be  estimated  by  numerically  solving  the  first  two  equations  of  (6. 11).  Similar  to  the 
analytical  solution  of  Equation  (6.16),  it  has  a  strong  dependence  on  the  pump  power 
coupled  into  the  fiber.  At  high  pump  powers,  the  Stokes  beam  grows  very  quickly  near 
the  front  of  the  fiber  resulting  in  a  short  effective  length.  At  low  pump  powers,  the 
Stokes  beam  has  a  slow  growth  in  the  backward  direction,  resulting  in  a  relatively  long 
effective  length.  The  relative  Stokes  power  distribution  in  the  fiber  along  with  the 
effective  length  of  the  fiber  for  various  pumping  conditions  are  plotted  in  Figure  55. 
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Fiber  Position  (m)  Pump  Input  Power  (mW) 

(a)  (b) 

Figure  55:  (a)  Relative  Stokes  power  distribution  in  a  fiber  for  a  20  mW  pump  beam 
(solid)  and  a  800  mW  pump  beam  (dashed),  (b)  Second-order  effective  length  as  a 
function  of  input  power  for  a  4.4  km  fiber. 

Using  this  for  the  effective  length  in  Equation  (6.14)  results  in  a  Stokes  power  at 
threshold  of  =  2.1  W  and  a  pump  power  at  threshold  of  2.4  W  .  If  the  Brillouin  gain 

coefficient  as  given  in  the  open  literature  is  used,  these  values  increase  significantly,  thus 
this  should  be  regarded  as  a  minimum  threshold. 

6.2.3.  Conclusions 

Even  in  the  presence  of  four-wave  mixing,  the  second  order  threshold  was  nearly 
fifteen  times  greater  than  the  first  order  threshold.  When  four-wave  mixing  was 
suppressed,  the  threshold  was  found  experimentally  to  be  at  least  27  and  theoretically  to 
be  more  than  130  times  the  first  order  threshold.  These  relatively  high  thresholds  are  due 
to  the  decreasing  effective  length  of  the  fiber,  ,  with  increasing  pump  power  as  seen 

by  the  second  Stokes  beam. 
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As  shown  in  Equation  (6.14),  the  second  order  threshold  will  proportionally  increase 
with  the  effective  area  of  the  fiber.  Thus  for  beam  combining,  cleanup,  and  phase 
conjugation  applications  that  require  a  large  fiber  core,  the  second  order  threshold  is  not 
expected  to  become  a  practical  limiting  factor.  The  useful  power  bandwidth  of  these 
devices  has  the  potential  to  stretch  from  the  first  order  threshold  to  the  second  order 

threshold,  which  was  calculated  to  be  43*^4^^  j^mW/ pm^^  in  the  absence  of  a  forward 

propagating  Stokes  beam.  In  the  50-micron  fiber  used  for  the  beam  cleanup  and 
combining  experiments  described  in  Chapters  4  and  5,  this  results  in  a  power  bandwidth 
of  more  than  80  kW.  In  almost  all  practical  cases,  it  can  be  concluded  that  second-order 
SBS  will  not  be  a  significant  limiting  factor  in  the  operational  dynamic  range  of  SBS 
fiber  beam  combiners. 

6,3.  SRS  Beam  Cleanup  and  Combining 

The  Brillouin  gain  bandwidth  is  very  narrow.  This  forces  the  pump  laser  to  be 
spectrally  pure  or  results  in  an  increase  in  the  Brillouin  threshold  [57,  58].  Either  of  these 
scenarios  can  be  problematic  for  practical  SBS  beam  cleanup/combining  schemes.  On 
the  other  hand,  the  Raman  gain  bandwidth  is  about  five  orders  of  magnitude  larger, 
eliminating  the  stringent  requirements  on  the  pump  lasers.  Spectrally  broad  beams,  such 
as  pulsed  sources  and  common  diode  lasers,  can  efficiently  excite  a  Stokes  beam. 
Unfortunately,  the  relaxed  requirements  come  with  a  price.  The  Raman  gain  coefficient 
is  about  two  orders  of  magnitude  smaller  than  the  Brillouin  gain  coefficient,  resulting  in  a 
greater  threshold.  Eor  extremely  high  power  applications  however,  the  threshold  is 
expected  to  be  a  small  percentage  of  the  overall  power  and  it  is  reasonable  that  stimulated 
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Raman  scattering  will  be  a  more  praetieal  nonlinear  meehanism  for  beam  eombining  and 
eleanup  than  SBS. 

Similar  to  stimulated  Brillouin  seattering,  SRS  has  been  used  for  phase  eonjugation 
[59,  60],  but  in  fibers  with  long  interaetion  lengths  the  phase  eonjugate  fidelity  is  lost  and 
the  resulting  Stokes  beam  is  generated  in  the  dominant  fiber  mode.  In  1992,  Chiang 
showed  that  when  SRS  was  exeited  in  a  liber,  fewer  fiber  modes  were  observed  in  the 
Stokes  beam  than  found  in  the  pump  beam  [25].  He  demonstrated  that  partieular  liber 
modes  eould  be  exeited  in  a  30m  fiber  by  varying  the  pumping  eonditions.  Longer  fiber 
lengths,  up  to  1  km,  resulted  in  a  mixture  of  up  to  two  fiber  modes  found  in  the  first  order 
Stokes  beam,  but  a  pure  LPqi  mode  found  in  eaeh  of  the  higher  order  modes. 

Raman  beam  eleanup  is  not  limited  to  fiber  based  systems.  Long  effeetive  interaetion 
lengths  ean  be  obtained  by  plaeing  a  Raman  aetive  material  in  a  eavity  resonant  at  the 
Stokes  frequeney.  These  Raman  lasers  have  theoretieally  and  experimentally  been  shown 
to  produee  a  high  quality  Stokes  beam  even  when  the  pump  is  a  low  quality  beam 
[61,  62].  SRS  has  also  been  used  to  amplify  a  high  spatial  quality  Stokes  beam  using  an 
aberrated  pump.  In  these  experiments,  the  spatial  profile  of  the  Stokes  beam  ean  be 
maintained  in  the  amplified  output  [63,  64]. 

This  seetion  deseribes  the  experimental  results  demonstrating  the  beam  eleanup 
eharaeteristies  of  stimulated  Raman  seattering  in  long  fibers.  The  experiment  is 
deseribed  in  Seetion  6.3.1  and  the  results  are  analyzed  and  explained  in  Seetion  6.3.2. 
Conelusions  are  made  in  Seetion  6.3.3. 
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6.3.1.  Experimental  Setup 

Stimulated  Raman  scattering  beam  cleanup  was  demonstrated  using  a  setup 
represented  by  the  schematic  shown  in  Figure  56.  A  frequency  doubled  Continuum 
Powerlite  Nd:YAG  laser  oscillator  was  used  to  pump  a  300  m  long,  50-micron  graded- 
index  fiber.  The  laser  produced  an  average  power  of  50  mW  delivered  in  10.5  ns  pulses 
at  30  Hz.  The  power  coupled  into  the  fiber  was  controlled  using  a  half-wave  plate  and  a 
polarizer.  A  lOx  microscope  objective  collimated  the  transmitted  beam. 


The  beam  was  then  focused  with  aim  focal  length  lens  and  a  long  pass  filter  was 
used  to  remove  the  pump  from  the  Stokes  beam.  The  evolution  of  the  beam  profile  was 
recorded  at  various  points  after  the  lens,  showing  the  spot  as  it  propagated  through  focus. 
Each  profile  was  normalized  to  have  an  area  of  one,  fit  to  a  Gaussian  intensity 


distribution,  — !-yexp 


and  the  spot  size  (half-width  at  the  l/ e^  point),  <j ,  was 


Tier  \  (j  j 

recorded.  This  process  was  repeated  with  the  laser  power  reduced  below  threshold  and 
the  long  pass  filter  removed,  showing  the  propagation  of  the  pump  beam  over  the  same 
region. 


The  transmitted  beam  was  then  characterized  spectrally.  First,  the  beam  spectrum 
was  analyzed  with  a  monochrometer  equipped  with  an  optical  multichannel  analyzer. 
Then  the  mirror  near  the  output  end  of  the  fiber  was  replaced  with  a  diffraction  grating 
and  the  back  fiber  facet  was  focused  directly  on  the  CCD.  In  this  manner  the  spatial 
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Figure  56:  Schematic  diagram  of  stimulated  Raman  scattering  beam  cleanup. 


profile  of  each  spectral  component  could  be  observed.  This  was  completed  at  various 
fiber  lengths  and  the  size  of  the  first  Stokes  component  measured  and  compared  to  the 
size  of  the  residual  transmitted  pump  beam. 

Fiber  lengths  ranging  from  50  m  to  300  m  at  25  m  increments  were  used  to  observe 
the  effect  of  fiber  length  on  the  output  spatial  characteristics  of  the  far- field  beam  profile. 
At  each  length,  65  successive  images  were  collected  under  four  separate  conditions.  The 
first  two  image  sets  showed  the  spectrally  dispersed  transmitted  beam.  These  pictures 
included  the  pump,  Stokes,  and  second  order  Stokes  beams  in  each  frame  and  were  used 
to  measure  the  diameter  of  each  spectral  component.  The  next  set  of  images  was 
collected  below  threshold,  verifying  that  even  without  the  excitation  of  the  Stokes  beam, 
the  pump  was  transmitted  in  multiple  fiber  modes.  Finally,  a  FleNe  laser  operating  at 
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543  nm  was  used  to  illuminate  the  fiber  end  and  another  set  of  images  was  eolleeted. 
This  set  showed  the  liber  eore  and  eladding  illuminated  uniformly  and  allowed  for  an 
aeeurate  ealibration  of  the  spotsize  measurements  made  with  the  CCD  eamera.  It  also 
ensured  that  the  liber  faeet  was  aeeurately  imaged  on  the  CCD. 


6.3.2.  Results 

The  spatial  quality  of  the  pump  and  Stokes  beams  was  measured  by  eomparing  their 
propagation  to  that  of  a  Gaussian  beam.  The  spot  size  of  a  Gaussian  beam  with  a  waist, 
Wq  ,  loeated  at  a  point,  ,  will  vary  with  distanee  aeeording  to  [52] 


w. 


(z)  = 


Wn 


1  + 


V  ) 


(6.18) 


If,  on  the  other  hand,  the  beam  is  not  a  perfeet  Gaussian,  the  spot  size  everywhere  will  be 
sealed  by  a  eonstant  factor,  M  [51] 

w(z)=MWg(z).  (6.19) 

While  this  factor  is  a  good  measure  of  beam  quality,  it  is  typically  reported  as  .  As 
can  be  found  by  a  manipulation  of  Equation  (6.19),  is  the  ratio  of  the  far- field 
divergence  of  a  real  beam  to  that  of  a  perfect  Gaussian  with  identical  waist  diameters. 

An  example  of  the  near-  and  far-field  images  of  the  pump  and  Stokes  beams  is  shown 
in  Figure  57.  Clearly  the  pump  near-held  image  is  much  larger  than  the  Stokes  near-held 
image.  If  both  of  these  beams  were  perfect  Gaussian  beams,  the  pump  beam  would  focus 
to  a  much  smaller  spot  than  the  Stokes  beam  as  predicted  by  Equation  (6.19).  The  fact 
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that  this  is  not  the  case  indicates  that  the  Stokes  beam  has  a  much  better  beam  quality 
than  the  pump. 


• 

Pump  Near-Field 

Pump  Far-Field 

• 

• 

Stokes  Near -Field 

_ 

Stokes  Far-Field 

57:  Near-  and  Far- field  images  of  the  pump  and  Stokes  beams  transmitted 
through  a  300  m  multimode  fiber. 


A  more  rigorous  analysis  can  be  completed  by  fitting  the  spot  size  of  the  real  beam  as 
a  function  of  position  to  Equation  (6.19).  The  spot  sizes  along  with  the  least  squares  best 
fit  are  shown  in  Figure  58.  This  analysis  is  complicated  by  the  fact  that  the  wavelength 
difference  between  the  first  and  last  observed  Stokes  order  is  nearly  100  nm  or  20%.  The 
conservative  approach  taken  here  assumed  all  Stokes  power  was  found  in  the  first  order 
beam  at  545  nm  thus  overestimating  the  value  of  .  This  analysis  results  in  an 
value  of  2.4  ±0.3  for  the  Stokes  beam  and  an  of  20.7  ±1.0  for  the  pump  beam.  The 
error  in  the  measured  values  takes  into  account  the  imperfect  fit  of  the  beam  size  to  the 
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measured  values.  It  negleets  the  error  in  the  measurement  of  the  pump  and  Stokes  beam 
diameters  and  the  error  indueed  by  assuming  that  the  Stokes  beam  is  composed  of  a 
single  frequency. 


Figure  59  shows  an  example  of  the  transmitted  spectrum.  Up  to  seven  Stokes  shifts 
were  observed,  each  separated  by  approximately  440  cm'\  Considerable  spectral 
broadening  was  noticed  with  each  higher  Stokes  shift.  The  bottom  portion  of  this  figure 
shows  the  spatial  characteristics  of  the  different  modes  as  imaged  after  the  diffraction 
grating.  The  spectral  broadening  of  successive  Stokes  shifts  makes  a  fit  to  a  two 
dimensional  Gaussian  beam  pointless,  but  the  beam  size  can  still  be  quantitatively 
analyzed  along  the  vertical  dimension. 

This  was  completed  for  far-field  images  of  the  pump  and  Stokes  beams  at  various 
fiber  lengths.  Each  image  was  examined  for  a  pump  and  Stokes  beam.  The  relative 
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Figure  59:  (Top)  The  spectrum  of  the  beam  transmitted  through  the  fiber.  (Bottom)  The 
spectrally  dispersed  near-field  image  of  the  transmitted  beam. 

power  in  each  spectral  component  was  measured  by  adding  the  intensity  in  each  pixel 
around  the  beam.  Images  in  which  the  Stokes  beam  had  less  than  25%  of  the  power  of 
the  pump  beam  were  discarded  to  eliminate  images  dominated  by  noise.  The  remaining 
images  were  evaluated  to  determine  the  diameter  of  the  pump  and  Stokes  beams.  The 
measurement  was  completed  by  integrating  the  power  along  the  entire  horizontal  axis  and 
integrating  up  and  down  from  the  center  of  the  beam  until  59%  of  the  power  was 
contained.  This  vertical  height  represents  the  full-width-half-maximum  if  the  beam 
profile  were  a  perfect  Gaussian. 

The  measured  far-field  beam  diameters  are  shown  in  Figure  60(a)  as  a  function  of  the 
fiber  length.  In  all  cases  the  Stokes  beam  is  considerably  smaller  than  the  pump  beam 
indicating  significant  beam  cleanup.  The  spot  size  of  the  far-field  image  is  the  Fourier 
transform  of  the  collimated  beam,  or  equivalently,  that  of  the  aberration  on  the  divergent 
spherical  wave  emanating  from  the  fiber  facet.  Thus  the  far-field  spot  size  is  a  measure 
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of  the  spatial  frequency  content  of  the  collimated  beam.  Consequently,  for  a  given  near¬ 
field  beam  size,  the  far-field  spot  size  is  a  meaningful  measure  of  the  beam  quality. 


Fiber  Length  (m) 
(a) 


■:  60:  (a)  Spot  size  of  the  pump  (circles)  and  first  order  Stokes  (squares)  beams  at 
e  fiber  facet.  (b,c)  Pump  (left)  and  Stokes  (right)  intensity  distributions  for  a  300 
m  and  75  m  fiber  respectively. 

At  short  fiber  lengths,  the  Stokes  beam  intensity  is  expected  to  better  match  the  pump 
intensity.  Although  there  is  no  indication  of  this  in  the  figure,  different  mode  shapes 
were  observed  in  the  Stokes  output  at  fiber  lengths  below  100  m.  Figure  60(b  and  c) 
shows  a  representative  close-up  of  the  pump  and  Stokes  beams  for  long  fibers  (-300  m) 
as  well  as  the  same  image  in  relatively  short  fibers  (~75  m). 

6.3.3.  Conclusions 

Since  a  real  beam  diverges  faster  than  a  Gaussian  with  the  same  waist  by  a  factor  of 
,  the  intensity  of  a  real  beam  drops  by  a  factor  of  compared  to  a  Gaussian  beam. 
This  demonstrates  how  important  beam  quality  is  to  the  propagation  of  high  intensity 
beams.  The  quality  of  the  input  beam  was  not  measured.  Flowever,  it  could  be  at  least  as 
aberrated  as  the  transmitted  pump  and  still  be  efficiently  coupled  into  the  fiber. 

Assuming  a  pump  with  this  poor  beam  quality  was  used,  the  results  of  the  SRS  beam 
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cleanup  experiment  show  that  the  Stokes  beam  eould  illuminate  a  target  with  an  intensity 
that  is  74  times  greater  than  that  expeeted  by  the  pump  beam  with  similar  power.  Of 
eourse  the  pump  to  Stokes  eonversion  is  not  100%  effieient.  However,  based  upon  the 
inerease  in  beam  quality,  the  effieieney  only  has  to  be  greater  than  1.3%  to  obtain  an 
inerease  in  the  intensity  on  target. 

The  results  shown  in  Figure  60(a)  illustrate  that  the  spot  size  of  the  Stokes  beam  does 
not  signifioantly  ehange  in  fibers  down  to  50  m  and  that  the  Stokes  beam  is  signifieantly 
smaller  than  the  pump  beam.  Fibers  longer  than  100  m  eonsistently  generated  Stokes 
beams  with  Gaussian  intensity  distributions.  This  provides  strong  evidenee  that 
stimulated  Raman  seattering  in  long  fibers  produees  a  elean  beam  similar  to  the  beam 
eleanup  seen  with  SBS.  However,  there  is  elearly  a  transition  point  between  long  and 
short  fibers  in  whieh  the  beam  eleanup  properties  are  lost  and  higher  order  modes  are 
exeited.  Nevertheless,  down  to  the  shortest  fiber  measured  (50  m),  the  Stokes  beam 
always  appeared  to  have  a  mueh  smaller  spot  size  than  the  pump  beam.  These 
observations  are  eonsistent  with  those  made  by  Chiang  in  whieh  several  fiber  modes  were 
observed  in  a  30  m  fiber  [25]. 

Forward  SRS,  as  used  in  this  experiment,  allows  the  use  of  short  laser  pulses.  But  it 
also  allows  four  wave  mixing  to  signifioantly  reduoe  the  seoond  order  Stokes  threshold. 
Although  this  did  not  appear  to  affeot  the  Stokes  beam  diameter,  it  does  broaden  the 
speotrum  of  the  Stokes  beam.  This  eould  possibly  be  reduoed  through  an  appropriate 
design  of  a  long  period  grating  or  an  absorbant  that  eould  eliminate  the  seoond  Stokes 
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beam  while  not  affecting  the  pump  or  first  order  Stokes  beams.  These  concepts  are 
explored  in  more  depth  in  the  chapters  on  SRS  reduction  in  fiber  lasers. 

Additionally,  the  forward  scattering  scheme  in  a  non-polarization  preserving  fiber 
provides  no  means  to  control  the  polarization.  Even  in  a  backward  scattering  design,  the 
greater  Stokes  shift  as  compared  to  SBS  increases  the  chromatic  dispersion  reducing  the 
maximum  fiber  length  that  could  be  used  while  maintaining  the  pump  polarization, 
although  this  point  requires  experimental  verification. 

It  can  be  concluded  from  the  experiments  described  above  that  stimulated  Raman 
scattering  could  be  used  as  a  substitute  for  SBS  in  an  operational  beam  cleanup  system. 
Although  it  introduces  its  own  set  of  problems,  this  would  significantly  reduce  the 
spectral  constraints  on  the  lasers  needed  in  a  beam  combining  or  beam  cleanup  design. 
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7.  SRS  Reduction  Using  Intercavity  Dopants 


lA.  Background 

Fiber  lasers  are  becoming  increasingly  popular.  The  active  region  is  so  small  that  it 
can  easily  handle  high  power  beams  without  significant  thermal  gradients  that  can  cause 
fracture  or  thermal  lensing.  The  low  loss  and  high  beam  confinement  can  also  reduce 
laser  threshold  allowing  for  novel  materials  with  low  gain  coefficients  to  be  used  in  laser 
devices  [11].  Using  carefully  chosen  fiber  geometries  can  also  ensure  that  high  spatial 
quality  beams  are  excited.  For  example,  double  clad  fibers  accommodate  highly 
aberrated  pump  beams  and  yet  ensure  single  mode  laser  beams.  In  these  fibers  the  inner 
core  is  doped  with  the  active  element,  often  ytterbium,  and  has  a  small  enough  diameter 
that  it  only  allows  the  LPqi  mode  to  be  transmitted  at  the  laser  wavelength  (see  Figure 
61).  The  first  cladding  layer  has  an  index  of  refraction  that  is  only  slightly  lower  than 
that  of  the  core.  This  region  is  used  to  transmit  the  pump  light.  It  is  large  enough  that  it 
supports  a  multimode  beam,  but  is  not  doped  with  the  active  laser  ion,  so  little  light  is 
absorbed  outside  of  the  core.  The  second  cladding  region  has  an  index  of  refraction  less 
than  that  of  the  first  cladding  region.  Its  main  purpose  is  to  guide  the  pump  light  down 
the  fiber. 


Figure  61 :  Double  clad  fiber 
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Although  this  complicated  fiber  geometry  ensures  that  the  laser  beam  oscillates  in  a 
high  spatial  quality,  the  single  mode  eore  is  very  small  resulting  in  high  intensities  and 
observable  nonlinear  effeets.  All  third  order  interactions  have  the  potential  to  be 
detrimental  to  a  fiber  laser  system,  but  stimulated  Raman  seattering  is  often  the  limiting 
nonlinear  effeet  [65].  Many  proeesses,  sueh  as  four-wave  mixing,  require  phase 
matehing  for  effieient  light  eonversion,  thus  in  most  situations,  unless  speeial  eare  is 
taken,  these  proeesses  will  be  insignifieant.  Stimulated  Brillouin  seattering,  on  the  other 
hand,  is  automatieally  phase  matched  in  the  baekward  direction.  It  even  has  a  relatively 
large  peak  gain  eoeffieient  eompared  to  SRS,  but  the  gain  bandwidth  is  very  narrow 
(~100  MHz)  and  typieal  fiber  lasers  are  speetrally  broad,  due  to  inhomogenous 
broadening  in  the  amorphous  glass  host  [11].  Figure  62  illustrates  the  typieal  speetral 
broadening  assoeiated  with  glass.  The  plot  shows  the  emission  speetra  for  Nd  in  YAG,  a 
erystal,  and  amorphous  SiOi.  In  eomparison  to  SBS,  SRS  has  a  very  broad  gain 
bandwidth  (~30  THz  [3]),  whieh  is  large  enough  to  eompensate  for  the  lower  peak  gain. 
The  entire  laser  speetrum  ean  work  together  to  pump  a  single  Raman  Stokes  beam 
eausing  SRS  to  typieally  have  a  lower  threshold  than  SBS. 

Sinee  fiber  lasers  are  typieally  not  extremely  long  (-'1-10  m),  it  takes  a  lot  of  power  to 
generate  SRS,  but  fiber  lasers  are  getting  more  powerful  every  year.  In  1997,  the  highest 
reported  cw  output  from  double-elad  fiber  lasers  was  35  watts.  In  1998  that  jumped  to  55 
watts  and  in  1999  a  fiber  laser  generating  110  watts  was  reported  [66].  At  power  levels 
this  high,  SRS  is  signifieant  even  though  the  fiber  is  short.  Lang  reported  significant 
output  in  the  Stokes  beam  at  power  levels  exeeeding  100  watts  when  presenting  details  of 
an  80-watt  fiber  laser  he  developed  [67]. 
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Figure  62:  Emission  Spectrum  of  Nd  in  a  crystal  and  in  glass  [11]. 

Pulsed  laser  systems  are  even  more  susceptible  to  stimulated  Raman  scattering. 
Typical  laser  radar  systems,  for  example,  only  require  a  few  watts  of  average  power,  but 
delivered  in  pulses  less  than  about  10  ns.  Even  at  pulse  rates  of  tens  of  kHz,  the  peak 
powers  observed  are  many  kilowatts.  Gapontsev  et  al.  [68]  have  been  investigating  this 
issue  with  very  good  success  in  a  master  oscillator  power  fiber  amplifier  configuration. 
They  were  able  to  reduce  nonlinear  effects  by  increasing  the  master  oscillator  spectral 
width  and  tapering  the  core  size  of  the  fiber  amplifier  such  that  as  the  pulse  increased  in 
power,  the  power  density  remained  low.  Yet  even  with  these  modifications  SRS  was  still 
noticeable  in  the  output  when  they  reached  their  maximum  peak  power  level  of  100  kW. 

The  negative  effects  of  SRS  are  clear  when  observing  how  it  depletes  the  pump  beam 
in  a  fiber.  The  coupling  between  the  Stokes  and  pump  beams  is  described  by  Equation 
(3.12),  repeated  here  for  simplicity. 
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Assuming  the  absorption  at  the  pump  and  Stokes  wavelengths  is  negligible  and  the 
Stokes  intensity  at  the  front  of  the  fiber  is  insignificant  relative  to  the  pump  intensity, 
IpQ  -  7^0  ~  7^0 ,  the  pump  beam  intensity  at  the  end  of  the  fiber,  z  =  L,  is 
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Even  a  relatively  small  initial  Stokes  intensity  (spontaneous  scattering  off  thermal 
phonons  for  example)  can  induce  significant  loss  in  high  power  beams. 

If  the  Stokes  beam  never  becomes  strong  enough  to  cause  significant  pump  depletion, 
the  solution  of  Equation  (7.1)  shows  that  the  Stokes  beam  intensity  distribution  in  the 
fiber  grows  exponentially. 


^.(2)  =  ^oexp 


^^(l-exp(-«^z))-«^z 


(7.3) 


This  is  only  valid  when  the  total  power  converted  into  the  Stokes  beam  is  small  relative 
to  the  pump.  Assuming  this  is  so.  Equation  (7.3)  can  be  treated  as  a  small  perturbation  to 
the  pump  beam  and  used  in  Equation  (7.1)  to  predict  the  effect  of  the  Stokes  generation 
on  the  pump  power.  If  the  absorption  at  the  pump  wavelength  is  negligible,  the  pump 
power  transmission  is 
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(7.4) 


r  =exp< 


Sr^sO^p 


l-exp[(g/^o-aJl_ 


From  Equation  (7.4),  it  becomes  apparent  that  as  the  attenuation  of  the  Stokes  beam 
increases,  the  transmission  of  the  pump  beam  increases.  This  is  particularly  obvious  as 
the  fiber  length  is  increased  to  infinity,  at  which  point  the  pump  transmission  becomes  a 
step  function,  transitioning  from  0%  when  the  Stokes  absorption  is  less  than  the  Raman 
gain,  gJpQ  <  ,  to  100%  when  the  Stokes  absorption  is  greater  than  the  Raman  gain. 

The  absorption  in  a  typical  fiber  is  shown  in  Figure  63.  This  shows  no  special 
absorption  features  near  1.12-1.15  microns  that  will  suppress  SRS  inNd  or  Yb  high 
power  lasers.  Novel  concepts,  however,  have  been  used  to  increase  the  loss  at  the  Stokes 
wavelength.  For  example,  Sylvestre  et  al.  used  a  second  laser  beam  to  induce  loss  at  the 
Stokes  wavelength  [69,  70].  When  the  second  pump  beam  was  frequency  tuned  to 
exactly  twice  the  Stokes  shift,  any  Stokes  radiation  generated  by  the  first  pump  beam  was 
quickly  converted  to  the  second  order  Stokes  beam,  increasing  the  amount  of  pump 
power  that  could  be  transmitted  through  the  fiber. 
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Figure  63:  Absorption  of  typical  optical  fiber  [13]. 

In  another  scheme,  stimulated  scattering  beyond  the  first  order  beam  was  suppressed 
by  bending  the  fiber  [71].  Bend  induced  loss  acts  as  a  short  pass  filter  with  the  cutoff 
wavelength  determined  by  the  bend  radius.  With  the  proper  bend  choice  the  second  order 
Stokes  beam  would  experience  high  attenuation  with  minimal  effect  on  the  first  order 
Stokes  and  the  pump  beam.  This  was  verified  experimentally  showing  that  the  threshold 
for  second  order  SRS  was  increased  by  greater  than  a  factor  of  two  and  the  power  in  the 
first  order  beam  was  increased  by  more  than  40%. 

Petrov  et  al.  modeled  and  demonstrated  second  order  SRS  suppression  using  a 
capillary  waveguide  [72].  A  capillary  waveguide  uses  a  liquid  core  with  a  glass  cladding. 
If  the  index  of  refraction  of  the  liquid  is  strongly  dispersive,  the  index  can  be  made  to  be 
higher  than  the  glass  cladding  for  the  pump  and  Stokes  beam,  but  lower  than  the  cladding 
for  the  second  order  Stokes  beam.  For  such  a  system,  the  second  order  beam  is  not 
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guided  and  will  have  a  very  high  threshold.  Experimentally,  a  30-35%  increase  in  the 
first  order  Stokes  beam  was  achieved  when  the  second  order  beam  was  suppressed. 

Another  method  that  has  shown  promise  in  suppressing  stimulated  Raman  scattering 
is  to  dope  the  nonlinear  material  with  an  element  that  has  an  absorption  line  near  the 
Stokes  order  that  must  be  suppressed.  Urata  et  al.  doped  praseodymium  into  a  lanthanum 
tungstate  crystal  [73].  Pr  has  an  absorption  feature  near  600  nm  which  is  very  close  to  the 
second  order  stokes  beam  generated  through  SRS  pumped  with  a  frequency  doubled 
Nd:YAG  laser  at  532  nm.  This  method  showed  that  the  second  order  Stokes  threshold 
increased  by  at  least  a  factor  of  two  without  a  significant  effect  on  the  first  order 
threshold. 

The  spectral  broadening  of  material  transitions  in  a  glass  fiber  can  complicate 
suppression  of  stimulated  Raman  scattering.  If  the  transition  is  not  sufficiently  far  away 
from  the  absorption  peak,  the  tail  in  the  absorption  line  will  reduce  transmission  of  the 
pump  beam,  negating  the  effects  of  SRS  suppression.  On  the  other  hand,  the  long 
interaction  lengths  in  a  fiber  relax  the  requirements  on  the  absorption  at  the  Stokes 
wavelength  to  significantly  affect  the  SRS  threshold. 

The  primary  interest  of  this  chapter  is  to  increase  the  power  that  can  be  generated  by 
an  ytterbium  or  a  neodymium  fiber  laser  by  limiting  the  adverse  affects  of  SRS.  Yb  and 
Nd  are  common  dopants  used  for  high  power  fiber  lasers  and  amplifiers  [74].  These 
lasers  operate  near  1  micron  with  a  peak  in  the  Stokes  gain  spectrum  that  is  shifted  about 
55  nm  to  longer  wavelengths.  They  are  pumped  near  920  or  980  nm  [10].  There  are  a 
couple  of  important  requirements  on  any  dopant  that  could  be  used  to  suppress  stimulated 


139 


Raman  scattering.  The  dopant  should  minimally  affeet  operation  of  the  laser.  This 
means  that  there  should  be  little  absorption  at  both  the  laser  pump  wavelength  and  the 
lasing  wavelength.  The  seeond  requirement  is  that  it  must  absorb  at  the  Stokes 
wavelength  to  inerease  the  SRS  threshold.  The  rest  of  this  ehapter  will  deseribe  the 
modeling  used  to  prediet  a  substantial  inerease  in  SRS  threshold  for  light  at  1.064  and 
1 .090  mierons  transmitted  through  a  silieate  glass  fiber  doped  with  Holmium. 

7,2.  Theory 

7.2.1.  Dopant  Characterization 

The  absorption  eharacteristics  of  holmium  in  silieate  glass  has  all  of  the  features 
required  for  suppression  of  the  first  order  Stokes  peak  assoeiated  with  neodymium  and 
ytterbium  doped  fiber  lasers.  The  ^Ig  -  ^le  transition  of  holmium  in  glass  is  eentered 
around  1.152  mierons  [75].  Although  this  is  not  perfectly  centered  on  the  Stokes  peak  for 
both  neodymium  and  ytterbium  lasers,  it  has  a  relatively  broad  tail  that  affects  the  Stokes 
beam  of  both.  Figure  64  shows  the  transmission  of  light  through  a  5  mm  silicate  glass 
window  doped  with  2  mol%  holmium.  The  inset  marks  the  loeation  of  the  pump  for  the 
fiber  laser,  laser  output,  and  the  first  order  Stokes  beam  for  an  Yb  doped  laser.  The 
window  was  not  polished  to  have  an  optieally  flat  surfaee,  thus  seattering  as  well  as  the 
Fresnel  refleetion  at  the  window  surfaees  aceount  for  part  of  the  deviation  of  the  pump 
and  laser  transmission  from  100%. 

The  output  of  a  neodymium  laser  is  eentered  near  1.06  mierons  with  a  first  Stokes 
gain  peak  at  1.117  pm.  Figure  65(a)  shows  the  holmium  absorption  speetrum  along  with 
the  Raman  gain  speetrum  showing  considerable  misalignment  of  the  two,  yet  the  material 
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Figure  64:  Transmission  of  a  5mm  thick  silicate  glass  window  doped  with  2  mol%  of 
Ho.  The  dashed  line  indicates  the  location  of  a  typical  laser  pump,  the  solid  line 
shows  the  wavelength  of  the  laser,  and  the  dotted  line  represents  the  first  order 

Stokes  beam. 


absorption  is  still  22.89  dB/(m  mol%)  greater  at  the  Stokes  wavelength  than  at  the  laser 
wavelength.  For  comparison,  a  standard  telecom  fiber  at  1 .064  pm  absorbs  less  than 
0.001  dB/m. 

More  effective  suppression  of  the  Stokes  radiation  may  be  expected  using  an  Yb 
laser.  Ytterbium  lasers  operate  from  1.06  pm  to  ~1.10  pm  [10].  Some  of  the  highest 
power  cw  fiber  lasers  to  date  have  been  operated  using  ytterbium  lasing  at  1.12  pm  [66]. 
This  pushes  the  Raman  gain  spectrum  further  into  the  infrared  where  it  overlaps  the 
holmium  absorption  curve  more  effectively.  The  Raman  suppression  will  be  most 
efficient  near  1.094  pm,  forcing  the  Stokes  gain  peak  to  overlap  the  holmium  absorption 
peak  at  1.15  microns.  The  Raman  gain  curve  and  Holmium  absorption  are  plotted  in 
Figure  65(b). 
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;  65:  Holmium  absorption  in  silicate  glass  (dotted  line)  along  with  the  Raman  gain 
(solid  line)  as  a  function  of  the  frequency  shift  from  a  neodymium  pump  beam  at 
1 .064  microns  (a)  and  as  a  function  of  the  frequency  shift  from  a  ytterbium  pump 
beam  at  1 .094  microns  (b). 


7.2.2.  Threshold  Calculation 

Smith  defines  a  convenient  form  of  the  SRS  threshold.  In  the  limit  of  no  pump 
depletion,  the  input  intensity  that  yields  an  equivalent  Stokes  beam  at  the  exit  facet  is 
identified  as  the  threshold  [18].  This  definition  is  relatively  simple  to  implement.  Since 
the  assumption  is  made  that  the  pump  power  in  the  fiber  is  only  affected  by  absorption, 
the  differential  equations  coupling  the  pump  and  Stokes  beams  can  be  solved 
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Ho  absorption(dB/(m-mol%))  Ho  absorption(dB/(m-mol%)) 


analytically.  Equation  (7.3)  gives  the  intensity  of  one  speetral  eomponent  of  the  Stokes 
beam  at  the  exit  faeet  of  the  fiber  for  a  given  pump  intensity.  Assuming  the  pump  beam 
is  monoehromatie,  the  only  terms  on  the  right  hand  side  with  any  speetral  dependenee  are 
the  Raman  gain  and  the  Stokes  absorption,  resulting  in  the  total  Stokes  power  of 


-  exp  (-«/))- (k) 


L 


dv 


(7.5) 


where  is  the  effeetive  area  of  the  fiber  and  the  integration  is  performed  over  the 

Raman  gain  bandwidth.  is  the  power  of  the  noise  Stokes  photons  that  give  rise  to  the 

stimulated  seattering.  It  is  approximated  by  the  intensity  of  a  single  photon  per  seeond 
transmitted  through  the  fiber.  The  integration  is  performed  over  all  frequeneies  for  whieh 
the  Raman  gain  is  greater  than  zero.  The  interpretation  of  the  remaining  terms  has  not 
ehanged  from  their  earlier  uses. 

Sinee  the  speetral  dependenee  of  the  Raman  gain  is  not  deseribed  by  an  analytieal 
funetion,  Equation  (7.5)  is  best  evaluated  using  a  numerieal  integration. 


^^^-fci^(l-exp(-«^E))-«,  (v) 


L 


(7.6) 


where  are  the  frequeney  eomponents  for  whieh  the  Raman  gain  is  tabulated  and 
A  u  =  -  K;  is  the  resolution  of  the  tabulation.  Stimulated  seattering  threshold  is  found 

by  evaluating  Equation  (7.6)  for  inereasing  pump  intensities,  ,  until  the  output  Stokes 
power  is  equal  to  the  input  pump  power. 
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This  analysis  was  completed  at  various  holmium  dopant  levels  for  a  fiber  pumped  at 
1.064  pm  and  at  1.09  pm.  A  100  m  length  of  Corning  SMF-28  fiber  with  an  8-micron 
mode  field  diameter  was  modeled  for  the  calculation.  The  undoped  threshold  was 
100  watts.  Adding  holmium  to  the  fiber  inereased  the  threshold  almost  linearly  resulting 
in  Figure  66. 


66:  The  threshold  for  SRS  generation  pumped  at  1.064  mierons  (dashed  line)  and 
at  1.09  microns  (solid  line)  calculated  versus  the  molar  %  of  holmium  doped  into 

the  fiber. 


7.3,  Conclusions 

Figure  66  shows  that  doping  a  fiber  with  holmium  can  significantly  affect  SRS 
threshold.  A  factor  of  two  increase  in  the  threshold  was  achieved  with  a  pump  beam  at 
1.064  (1.094)  pm  using  a  dopant  level  of  0.053  (0.012)  molar  percent  of  holmium. 
However,  the  effeets  of  the  dopant  on  the  pump  beam  have  not  been  tested.  The  poor 
optical  quality  of  the  holmium  doped  silicate  windows  prevented  a  quantitative  measure 
of  absorption  at  the  pump  wavelengths.  Although  absorption  is  4.95  dB/(m  mol%) 
greater  at  1.09  pm  than  at  1.064  pm,  the  greater  absorption  at  the  Stokes  peak 
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compensates  for  the  extra  loss  at  1 .09  mierons.  Absorption  at  the  Raman  gain  peak  is 
79.4  dB/(m  mol%)  greater  when  the  laser  operates  at  1.094  pm. 

The  slope  of  the  curves  in  Figure  66  is  determined  by  the  ratio  of  the  SRS  gain  to  the 
holmium  absorption  at  the  dominant  Stokes  wavelength.  At  low  dopant  levels,  most  of 
the  Stokes  power  is  contained  in  frequencies  near  the  peak  of  the  Stokes  wavelength,  but 
at  high  dopant  levels,  the  peak  of  the  Raman  gain  is  suppressed  and  the  Stokes  beam  is 
generated  at  frequencies  near  the  edge  of  the  gain  spectrum.  The  kink  in  the  slope  is  due 
to  a  discrete  shift  in  the  dominant  SRS  peak. 

Several  factors  that  affect  integration  of  this  Raman  suppression  technique  in  a  fiber 
laser  were  ignored  in  the  analysis  presented  here.  First,  the  thresholds  reported  in  Figure 
66  do  not  represent  the  maximum  power  achievable  in  a  fiber  laser  before  the  effects  of 
SRS  are  visible.  Depending  on  the  reflectivity  of  the  output  coupler  and  the  fiber  length, 
the  oscillating  power  in  the  cavity  could  vary  considerably  along  the  fiber  length  and 
eould  be  significantly  stronger  than  the  power  exiting  the  fiber;  hence  there  is  a  strong 
cavity  dependence  that  will  affect  the  power  level  that  can  be  achieved  before  the  onset 
of  SRS. 

Secondly,  although  the  attenuation  at  the  pump  and  laser  frequencies  was  eonsidered, 
the  addition  of  a  dopant  to  the  laser  cavity  can  have  other  unforeseen  effects.  Ytterbium 
has  been  used  to  sensitize  holmium  to  InGaAs  laser  diodes  near  940  nm.  This  transition 
in  ytterbium  is  broad  and  once  excited  to  that  level,  it  easily  transfers  energy  to  the 
holmium  ion.  This  energy  transfer  process  is  used  to  increase  holmium  luminescence  in 
the  green  [76-78]  and  to  generate  lasers  operating  at  2.1  microns  [79,  80].  The  rate  of 
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energy  transfer  from  the  ytterbium  to  the  holmium  depends  on  the  eoneentration  of  both 
ions.  Aeeording  to  Zhang  et  al.  [77,  78],  the  rate  inereases  linearly  with  holmium 
eoneentration  and  quadratieally  with  ytterbium  eoneentration,  so  even  if  the  pump 
absorption  does  not  limit  the  holmium  dopant  level,  the  energy  transfer  between  the 
holmium  and  ytterbium  laser  will. 

Finally,  the  laser  dopants  as  well  as  the  holmium  will  affeet  the  Raman  gain 
speetrum.  For  low  doping  levels,  the  silieate  glass  will  dominate  the  Raman  speetrum. 
However,  at  high  dopant  levels,  new  resonant  Raman  peaks  may  appear  that  are  not 
neeessarily  suppressed  by  holmium.  These  peaks  would  require  eonsideration  in  a  laser 
system. 

Clearly  there  are  many  issues  that  must  be  addressed  before  SRS  suppression  with  a 
holmium  intereavity  dopant  is  feasible  in  a  fiber  laser  system.  However,  some  of  the 
theoretieal  groundwork  has  been  laid  showing  a  remarkable  improvement  in  SRS 
threshold.  At  0.3  molar  %  doping  level,  a  100  m  long  ytterbium  laser  is  predieted  to  be 
able  to  operate  at  19  times  the  power  level  that  an  undoped  fiber  eould  support. 
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8.  SRS  Reduction  Using  Long  Period  Gratings 


Although  introducing  a  dopant  into  a  fiber  laser  was  shown  in  Chapter  7  to 
theoretically  increase  the  Stokes  threshold,  it  is  not  a  flexible  solution.  Figure  65  showed 
that  a  neodymium  laser  operating  at  1 .064  mierons  has  a  much  lower  Stokes  absorption 
than  an  ytterbium  laser  operating  at  1.094  microns.  Additionally,  operation  of  the 
ytterbium  laser  at  longer  wavelengths  would  result  in  severe  absorption  of  the  laser  beam. 
A  long  period  grating  on  the  other  hand  can  be  tailored  to  eliminate  the  unwanted  Stokes 
beam  while  minimally  affecting  the  pump  beam. 

Long  period  gratings,  LPGs,  are  similar  to  fiber  Bragg  gratings  in  that  they  are 
periodic  variations  in  the  index  of  refraction  that  affect  the  propagation  of  light  in  the 
fiber.  The  results  of  the  structures  are  significantly  different  though.  Fiber  Bragg 
gratings  reflect  light  over  a  very  narrow  bandwidth  (<  1  nm)  while  long  period  gratings 
couple  light  from  a  speeific  eore  guided  mode  to  a  cladding  guided  mode  and  have  much 
larger  bandwidths.  In  most  fibers,  cladding  guided  modes  are  very  lossy  thus  the  long 
period  grading  essentially  beeomes  a  frequency  dependent  loss  in  the  fiber  (Figure  67). 


Figure  67:  A  long  period  grating  couples  light  out  of  the  core  and  into  the  fiber  cladding 
where  high  losses  account  for  the  reduction  in  transmission. 
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This  method  for  increasing  the  loss  at  the  Stokes  wavelength  is  investigated 
theoretically  in  this  chapter.  The  spectral  characterization  of  the  Raman  gain  is  examined 
in  Section  8.1.  In  Section  8.2  the  design  of  the  long  period  grating  is  explained,  and 
conclusions  are  made  in  Section  8.3. 

8,1.  Raman  Spectrum  Measurement 

An  intense  ultraviolet  light  source  can  generate  a  permanent  change  in  the  refractive 
index  of  a  photosensitive  fiber.  This  effect  is  responsible  for  long  period  gratings  and  can 
be  improved  with  the  addition  of  certain  dopants  to  the  fiber,  mainly  germanium.  To 
maximize  the  effect  of  the  LPG  on  the  SRS  threshold,  a  highly  photosensitive  fiber  from 
StockerYale,  PS-SMF-30,  was  purchased.  Jain’s  group  at  the  University  of  New  Mexico 
measured  the  index  change  in  this  fiber  as  well  as  in  a  hydrogen  loaded  standard  single¬ 
mode  fiber  and  found  that  it  was  a  factor  of  170  greater  in  the  photosensitive  fiber.  The 
additional  dopants,  however,  can  affect  the  Raman  gain  spectrum,  forcing  a  measurement 
of  this  spectrum  before  a  final  grating  design  is  possible. 

8.1.1.  Experimental  Setup 

The  Raman  spectrum  was  measured  using  the  experimental  setup  shown  in  Figure  68. 
A  frequency  doubled  Q-switched  Continuum  Powerlite  NdiYAG  laser  was  used  as  the 
laser  source.  Delaying  the  Q-switch  relative  to  the  fiashlamp  and  transmitting  the  beam 
through  a  half-wave  plate  and  a  polarizer  controlled  the  laser  power.  The  beam  was  then 
coupled  into  the  photosensitive  fiber.  At  the  output  end  of  the  fiber,  a  microscope 
objective  collimated  the  beam,  which  then  propagated  through  a  Raman  notch  filter  that 
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reduced  the  pump  power.  The  Stokes  beam  was  then  coupled  into  a  fiber  bundle  attached 
to  a  spectrum  analyzer  and  an  optical  multichannel  analyzer,  with  a  resolution  of  14  cm'*. 


Figure  68:  Schematic  for  Raman  spectrum  measurement. 


The  experiment  was  performed  at  various  fiber  lengths  between  0.5  and  200  m.  The 
pump  power  was  also  varied.  However,  it  could  not  be  controlled  well  because  the  laser 
was  operated  with  the  Q-switch  delay  not  optimized  to  reduce  laser  power.  To  minimize 
the  random  variations  from  one  measurement  to  the  next,  100  spectral  scans  were  taken 
at  each  power  level  and  averaged.  The  measurement  was  further  complicated  by  the  high 
power  of  the  laser.  A  small  percentage  of  the  beam  could  easily  damage  the  front  fiber 
facet. 
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8.1.2.  Results 


A  complicated  Stokes  spectrum  is  revealed  in  Figure  69(a).  This  data  was  collected 
using  a  0.5  m  long  fiber.  At  low  power  levels,  peaks  in  the  Stokes  spectrum  are  found  at 
265  cm'^  and  480  cm'\  As  the  power  increases,  these  peaks  shift  to  longer  wavelengths, 
290  cm'^  and  540  cm’'  respectively,  and  another  peak  at  440  cm’'  appears  in  the 
spectrum.  Figure  69(b),  shows  the  Stokes  spectrum  after  transmission  through  a  slightly 
longer  fiber,  1  m.  This  inerease  in  interaction  length  reduces  the  threshold  for  SRS.  The 
results  show  that  as  the  440  cm'^  line  grows,  the  290  and  540  cm"'  peaks  remain 
relatively  constant.  Even  longer  fibers  have  confirmed  this  trend.  The  Stokes  spectrum 
of  a  200m  fiber,  shown  in  Figure  69(c),  demonstrates  that  the  440cm‘'  peak  dominates 
the  others. 
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Figure  69:  Stokes  spectra  at  different  pump  power  levels  in  (a)  a  0.5  m  fiber,  (b)  a  1  m 
fiber,  and  (c)  a  200  m  fiber.  The  dashed  lines  in  (a)  and  (b)  indicate  the 
transmission  spectrum  of  the  Raman  notch  filter  used  in  those  experimental  setups 

to  attenuate  the  pump  beam. 
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The  growth  of  a  frequency  component  at  440  cm'^  is  expected  based  upon  the  Raman 
gain  spectrum  of  both  SiOi  and  fused  GeOi  (see  Figure  70).  Conversely,  the  Stokes 
peaks  at  290  and  540  cm'^  and  their  variation  with  pump  power  is  unexplained  by  the 
Raman  gain  spectrum  alone.  Qualitatively,  however,  sidelobes  of  the  pump  beam 
generated  through  modulation  instability  (MI)  have  many  of  the  same  characteristics  as 
these  two  Stokes  peaks.  Modulation  instability  is  a  four-wave  mixing  phenomenon  that 
is  phase-matched  through  self-  or  cross-phase  modulation  [3].  In  single  mode  fibers,  MI 
is  typically  observed  in  the  anomalous  group-velocity  dispersion  region  above 
-1.3  microns.  However,  in  multimode  fibers  the  pump,  Stokes  and  anti-Stokes 
components  can  travel  in  different  fiber  modes.  The  difference  in  the  effective  indices  of 
these  modes  can  offset  the  material  and  nonlinear  dispersion  [81,  82].  Because  the 
nonlinear  dispersion  is  intensity  dependent,  the  frequencies  for  which  the  beams  are 
phase  matched  changes  with  pump  power,  thus  explaining  the  variation  in  the  Stokes 
spectrum  at  different  power  levels. 


Figure  70:  Raman  gain  spectrum  for  various  types  and  concentrations  of  dopants  [83]. 
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Because  of  these  complications  with  the  measured  spectrum  and  due  to  the  theoretical 
nature  of  the  research  described  here,  the  measured  Raman  gain  spectrum  was  not  used  in 
the  analysis  described  below.  The  PS-SMF-30  fiber  is  known  to  have  a  high 
concentration  of  germanium  thus  for  the  purpose  of  theoretically  demonstrating  the  affect 
of  a  long  period  grating  on  stimulated  Raman  scattering,  the  gain  spectrum  of  Si02  doped 
with  25  molar  %  of  Ge02  is  used. 

8,2.  Long  Period  Grating  Design 

8.2.1.  Theory 

The  LPG  design  hinges  upon  a  correct  calculation  of  the  transmission  characteristics 
of  a  long  period  grating.  The  steps  laid  out  by  Erdogan  [12,  16]  were  used  to  determine 
the  LPG  transmission.  These  calculations  are  clearly  explained  and  quite  intricate.  For 
these  reasons,  only  the  critical  highlights  are  mentioned  here.  The  results  of  these 
calculations  will  be  used  to  design  an  appropriate  LPG  to  reduce  SRS  and  theoretically 
demonstrate  its  expected  effect  on  the  SRS  threshold. 

As  mentioned  previously,  a  long  period  grating  couples  a  core-guided  fiber  mode  to  a 
cladding  guided  fiber  mode.  In  essence,  it  is  a  volume  hologram  in  the  fiber  core  with 
the  core-guided  mode  as  the  reference  wave  and  the  cladding-guided  mode  as  the  object 
wave  [84].  Thus  the  grating  should  have  an  index  variation  that  precisely  matches  the 
interference  pattern  of  the  two  modes.  Such  an  interpretation  immediately  demonstrates 
that  the  LPG  should  have  a  period  equal  to  the  beat  frequency  between  the  core-  and 
cladding-guided  modes.  Since  both  of  these  waves  travel  in  the  same  direction,  the 
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grating  has  a  long  period  relative  to  that  of  a  Bragg  grating  in  whieh  the  reference  wave 
propagates  in  the  opposite  direction  as  the  object  wave. 

Depending  on  the  cladding  mode  that  the  beam  is  coupled  into,  the  period  of  the 
grating  will  change.  Figure  4,  page  17  of  Section  2.2,  shows  the  cladding  mode 
dispersion  relation.  The  independent  variable  in  this  figure  is  the  effective  index  and  it  is 
apparent  from  this  graph  that  over  the  region  shown,  the  difference  between  adjacent  odd 
modes  increases  with  mode  number.  Thus  the  resonance  peaks  for  higher  order  modes 
will  be  much  further  apart  than  those  for  lower  order  modes.  This  is  most  significant 
when  wavelength  regions  where  high  transmission  is  required  surround  the  desired 
resonance.  For  the  case  at  hand,  the  Stokes  light  must  be  eliminated  without  inducing 
significant  loss  at  the  pump  wavelength.  Thus  using  certain  cladding  modes  may  be 
problematic  if  one  of  the  lower  order  sidelobes  overlaps  the  pump  beam.  On  the  other 
hand,  apodization  of  the  LPG  can  reduce  unwanted  sidelobes  [84],  but  would  complicate 
the  analysis  and  is  not  attempted  here. 

The  LPoi  fiber  mode  will  efficiently  couple  to  only  the  HEin  fiber  modes,  where  n  is 
an  odd  number.  All  of  the  HEin  cladding  modes  are  circularly  symmetric,  a  characteristic 
that  is  required  to  produce  an  interference  pattern  with  the  EPoi  core  mode  that 
approximates  a  fiat  grating.  Additionally,  as  can  be  seen  by  Eigure  5,  page  18  of  Section 
2.2,  the  odd  modes  have  local  maxima  at  the  fiber  core  and  overlap  the  fundamental  core- 
guided  mode  much  better  than  the  even  modes.  As  the  mode  number  increases,  this 
generalization  loses  its  validity,  but  it  is  still  important.  The  overlap  of  the  core  and 
cladding  modes  is  closely  tied  to  the  coupling  constant,  which  is  a  proportional  to  the 
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efficiency  with  which  the  core  mode  is  coupled  to  the  cladding  mode.  Figure  71  shows 
the  coupling  constant  from  the  core  mode  to  various  cladding  modes.  As  this  plot  shows, 
up  to  relatively  high  cladding  mode  numbers,  the  beam  is  coupled  to  the  odd  modes  much 
stronger  than  to  the  even  modes.  Additionally,  it  becomes  apparent  that  the  peak  in  the 
coupling  constant  of  the  odd  modes  is  not  when  the  LPqi  core  mode  is  coupled  into  the 
HE  11  cladding  mode,  but  rather  into  a  higher  odd  mode.  Although  using  this  mode  would 
allow  for  a  grating  with  a  wider  spectral  bandwidth,  complications  result  from  lower 
order  modes  overlapping  the  pump  beam. 
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Figure  71:  Proportionality  constant  for  the  efficiency  of  the  coupling  between  the  LPqi 
and  cladding  modes  [12].  The  odd  cladding  modes  are  designated  with  triangles. 
The  even  cladding  modes  are  designated  with  squares. 


The  transmission  characteristics  of  a  long  period  grating  are  sensitive  to  the  depth  of 
the  index  grating  written  into  the  fiber.  This  can  be  problematic  since  the  central 
wavelength  must  overlap  the  peak  of  the  Raman  gain  curve  to  strongly  affect  the  SRS 
threshold.  However,  after  the  grating  is  written  into  the  fiber,  the  transmission  can  be 
adjusted  slightly.  It  has  been  shown  that  by  etching  the  cladding  slightly,  thus  slightly 
changing  its  radius,  the  resonance  wavelength  will  increase  [85]. 
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Another  complication  with  the  design,  as  implemented  here,  is  that  the  analysis 
neglects  the  fact  that  the  PS-SMF-30  fiber  supports  many  core-guided  modes  at  1.064 
microns.  This  will  clearly  change  the  results,  since  the  LPG  is  designed  only  to  induce 
coupling  between  the  LPoi  mode  and  the  cladding  modes.  There  will  be  accidental 
resonances  between  higher  order  core-guided  modes  and  the  cladding.  However,  the 
position  of  these  resonances  has  not  been  calculated  and  it  is  unlikely  that  any  of  these 
will  be  useful  in  the  suppression  of  SRS.  Under  most  situations,  however,  the  fiber  laser 
or  amplifier  in  question  will  only  support  a  single  fiber  mode  to  ensure  a  high  quality 
laser  beam.  Additionally,  as  demonstrated  in  Section  6.2,  even  in  multimode  fibers,  SRS 
was  generated  predominantly  in  the  LPqi  mode.  Thus  the  application  of  a  long  period 
grating  to  this  problem  is  appropriate. 

Once  the  transmission  of  the  long  period  grating  and  the  Raman  gain  spectrum  is 
determined,  they  can  be  used  to  calculate  the  threshold  for  SRS  in  a  way  that  is  similar  to 
that  outlined  in  the  previous  chapter.  Differences  between  the  two  cases  arise  from  the 
fact  that  rather  than  absorption  along  the  fiber,  the  loss  is  localized  to  discrete  points.  To 
account  for  this,  each  section  of  fiber  can  be  treated  independently.  The  power  in  the 
Stokes  beam  will  initially  grow  unimpeded,  seeded  by  the  spontaneous  scattering  in  the 
fiber.  The  beam  will  then  pass  through  the  long  period  grating  and  power  will  be  lost. 
This  reduced  Stokes  power  is  then  the  seed  for  the  next  stretch  of  unimpeded  Stokes 
growth.  This  pattern  continues  until  the  beam  exits  the  fiber  and  from  Equation  (7.6)  the 
total  Stokes  power  is  given  by 
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(8.1) 


gr[yi)lp<> 


a 


(l- exp  (-«£))-«£ 


In  this  equation  the  transmission  of  the  grating  at  a  frequeney  v  is  given  by  T{v^  ,  n 

is  the  number  of  gratings  in  the  fiber,  and  the  absorption  at  the  pump  and  Stokes 
wavelengths  is  assumed  to  be  eonstant,  «  =  •  The  summation  over  i  is 

eompleted  to  obtain  a  numerieal  integration  over  the  frequeneies  for  whieh  there  is  gain. 
Finally,  the  total  Stokes  power  is  eompared  to  the  pump  power.  In  the  same  manner  as 
deseribed  by  Smith  [18]  and  in  Chapter  7,  threshold  is  identified  as  the  pump  power  that 
produees  an  equally  powerful  Stokes  beam. 

The  optimal  grating  design  maximizes  the  SRS  threshold  while  maintaining  a  high 
transmission  at  the  pump  wavelength.  For  the  design  performed  here,  the  produet  of  the 
Raman  threshold  and  the  transmission  at  the  pump  wavelength  is  maximized  relative  to 
the  grating  period  and  the  grating  length.  The  grating  period  determines  the  eenter 
wavelength  and  mode  order  as  deseribed  above.  The  grating  length  is  also  eritioal. 

Unlike  a  distributed  Bragg  refleetor,  a  longer  grating  does  not  neeessarily  result  in  a 
stronger  eoupling.  At  some  point,  power  will  begin  to  flow  baek  from  the  cladding  mode 
into  the  core  mode.  Thus  the  LPG  should  not  be  greater  than  this  length.  For  the 
purposes  described  here,  the  grating  is  slightly  shorter  than  the  length  resulting  in  best 
coupling.  A  shorter  length  increases  the  spectral  bandwidth  for  a  better  overlap  with  the 
Raman  gain  curve. 
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8.2.2.  Results 


Using  a  single  grating  to  maximize  the  figure  of  merit  described  above  results  in  the 
transmission  curve  shown  in  Figure  72(a).  The  transmission  at  1.064  microns  is  nearly 
98%.  The  Raman  gain  curve  is  overlaid  to  demonstrate  the  position  of  the  gain  peak. 

The  grating  characteristics  are  summarized  in  Table  8.  The  SRS  threshold  was  calculated 
to  be  42. IW.  This  is  a  factor  of  19%  greater  than  the  calculated  threshold  when  no 
gratings  are  present  in  the  fiber,  35.4  W.  The  low  threshold  relative  to  that  calculated  in 
the  previous  chapter  with  no  dopant  in  the  fiber  is  a  direct  result  of  the  increased  gain  of  a 
germanium  doped  fiber  (see  Figure  70). 


Frequency  (Ay)  Frequency  (Ay) 

(c)  (d) 

:  72:  LPG  transmission  (solid)  and  Raman  gain  (dashed)  as  a  function  of  the 
~  "equency  difference  from  the  pump  beam  at  1 .064  microns,  (a)  A  single  long 
period  grating  with  the  HE^  mode  centered  at  440  cm  ',  (b)  A  single  grating 
design  placed  at  three  points  in  the  fiber,  (c)  One  set  of  three  different  gratings,  (d) 
Three  sets  of  three  different  gratings. 
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Table  8:  Critical  parameters  for  several  LPG  designs. 


Gratings  per  set 

0 

1 

1 

3 

3 

#  of  sets 

0 

1 

3 

1 

3 

Center  Grating  Period  (qm)  -  A 

NA 

236.6 

236.7 

236.7 

236.9 

Grating  Length  (cm)  -  L 

NA 

1.12 

1.09 

1.10 

1.08 

AA  Between  Gratings  (qm) 

NA 

NA 

NA 

0.60 

1.07 

Threshold  Power  (W) 

35.4 

42.0 

51.0 

52.1 

72.3 

Threshold  Increase  (%) 

0 

19 

44 

47 

104 

Tx  @  1.064  qm 

100 

97.7 

95.6 

95.0 

84.4 

Simply  adding  more  of  these  gratings  will  further  increase  the  SRS  threshold. 
However,  for  optimal  performance,  a  slightly  different  design  should  be  implemented 
when  multiple  gratings  are  used.  Optimizing  the  same  figure  of  merit  requires  a  grating 
with  a  slightly  shorter  length  and  a  shifted  center  frequency.  The  shorter  length  widens 
the  spectral  width  of  the  grating,  which  suppresses  SRS  over  the  full  gain  curve  better. 
The  depth  of  the  overall  transmission  dip  is  maintained  by  adding  the  extra  gratings  to  the 
fiber.  The  result  of  optimizing  the  design  for  three  gratings  in  the  fiber  is  shown  in 
Figure  72(b)  and  again  the  critical  characteristics  are  summarized  in  Table  8. 

Although  it  is  a  simple  matter  to  just  add  more  of  a  single  grating  to  a  fiber,  the 
Raman  gain  curve  could  be  better  matched  by  a  set  of  LPG’s,  each  with  a  different  center 
frequency.  To  test  this  idea,  a  single  set  of  three  gratings  were  designed  that  would 
maximize  threshold.  To  simplify  the  optimization,  all  gratings  were  assumed  to  be  the 
same  length,  the  grating  periods  were  constrained  to  be  equally  spaced  (A  -  AA ,  A ,  and 
A  +  AA ),  and  there  is  assumed  to  be  no  interaction  between  the  gratings.  The  resulting 
transmission  curve  is  shown  in  Figure  72(c)  and  the  grating  characteristics  are 
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summarized  in  Table  8.  This  design  inereased  SRS  threshold  by  a  factor  of  47%  to 
52.1  W,  while  maintaining  a  relatively  high  transmission  at  1.064  microns  of  95%. 

For  completeness,  another  case  was  investigated  with  three  sets  of  gratings  in  the 
fiber.  Each  set  had  three  gratings  with  different  periods  for  a  total  of  nine  gratings 
written  into  the  fiber.  The  transmission  curve  and  grating  characteristics  are  shown  in 
Figure  72(d)  and  Table  8  respectively.  This  design  resulted  in  an  increase  of  104%  in 
SRS  threshold  to  72.3  W. 

8,3.  Conclusions 

Based  upon  the  theoretical  results  described  above,  a  set  of  long  period  gratings  could 
clearly  increase  SRS  threshold  significantly.  However,  there  were  many  simplifications 
made  to  make  the  problem  tractable.  The  expected  effects  of  these  simplifications  are 
described  below. 

First,  the  pump  beam  was  only  assumed  to  pass  through  the  fiber  one  time.  Although 
this  appropriately  models  a  fiber  amplifier,  this  is  obviously  not  the  case  for  a  fiber  laser 
where  feedback  is  required.  A  laser  may  be  affected  by  the  grating  in  a  remarkably 
different  way.  Depending  on  the  cavity,  the  laser  beam  would  bounce  many  times  back 
and  forth  through  the  cavity,  thus  the  transmission  of  the  pump  beam  would  become 
extremely  important.  It  would  be  appropriate  to  increase  the  importance  of  its 
transmission  in  the  figure  of  merit  used  to  optimize  the  grating.  Additionally,  if  there  is 
any  feedback  of  the  Stokes  beam  by  the  cavity  mirrors,  SRS  threshold  would  be 
remarkably  affected. 
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It  was  also  assumed  that  the  index  variation  was  perfeetly  sinusoidal  over  the  grating 
length.  This  simplified  the  analysis,  but  any  variations  in  the  grating  period  or  index 
depth  will  smear  the  transmission  eurve.  Additionally,  an  intentional  apodization  of  the 
index  variation  eould  eliminate  eoupling  between  the  pump  and  lower  order  eladding 
modes  [84],  thus  inerease  the  transmission  at  1.064  mierons.  This  would  also  allow  for 
efficient  use  of  higher  order  cladding  modes  that  have  a  stronger  overlap  with  the  LPqi 
core -guided  mode. 

Along  the  same  lines,  a  transmission  band  can  be  generated  within  a  long  period 
grating  using  a  n  phase  shift  in  the  center  of  the  grating  [86].  Thus  by  centering  the 
transmission  band  at  the  pump  wavelength,  a  high-order  cladding  mode  with  a  broad 
spectral  bandwidth  could  be  used.  Again,  this  could  better  match  the  Raman  gain  curve 
while  maintaining  the  required  high  pump  transmission. 

Finally,  the  analysis  performed  here  neglects  the  interaction  between  adjacent 
gratings  and  coupling  from  the  core  into  multiple  cladding  modes.  Since  the  calculation 
assumed  that  only  the  core-guided  mode  and  a  single  cladding  guided  mode  were 
coupled,  an  overlapping  third  mode  would  clearly  change  the  results  [12].  This  is 
expected  to  affect  the  transmission  near  200  cm'^  where  coupling  into  the  HEn  and  the 
HEi3  modes  nearly  overlap,  but  there  should  not  be  a  significant  problem  at  the  Raman 
gain  peak  unless  successive  gratings  in  a  set  are  closely  spaced.  This  problem  would 
have  to  be  characterized  to  fully  appreciate  its  effect  on  the  transmission  curve. 
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9.  Conclusions 


9.1,  Technical  achievements 

The  goal  of  the  researeh  examined  in  this  doeument  was  to  generate  high  intensity 
laser  beams  in  optieal  fibers.  This  problem  was  investigated  from  a  variety  of  angles.  In 
Chapter  4,  the  spatial  eharaeteristics  of  stimulated  Brillouin  seattering  in  a  long 
multimode  fiber  were  examined.  It  was  clearly  demonstrated  in  this  chapter  that  long, 
multimode  fibers  generated  a  Gaussian-like  Stokes  beam  even  when  the  pump  beam  was 
highly  aberrated.  This  feature  of  SBS  was  also  examined  theoretically.  It  was  found  that 
for  many  pump  conditions,  the  finite  Stokes  shift  due  to  SBS  did  not  allow  for  phase 
conjugation  in  long  fibers  and  resulted  in  a  stronger  gain  for  the  LPqi  fiber  mode.  On  the 
other  hand,  the  polarization  states  of  the  beam  were  preserved  by  Brillouin  reflection.  An 
elliptically  polarized  pump  beam  was  found  to  reflect  into  an  elliptically  polarized  Stokes 
beam  with  the  same  handedness,  identical  to  the  reflection  properties  of  a  phase 
conjugate  mirror.  These  results  were  found  to  agree  with  a  theoretical  model  that 
assumes  orthogonal  polarization  states  reflect  off  the  same  density  grating  in  the  fiber. 
Additionally,  this  chapter  examined  the  characteristics  of  higher  order  modes  excited  in 
the  fiber.  It  was  found  that  under  specific  pumping  conditions,  an  LPn  fiber  mode  could 
be  excited.  However,  small  perturbations  to  the  gain,  such  as  launching  a  second  pump 
beam  at  the  same  frequency  into  the  fiber,  would  result  in  the  excitation  of  the  LPoi 
mode.  Finally,  the  chapter  was  concluded  with  an  analysis  of  SBS  beam  cleanup  in  a 
ring  cavity  design.  When  there  was  strong  coupling  of  the  Stokes  beam  back  into  the 
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fiber,  the  beam  quality  suffered.  On  the  other  hand,  a  small  seed  was  found  to  inerease 
the  Stokes  output  with  a  minimal  affeet  on  the  beam  quality. 

Beam  eombining  using  stimulated  Brillouin  seattering  was  examined  in  Chapter  5. 
Sinee  the  fundamental  fiber  mode  was  exeited  regardless  of  the  pumping  eonditions,  two 
pump  beams  eoupled  into  the  fiber  at  different  positions  eould  be  spatially  eombined. 

The  effeet  was  quantitatively  analyzed  using  lasers  at  the  same  frequency  that  worked 
together  to  produce  two  temporally  coherent  Stokes  beams  as  well  as  using  two  lasers  at 
different  frequencies  that  produced  two  temporally  incoherent  Stokes  beams.  In  both 
cases,  however,  the  Stokes  beams  were  found  to  have  nearly  perfect  spatial  overlap. 
Coherent  beam  combining  has  the  benefit  of  increased  overall  efficiency  over  that  of 
incoherent  beam  combining.  However,  the  frequency  requirements  on  the  two  pump 
beams  are  very  stringent,  thus  incoherent  combining  is  a  more  practical  solution  to 
combining  multiple  lasers  in  many  situations.  During  the  course  of  this  work,  the 
efficiency  was  measured  as  a  function  of  the  coupling  position  and  found  to  fall  off  as  the 
beam  was  coupled  into  the  fiber  edges.  Additionally,  the  polarization  state  of  the  Stokes 
beam  was  found  to  mimic  that  of  the  pump  in  the  same  manner  as  it  did  for  single  beam 
cleanup.  Finally,  the  spatial  coherence  of  the  beams  was  demonstrated  by  the  linear 
interference  pattern  generated  when  different  portions  of  the  Stokes  beam  overlapped  on 
an  observation  screen. 

In  Chapter  6,  the  limitations  to  stimulated  Brillouin  scattering  beam  cleanup  and 
combining  were  discussed  along  with  several  methods  that  could  be  used  to  overcome 
these  limitations.  Section  6.1  examined  the  requirements  on  an  optical  system  designed 
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to  combine  more  than  two  beams  and  effieiently  eollect  the  Stokes  output  into  a  usable 
beam.  The  ideas  presented  here  showed  that  eombining  four  sources  would  clearly  be 
possible.  It  also  demonstrated  a  path  to  incorporate  more  pump  beams  into  the  design. 

The  following  section  investigated  the  limitation  that  second  order  SBS  would  have 
on  a  beam  combining  or  coherent  beam  cleanup  system.  Using  a  relatively  small  fiber, 
the  Stokes  output  was  monitored  to  find  the  seeond  order  threshold.  This  was  then 
modeled,  taking  into  aecount  the  forward  propagating  pump,  backward  propagating 
Stokes,  the  portion  of  the  Stokes  beam  that  reflects  off  the  front  fiber  facet  and 
propagates  forward  in  the  fiber,  and  the  seeond  order  Stokes  beam.  These  beams  were 
coupled  together  through  stimulated  Brillouin  scattering  and  four-wave  mixing.  Using 
this  model,  it  was  found  that  without  four-wave  mixing  to  seed  second  order  SBS,  the 
seeond  order  threshold  would  be  much  higher  and  would  not  be  a  significant  limiting 
factor.  This  was  verified  to  some  extent  experimentally  by  angle  cleaving  the  front  fiber 
facet,  thus  eliminating  the  forward  propagating  Stokes  beam.  In  this  experiment  no 
second  order  SBS  was  observed  at  power  levels  much  higher  than  the  original  seeond 
order  threshold. 

Finally,  in  Section  6.2,  SRS  was  investigated  as  an  alternative  to  SBS  for  beam 
eleanup  applieations.  SRS  has  the  benefit  of  an  increased  gain  bandwidth,  thus  relaxing 
the  spectral  requirements  on  the  pump.  In  this  section,  the  beam  quality  of  the  Stokes 
power  was  investigated  and  found  to  be  nearly  nine  times  better  than  that  of  the 
transmitted  pump  beam.  This  has  the  potential  of  increasing  the  power  on  target  by  a 
factor  of  74. 
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The  remaining  ehapters  shift  from  using  the  benefits  of  stimulated  scattering  to 
eliminating  stimulated  Raman  scattering  in  optical  systems  for  which  SRS  can  be 
detrimental.  In  Chapter  7,  the  effect  of  an  intercavity  dopant  is  modeled  in  a  liber.  The 
dopant,  holmium,  is  found  to  have  a  strong  absorption  feature  at  the  Raman  gain  peak  for 
lasers  centered  at  1.094  microns.  This  wavelength  is  well  within  the  gain  curve  of 
ytterbium  doped  fibers  and  SRS  threshold  could  be  increases  by  a  factor  of  two  with  a 
relatively  small  amount  of  dopant,  0.012  molar  percent.  Even  when  pumped  by  a 
neodymium  laser  at  1.064  microns,  the  absorption  of  holmium  is  found  to  have  a 
significant  effect  on  the  SRS  threshold.  The  SRS  threshold  could  be  doubled  with  0.053 
molar  percent  of  holmium.  The  effect  of  this  dopant  on  the  transmission  at  the  pump 
wavelength  was  only  qualitatively  analyzed.  However,  it  was  expected  to  be  relatively 
low. 

Finally,  in  Chapter  8,  SRS  suppression  using  a  long  period  grating  was  investigated. 
Although  the  results  were  not  as  dramatic  as  those  using  an  intercavity  dopant,  the  long 
period  grating  did  provide  a  significant  increase  in  SRS  threshold.  Using  a  series  of  nine 
gratings  in  a  100  m  fiber,  the  SRS  threshold  was  more  than  doubled.  Additionally, 
relatively  little  effort  would  be  required  to  apply  this  solution  to  lasers  at  different 
wavelengths. 

9.2,  Research  impact 

The  research  described  in  this  dissertation  provides  a  tool  that  could  be  used  in  a 
variety  of  applications  to  increase  the  state  of  the  art  in  laser  brightness.  Several  of  those 
applications  are  examined  in  the  following  discussion. 
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9.2.1.  Laser  Cleanup  and  Combining 

Diode  lasers  are  extremely  useful  for  their  unmatehed  effieieney  and  compact  design, 
but  they  are  limited  in  output  power  and  spatial  quality.  Current  state  of  the  art,  high 
power  diode  lasers  are  capable  of  producing  more  than  8W  [87],  but  the  typical  beam 
quality  of  diode  lasers  is  on  the  order  of  50  times  diffraction  limited  [88].  Although 
individual  devices  have  been  designed  with  high  quality  output  using  novel  cavity 
designs  such  as  the  tapered  laser,  z-laser,  and  alpha-DFB,  these  will  do  little  to  overcome 
the  power  limit  imposed  by  facet  damage.  Arrays  of  these  lasers  can  produce  higher 
power  than  individual  emitters,  but  unless  the  beams  can  be  spatially  combined  and 
cleaned,  the  increased  power  does  not  result  in  a  significant  increase  in  overall  intensity. 
Furthermore,  heat  dissipation  remains  a  limitation  to  the  power  generated  by  diode  arrays 
[88].  Physically  separating  the  arrays  and  combining  the  overall  diode  laser  power  in  a 
fiber  using  stimulated  scattering  can  overcome  this  limitation,  increasing  the  beam 
quality  dramatically  and  increasing  the  overall  power  by  an  order  of  magnitude. 

Although  special  consideration  is  required  to  force  a  diode  laser  to  operate  at  a  single 
frequency  for  use  with  an  SBS  beam  combiner,  narrow  linewidth  devices  have  been 
designed  and  could  easily  be  combined  [89,  90].  A  portion  of  the  research  conducted  in 
the  initial  stages  of  this  dissertation  was  completed  using  a  500  mW  alpha-DFB  laser 
designed  by  SDL,  Inc.  Assuming  careful  optical  engineering  would  result  in  optimal 
coupling  of  these  diffraction  limited  devices  into  a  long  optical  fiber  (perhaps  using  two 
concentric  pump  rings  in  the  arrangement  shown  in  Figure  48),  a  total  of  18  pump  lasers 
could  be  combined  producing  a  diffraction  limited,  spatially  coherent  diode  laser  source 
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with  nearly  the  same  output  power  as  the  highest  reported  diode  laser.  Yet  this  deviee 
would  be  far  from  the  second  order  threshold,  predicted  to  be  more  than  80  W. 

Commercially  available  Diode  bars,  on  the  other  hand,  produce  more  than  50  W  cw 
power  and  have  been  demonstrated  in  the  lab  to  produce  270  W  [88].  With  linewidths  on 
the  order  of  1  THz  [91],  using  SBS  to  spatially  combine  the  individual  emitters  is 
impractical.  However,  SRS  beam  combining  is  a  feasible  alternative  for  sources  over  the 
SRS  threshold  of  nearly  100  W.  As  the  power  per  bar  becomes  closer  to  the  theoretical 
limit  of  10  kW  [88],  the  total  efficiency  of  an  SRS  beam  combiner  becomes  reasonable. 

If  the  diode  laser  bar  had  a  beam  quality  on  the  order  of  50  times  the  diffraction  limit, 
Raman  beam  cleanup  could  increase  beam  brightness  by  a  factor  of  about  2000. 

In  particular,  these  ideas  could  be  applied  to  current  military  requirements. 
Specifically,  an  Air  Force  topic  found  in  the  Department  of  Defense  small  business 
innovation  research  program  solicitation  for  2002  seems  well  suited  for  SBS  beam 
combining  and  cleanup.  Topic  AF02-006  describes  a  need  for  an  efficient  0.9  micron 
laser  for  space-based  lasers.  The  laser  must  produce  1-10  watts  cw  power  and  have  a 
high  spatial  quality  output.  These  characteristics  are  optimally  suited  for  SBS  beam 
combining  using  frequency  stabilized  GaAs  or  AlGaAs  diode  lasers. 

For  other  DOD  applications,  high  power  may  mean  several  kilowatts  average  power. 
The  joint  technology  office  has  funded  research  by  Boeing  to  investigate  high  power 
devices  generated  through  SRS  beam  combining  of  high  power  Nd:YAG  lasers  [92]. 

This  program  was  explored  in  part  due  to  the  initial  results  gathered  and  presented  as  a 
part  of  this  dissertation  research. 
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9.2.2.  SRS  Reduction 


The  development  of  high  power  fiber  lasers  grew  remarkably  from  1997-1999.  Over 
two  years,  the  peak  CW  power  reported  in  those  deviees  rose  by  a  faetor  of  3  [66].  This 
growth  eame  to  a  halt  partly  due  to  the  observation  of  SRS  in  the  laser  output.  The 
applieation  of  an  intercavity  dopant  or  long  period  grating  to  the  problem  of  stimulated 
scattering  in  fiber  lasers  removes  or  at  least  delays  this  limitation.  At  dopant  levels  of 
0.3  molar%,  more  than  an  order  of  magnitude  increase  in  SRS  threshold  is  expected, 
allowing  for  an  increase  of  the  single  mode  fiber  laser  output  from  1 10  W  to  1 .7  kW. 
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Appendix  A.  Fiber  Cladding  Modes 


The  fiber  cladding  mode  dispersion  relation  is  critical  in  the  design  and 
characterization  of  long  period  gratings.  Using  the  notation  of  Erdogan,  the  dispersion 
relation  is  shown  below  [12,  16]. 
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(A.5) 


Ji 

u^Ji  {u^a^) 
Ki'{w^a2) 
w^K,[w^a^) 


Pi  (r)  =  Ji  {u^r^Ni  {u^a^^-Ji  A;  {u^r^ 

qi{r)  =  Ji[u^r)Ni'{u^ai)-Ji'{u^ai)Ni{u^r) 

r,  (r)  =  Ji  '{u2r')Ni  {u^a^^-Ji  A,  '{u^r^ 

$1  (r)  =  Jj  '{u^r^Ni  '{u2ai)-Ji  '(wjaj)  '(^2^) 


(A.6) 


(A.7) 


where  is  the  radius  of  the  fiber  eore  and  is  the  radius  of  the  fiber  eladding,  is 

the  effective  index  of  the  cladding  mode,  and  A.  is  the  wavelength,  n^,  n^,  and  173  are 
the  indices  of  refraction  for  the  core,  cladding,  and  jacket  respectively,  and  Ji,  K,,  and 
Ni  are  the  Bessel  function,  modified  Bessel  function,  and  Neumann  functions 
respectively. 

Once  is  determined  from  the  dispersion  relation,  the  field  distributions  can  be 

calculated.  Since  the  rest  of  this  document  is  only  concerned  with  the  /  =  1  cladding 
modes,  these  are  the  only  modes  explicitly  shown  here.  In  the  fiber  core,  the  r 
component  of  the  electric  field  is 
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Ef  {r,(l),z)  =  J2  {u,r)  +  J^  {u,r)  -  J 

exp  (z^zJ) exp [z  [fiz-  zy^)] 


(A.8) 


The  (j)  -eomponent  of  the  electric  field  is 
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And  the  z  -component  of  the  electric  field  is 
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where  P  is  the  propagation  constant  of  the  mode  and  E^p  is  the  amplitude  of  the  field  in 
that  mode.  The  r  -component  of  the  magnetic  field  is 


Hp‘  (r,  (p,z)  =  eJ  ^{zcTi  [ J2  {uir)  -J^  {u,r)\  +  z^q  \_J^  {u,r)  +  (zZir)]}. 

exp  {ic/))  exp  [z  [Pz-  zyt)] 


.  (A.  11) 


The  (j)  -component  of  the  magnetic  field  is 
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And  the  z  -component  of  the  magnetic  field  is 
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The  fields  in  the  cladding  region  are 
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Finally,  the  fields  in  the  jacket  are 
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From  these  field  distributions,  the  intensity  pattern  of  the  modes  can  be  calculated 
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Appendix  B.  Local  Effective  Brillouin  Gain 


Section  4.2  uses  the  effective  Brillouin  gain  to  explain  the  origin  of  beam  cleanup. 
This  appendix  fills  in  the  steps  between  the  definition  of  the  effective  gain, 

Equation  (4.3),  and  Equation  (4.4).  The  effective  Brillouin  gain  in  a  fiber  is  repeated 
here. 
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The  pump  and  Stokes  electric  fields  can  be  expanded  in  terms  of  the  optical  fiber 
modes  described  in  Chapter  2. 
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where  the  summations  are  over  all  fiber  modes  and  the  m  *  mode  has  a  transverse 
distribution,  y/^  ,  an  amplitude,  (z) ,  a  propagation  constant,  ,  and  a  phase, 

(j)^  .  After  labeling  each  of  these  terms  p  or  s  for  the  pump  or  Stokes  beam,  this 
representation  for  the  intensity  can  be  substituted  into  Equation  (B.l) 
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The  variables  in  Equation  (B.3)  are  defined  by 
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It  is  useful  to  expand  the  summations  in  Equation  (B.3)  as  follows. 

E-  E  +  E  +  E  +  E  + 

m,q,j,v  m,q<m,j,v<j  m,q=m,j,v<j  m,q>mj,v<j  m,q<m,j,v=j 

E  +  E  -  E  -  E  +  E 

m^q=mj,v=j  m,q>mj,v=j  m,q<m,j,v>j  m,q=m,j\v>j  m,q>m,j,v>j 


(B.5) 


When  v  =  j{q  =  m),  Ak  =  0  ( Ak  =  0 )  and  =  0  ( AS  =  0 ),  which  leads  to 


SBeffX^)  p(>y\p  (  1  "  sm,sq,pj,pv  J  /  sm,sq,pj,pv 

\^}^p  \^}  ym,q=mj,v=j 


S  Am,sp,pj,pv{^)ysm,sp,pj,pv^^v\i[{^smp  +  ^  pjv)  ^  +  {^Amp  +  )) 

m^q<mj,v<j 

ZA  (z)y  expl /(AA:  .  z  +  .  )  I  + 

sm,sq,pj,pv  \  //  sm,sq,pj ,pv  y  pjv  '  PJV  J 

=m,j\v<  j 

Z  Am,sp,pj,p.  {^)rsm,sp,pj,pv  ^  +  {^Amp  +  )) 

Z  Am,sp,pj,pv  {^)rsm,sp,pj,pv  ^XP )] 

Z  Am,sp,pj,pv  (^)rsm,sp,pj,pv  ^Xp  [z  (AAt^^^Z  +  A^^^^  )_ 

Z  Am,sp,pj,pv  {^)ysm,sp,pj,pv  ^xp [z  (( +  Ak )  Z  +  (  A^,„^  +  A^,  )) 

Z  Am,sp,pj,pv  {^)ysm,.p,pj,pv  ^XP  [z  (Ak^j^Z  +  A^.^  )] 

^sm,sq,pj\pv  y sm,sq,pj ,pv  ^Xp 


m^q=mj  ^v<  j 


m,q>mj  ,v<  j 


m^q<mj  ,v= j 


m,q>myj\v= j 


m^q<m,j  ^v>  j 


m^q=m,j  ,v>  j 


+ 


+ 


+ 


+ 


+ 


+ 


m^q>mj,v>  j 


[[^Kmq  +  AA:^,v )  ^  +  ( A^.™,  +  ^(l>pp)) 
Two  of  the  nine  summation  terms  in  Equation  (B.6)  are  shown  below. 
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m,q<mj  yV= J 


Am,sq,pj,pv  (^)  are  invariant  when  the  indiees  q  and  m  are  interehanged. 

,  on  the  other  hand,  are  antisymmetrie  and  ehange  sign  when  q  and  m 
are  interehanged.  Therefore,  Equation  (B.8)  can  be  simplified  as  follows 
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Applying  these  same  steps  to  other  pairs  of  sums  in  Equation  (B.6)  results  in 


E  4 


■  sm,sq,pj ,pv 


Z2A  cos(AA:  .  .  W 

sniySq^pj^pv  V  //  sm,sq,pj ,pv  y  pjv  '  PJ'^  J 

m,q=mj,v<  j 

Z  '^Am,sq,p],pv  (  ^  )  Ysm,sq,pj,pv  ^OS  (  Z  +  A^^^^  ]  + 

:.mj,v= j 

Z  ^Am,sq,pj,pv{^)rsm,sq,pj,pv^^y(^k^^^+Akpj^)z  +  [A(I)^^^+A(l)pj^)_  ■ 


m,q<m,j  yV= j 


(B.IO) 


m,q<m,j  ,v<  j 


Z  ^Am,sq,pj,pv{^)rsm,sq,pj,pvCOs\j^Ak^^^  -  Ak  pj^)  Z  +  -  ^^pjv) 


m^q<m,j  ,v<  j 


This  is  the  form  of  the  effective  Brillouin  gain  used  in  Equation  (4.4). 
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